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Создан в целях информирования читательской аудитории о новейших достижениях и перспективах в области 
механики, машиностроения, информатики и вычислительной техники. Издание является форумом для 
сотрудничества российских и иностранных ученых, способствует сближению российского и мирового научно- 
информационного пространства. 


Журнал включен в перечень рецензируемых научных изданий, в котором должны быть опубликованы 
основные научные результаты диссертаций на соискание ученой степени кандидата наук, на соискание 
ученой степени доктора наук (Перечень ВАК) по следующим научным специальностям: 


1.1.7 — Теоретическая механика, динамика машин (технические науки) 

1.1.8 — Механика деформируемого твердого тела (технические, физико-математические науки) 

1.1.9 — Механика жидкости, газа и плазмы (технические науки) 

1.2.2 — Математическое моделирование, численные методы и комплексы программ (технические науки) 

2.3.1 — Системный анализ, управление и обработка информации, статистика (технические науки) 

2.3.3 — Автоматизация и управление технологическими процессами и производствами (технические науки) 

2.3.5 — Математическое и программное обеспечение вычислительных систем, комплексов и компьютерных сетей (технические науки) 
2.3.7 —- Компьютерное моделирование и автоматизация проектирования (технические, физико-математические науки) 
2.3.8 — Информатика и информационные процессы (технические науки) 

2.5.2 — Машиноведение (технические науки) 

2.5.3 — Трение и износ в машинах (технические науки) 

2.5.5 — Технология и оборудование механической и физико-технической обработки (технические науки) 

2.5.6 — Технология машиностроения (технические науки) 

2.5.8 — Сварка, родственные процессы и технологии (технические науки) 

2.5.9 — Методы и приборы контроля и диагностики материалов, изделий, веществ и природной среды (технические науки) 
2.5.10 — Гидравлические машины, вакуумная, компрессорная техника, гидро- и пневмосистемы (технические науки) 
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Анализ напряженно-деформированного состояния цилиндра 
с переменными модулями упругости на основе трехмерных уравнений 


теории упругости 
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Аннотация 

Введение. Функционально-градиентные материалы находят большое применение, т.к. неоднородность свойств 
позволяет управлять прочностью и жесткостью контрукций. Этим вызван большой интерес к данной теме в 
мировой научной литературе. Построение решения таких задач существенно зависит от типа граничных условий. 
В настоящей работе рассматривается равновесие тонкостенного кругового цилиндра, механические свойства 
которого заменяются вдоль радиуса. На цилиндрических поверхностях заданы однородные граничные условия, 
которые до этого не рассматривались, воздействие оказывается на торцах. Математическая постановка задачи 
осуществляется в линейной теории упругости в рамках осесимметричной деформации. В работе построены 
выражения для компонент напряженно-деформированного состояния цилиндра, в которых некоторые 
коэффициенты находятся из решения полученной системы линейных алгебраических уравнений. 

Материалы и методы. Материал цилиндра является линейно упругим, модуль упругости которого линейно 
зависит от радиальной координаты. Основным методом исследования является ассимптотический метод, в 
котором в качестве малого параметра выступает половина логарифма отношения внешнего и внутреннего 
радиусов. Для построения характеристик напряженно-деформированного состояния цилиндра применены 
итерационные процессы. 

Результаты исследования. Для линейно-упругого функционально-градиентного полого тонкостенного 
цилиндра получены однородные решения краевой задачи. Анализ этих решений позволяет раскрыть характер 
напряженно-деформированного состояния в стенке цилиндра. С этой целью проведен асимптотический анализ 
решений, получены соотношения для перемещений и напряжений. Установлено, что эти решения соответствуют 
потраничному слою, при этом их первые члены определяют краевой эффект Сен-Венана, аналогичный теории плит. 
Обсуждение и заключение. Построенное с помощью асимптотического разложения аналитическое решение 
задачи о равновесии неоднородного по радиусу тонкостенного цилиндра может быть использовано для 
численного решения конкретной задачи. Для этого нужно решить полученные системы линейных 
алгебраических уравнений и определить соответствующие коэффициенты. Полученные асимптотические 
представления позволяют анализировать трехмерное напряженно-деформированное состояние. Выбор 
количества членов разложения позволяет рассчитать перемещения и напряжения с заданной степенью точности. 
Этот анализ можег быть полезен при оценке адекватности прикладных методов расчета, применяемых в 
инженерной практике. 


Ключевые слова: линейная теория упругости, функционально-градиентный материал, тонкостенный полый 
цилиндр, однородные решения, пограничный слой, вариационный принцип 
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Введение. Функционально-градиентные материалы находят широкое применение в различных конструкциях. 
Благодаря зависимости механических свойств от координат можно управлять напряженно-деформированным 
состоянием (НДС) деталей. Примером использования такой неоднородности может служить цилиндр, 
механические свойства которого зависят от радиуса. При этом интерес может представлять цилиндр как 
отдельная конструкция, так и являющийся подтелом составного тела, например, соединяющий две среды с 
сильно различающимися свойствами. При расчете НДС тонкостенного цилинда используются некоторые 
прикладные теории. При этом важно знать оценку их адекватности. Особенно в случае неоднородных свойств 
она может быть осуществлена с помощью компьютерного моделирования или асимптотического анализа на 
основе трехмерной постановки. Последним определяется актуальность настоящего исследования. 


Исмайылова Д.Д. Анализ напряженно-деформированного состояния цилиндра с переменными модулями упругости 


Изучению НДС полых цилиндрических тел в рамках линейной теории упругости посвящен ряд исследований. 
В [1, 2] на основе метода сплайн-коллокации и метода конечных элементов изучено механическое поведение 
радиально-неоднородного цилиндра в трехмерной постановке. В [3] изучено НДС цилиндра, свойства которого 
зависят от радиуса, нагруженного равномерным внутреннем давлением. В [4] проведено аналитическое 
исследование для функционально-градиентного пьезоэлектрического цилиндра. В [5] построено точное решение 
для радиально-неоднородного полого цилиндра с экспоненциальным модулем Юнга, постоянным 
коэффициентом Пуассона и степенным модулем Юнга. В [6,7] с помощью метода прямого интегрирования 
получено аналитическое решение осесимметричной задачи термоупругости для сплошного цилиндра, когда 
коэффициент линейного теплового расширения является произвольной функцией от радиуса. В [8] разработана 
общая асимптотическая теория трансверсально-изотропного однородного полого цилиндра. Для трансверсально- 
изотропного однородного цилиндра получены новые группы решений. Приведено сравнение построенных 
решений с решениями, построенными с помощью прикладных методов расчета. В [9, 10] изучены некоторые 
краевые задачи теории упругости для функционально-градиентного изотропного и трансверсально-изотропного 
(плоскость изотропии перпендикулярна оси) цилиндра, в случае, когда модули упругости являются 
произвольными непрерывными функциями от радиуса цилиндра. В [11] проведен анализ задачи изгибной 
деформации для радиально-неоднородного цилиндра. Анализ этих работ показывает, что не для всех типов 
граничных условий на цилиндрических поверхностях имееются асимптотические представления решений. 

В настоящей работе на основе асимптотического анализа трехмерных уравнений теории упругости изучаются 
особенности НДС тонкостенного цилиндра, свойства которого линейно меняются вдоль радиуса. При этом 
внутренняя граница закреплена в осевом направлении и свободна в радиальном направлении. 

Достижение этой цели было получено с помощью нескольких шагов: асимптотическое интегрирование 
дифференциальных уравнений и построение однородных решений; вывод формул для компонент вектора 
перемещений и тензора напряжений; учет краевых условий на торцевых поверхностях. 

Материалы и методы. Рассматривается радиально-неоднороднополый тонкостенный цилиндр 


Г= { Е в фе [0; 2*|, = [-4;% ] в цилиндрической системе координат г,ф, < с началом на его оси. Задача 


о его равновесии, в случае закрепления цилиндрических поверхностей вдоль оси и нулевых нормальных 
напряжений, решается в осесимметричной постановке при действии напряжений на его торцах. 
Краевая задача состоит из уравнений равновесия [8]: 


ых ха +— 2 =0, (1) 
ди 05 г 
05. + 05. + 0, = 0, (2) 
дг 05 г 
где о„,о„,О,б„— компоненты тензора напряжений. 
Определяющие соотношения [8]: 
а. = 6+4) 7% [9 (3) 
дг г 
с. - (26+%) 9% +29 = (4) 
ы 05 д г 
и ди. ди 
ен =(2“+^)— +7 —+— |, р 
т ) г 0 0 © 
с. = с“ + = (6) 
. 0% бк 
Здесь и, =и,(г,5), и, =и, (г, <) — компоненты вектора перемещений. 
Параметры Ламе изменяются по линейному закону вдоль радиуса: 
С(г)=С.т, А(г) =, (7) 
где сС.,Л. — константы. 
После подстановки (3)—(7) в уравнения (1), (2) безразмерная система уравнений принимает вид: 
0?и ди, д?и. ди. 0?и, 
(26, +) +8—2 |+=(б, +2, ее — = +А.Е?е® — + Су ?е2® —_^ — 2СуЕ?и, = 0, (8) 
др? д дрдё 06 08? 
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ди, би. б?и. ди д?и 
С, +—* |+(26, +, )=?| е?® че —^ |+(б,+^, вет —2=0. (9) 
др? др ОЕ? 08 дрдё 
Здесь: 
1 т м 1 Г, ‚. 
р=-Шш| —|, &=— — новые безразмерные координаты; &=-—Ш|-— | — малый параметр в случае 
5 п ”® 2 п 
| А. С. 
тонкостенности; п =./, ре [181], 6 =[-[;1, 1=—°; и, = и, ‚ = а Аве ло ‚ б= ло ; С 
® п п С, С, 


некоторый параметр, имеющий размерность напряжения. 
Рассмотрим задачу, в которой на боковых поверхностях цилиндра заданы однородные граничные условия: 


и: =0, (10) 
би =0. (11) 
К торцам цилиндра приложены напряжения: 
бы =» (В), (12) 
би =, (р), (13) 
(+=62) 
[ [} 


УГ 72 ыы 55 


с = с’ 6 = с’ б= С безразмерные напряжения. 
1 1 1 


Компоненты вектора напряжений 1, (р),6,, (р), (5=1;2) удовлетворяют условиям равновесия. 


Для построения однородных решениий компоненты вектора смещений будем искать в виде: 


и, (р;5) =и(р)е“е, и (р;&) = и(р)е*. (14) 
Подставляя представления (14) в систему (8)—(11), получим: 
(26, +^,) (и" (р) + ги’ (р))+ ве ((б,+7, м (р)+.=и(р))+ =2С, (а?е?> —2)и (р) = 0, (15) 
5. (и"' (р) + вм" (р)) +(2С, +^., )=? (аети (р) + че? (р) +8(С,+^. )аети’ (р) = 0, (16) 
У. =0, (17) 


= 0. (18) 


р=-+1 


[(26, + м (р)-+=^. (и(р)+ чету (р)) | 


Исследуем краевые задачи (15)-(18) при &—0. Для решения (15)—(18) при &-—0 воспользуемся 
асимптотическом методом [9-13]. 
Ненулевые решения (15)—(18) соответствуют третьему итерационному процессу, компоненты вектора 
смещений будем искать в виде разложений по малому параметру: 
и(3) (р) = =(и (р) + 5И.1 (р) я .), 


из (р) ==(и (р)-+ =, (р)+---), (19) 
а = =! (В +В, +---). 


После подстановки разложений (19) в уравнения (15)—(18) для членов первого порядка имеем: 


(26, +2 (р)+Вь (6, +2) (в) бы (р) =0, со) 
Сом (р) + В, (С, + изо (рР)+ (26, +0 Ви (р) =0, (20 
У (Ри. =0, (22) 


= 0. (23) 


р=+1 


(26+ )мь (р) + Вим (р) 


Следуя [13], спектральная задача (20)—(23) соответствует потенциальному решению для плиты. 
Таким образом решения представляются в виде: 


а) ибо (р; &) = =>. т, (-26. В, эпВих $ (Вор) = 0(=)) й оч (Ви +В, +.- .) :] ь (24) 
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> : 1 
м" (р; 5) = => (2.В, зшВь, со$ (Вир) + 0(=)) * вар (ры +аВ,, +.- .) :) . (25) 
Здесь В.„ является решением уравнения: 
соз Вх = 0. (26) 
Напряжения, соответствующие решениям (24) и (25), имеют вид: 
ов = т (- —4С2Вз, зтВих с0$(Вир) + О (Е (одуееме[ Вох + ЕВь +8) , (27) 
0. = => Т, (4 ?В, п Во т (Вир) + О(= х еек Вик +=Вь, +-- 3) ь (28) 
& 
е 1 
бе 9 = =>, (4 СоВа 91п Вох соь(Выр)-+0())екр] (ры +В, +=) (29) 
5 = 0(:). (30) 
[9] 1 
Ь) из (р;Е) = =>. :(2С.В:, созВ., соз (Вир) +О(= =))х ор +В, +.. 3} (31) 
Э Ь : 1 
иг” (р; 5) — => Е (26,88, созВь, п (Вир) ы: 0(=)) 2 ее: (Во +=В,, +-- .) :] . (32) 
Здесь Во, является решением уравнения: 
зшВь; = 0. (33) 
Напряжения, соответствующие решениям (31) и (32), имеют вид: 
[2] . т 
0%? =>. Е (-403В% созВо (Вир) + ое +аВ, +. 3} (34) 
со 1 
6%? = =>Е ‚ (4С2В, созВо, соз (Вир) +О (= ево + В, +] (35) 
Ьы : 1 
од = 52 (4010 воз, м (Вьр)+ (6) вы ару += в6) 
68? =0(=). (37) 


Общее решение (15)—(18) будет суперпозицией решений (24), (25), (31), (32): 
и, (р;5) = =>. Т, (-26.В&, зпВо, т (Вир) + 0(=))х 


С] 


ке] (ры +В» ++ =>. Е :(2С,В&, созВ; соз(Выр) + обе] (р + &В,, + =) (38) 


и, (р;5) — =>, (26,В, зп Ви соз(Вир) + 0(=))х 


кем (Вы + ёВ, + ...) в (26,В3, созВ, чт (Вор) +О(=)х окр + В, + ...) :] (39) 


Решения (24), (25), (31), (32) имеют характер пограничного слоя. При удалении от торцов решения (24), (25), 
(31), (32) экспоненциально убывают. 
Чтобы определить константы Т., ЕЁ, воспользуемся вариационным принципом Лагранжа. Вариационный 


7 


принцип принимает вид [8]: 


21 
> ] [(с„. 1, ) би, + (5 —ь,, би. | е2 р = 0. (40) 
— = 
Подставляя (24—36) в (40), имеем: 
Х Мать = Ро, 41 
К=1 
Зо, в) о 
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Здесь: 


М „ =16б0 ВВ (Во, —Ва» ) зп Вих Во, чп (Во, — В» )х =з- А беты › (при а к) 
‘ = = 


, 28, 1 ов, 
М ; =16С3В8, п? Ву, она (при =К) 


| ы 
Ро = 2СВ5, яп Во; | (> (р)соз (Выр)- и бррр уар-ече е } 


+1 (6. (р)есз (Вир) и. буть очен 


г ог +Во; )1 + Во; |1 
О; =166. В Ви (Во В, ) созВь, ; созВ,, чл (Во; — В» ) эВ ыы ‚ (при #* 1) 
, - | 
О; =16С5Вь, со? Во, |= т ны (при 1=] ) 
| = & 


(ею (Выр) + бумер }ар-ено м, 
7. 


Постоянные Т,,Ё, (р=1,2,...) находятся из систем линейных алгебраических уравнений (41), (42), 


Ру = 2С.Вь; созВь, | 


Во 


7 
= 


+ (6. (р)соз(Вор) + (р)з (Выр))ар . 5 


Т, =То +ЕТи+.-*, 
Е =ЕВ +ЕЁ, +... 


аналогичные которым изучены в [13]. 

Результаты исследований. В работе в осесимметричной постановке рассмотрено решение задачи линейной 
теории упругости для функционально-градиентного полого тонкостенного цилиндра, свойства которого 
изменяются по толщине по линейному закону. На боковых поверхностях цилиндра заданы однородные 
перекрестные граничные условия, на торцах задан вектор напряжений. Построенные однородные решения 
удовлетворяют граничным условиям на цилиндрических поверхностях. Для их построения использован 
асимптотический подход на основе раздожения по малому параметру, характеризующему относительную 
толщину цилиндра. Для учета неоднородных граничных условий на торцах получены системы линейных 
алгебраических уравнений, аналогичные изученым в литературе. Показано, что построенные решения НДС 
имеют погранслойный характер, который соответствует краевому эффекту, аналогичному теории неоднородных 
плит, носящему имя Сен-Венана. 

Обсуждение и заключение. Обычно при изучении НДС тонкостенных конструкций строятся прикладные 
методы расчета, снижающие размерность задачи. В этой связи актуальной является задача определения 
диапазона геометрических и механических параметров, в которых эти методы дают приемлемую точность. 
Построенные в работе на основе асимптотического анализа решения трехмерных уравнений позволяют оценить 
адекватность таких прикладных теорий с заданным наперед порогом точности. Кроме того, эти решения могут 
найти применение при оценке численных решений задач для конструкций с функционально-градиентными 
материалами. 
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Реализация базовых операций для разреженных матриц при решении 
обобщенной задачи на собственные значения в комплексе АСЕГАМ-СОМРОЗ$ 
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Аннотация 

Введение. Широкое использование пьезоматериалов в различных отраслях стимулирует изучение их физических 
характеристик и обусловливает актуальность таких изысканий. В рассматриваемом случае модальный анализ 
позволяет определить рабочую частоту и коэффициент электромеханической связи пьезоэлементов различных 
устройств. Эти индикаторы представляют серьезный теоретический и прикладной интерес. Цель исследования — 
разработка численных методов для решения задачи определения частот резонанса в системе упругих тел. Для 
достижения цели нужны новые подходы к дискретизации задачи на основе метода конечных элементов и 
выполнение программной реализации выбранного метода на языке С# на платформе „пе. Актуальные решения 
созданы в контексте библиотеки классов комплекса АСЕГАМ-СОМРО$З. Основанные на обращении матриц 
известные методы решения обобщенной задачи на собственные значения неприменимы к матрицам большой 
размерности. Для преодоления этого ограничения в представленной научной работе реализована логика 
построения матриц масс и созданы программные интерфейсы для обмена данными о задачах на собственные 
значения с модулями пре- и постпроцессинга. 

Материалы и методы. Для реализации численных методов задействовали платформу .пеё и язык 
программирования СЯ. Валидация результатов исследования проводилась путем сравнения найденных значений 
с решениями, полученными в известных САЕ-пакетах (англ.  сотрщег-аа4е4 епошеешпе — 
компьютеризированная инженерия). Созданные подпрограммы оценивались с точки зрения производительности 
и применимости для задач большой размерности. Проводились численные эксперименты с целью валидации 
новых алгоритмов в задачах малой размерности, которые решаются известными методами в МАТЕАВ. Далее 
подход тестировали на задачах с большим числом неизвестных и с учетом распараллеливания отдельных 
операций. Чтобы избежать нахождения обратной матрицы, программно реализовали модифицированный метод 
Ланцоша. Рассмотрели форматы хранения матриц в оперативной памяти: триплеты, СВ, С$С, ЗКуПпе. Для 
решения системы линейных алгебраических уравнений (СЛАУ) задействовали итерационный симметричный 
метод ГО, адаптированный к этим форматам хранения. 

Результаты исследования. Разработаны новые расчетные модули, интегрированные в библиотеку классов 
комплекса АСЕГАМ-СОМРО$. Проведены расчеты для определения применимости различных форматов 
хранения разреженных матриц в оперативной памяти и различных методов реализации операций с разреженными 
матрицами. Графически визуализирована структура матриц жесткости, построенных для одной и той же задачи, 
но с различной перенумерацией узлов конечноэлементной сетки. Применительно к задаче теории 
электроупругости обобщены и представлены в виде таблицы данные о времени, необходимом на выполнение 
базовых операций с матрицами жесткости в различных форматах хранения. Установлено, что перенумерация 
узлов сетки дает существенный прирост производительности даже без изменения внутренней структуры матрицы 
в памяти. С учетом поставленных задач исследования названы преимущества и слабые стороны известных 
форматов хранения матриц. Так, СЗК оптимален при умножении матрицы на вектор, ЗК $ — при обращении 
матрицы. В задачах с числом неизвестных порядка 103 выигрывают в скорости итерационные методы решения 
обобщенной задачи на собственные значения. Оценивалась производительность программной реализации метода 
Ланцоша. Измерялся вклад всех операций в общее время решения. Выяснилось, что операция решения СЛАУ 
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занимает до 95 % от общего времени работы алгоритма. При решении СЛАУ симметричным методом ГО 
наибольшие вычислительные затраты нужны для умножения матрицы на вектор. Для увеличения 
производительности алгоритма прибегли к распараллеливанию с общей памятью. При использовании восьми 
потоков производительность выросла на 40-50 %. 

Обсуждение и заключение. Полученные в рамках научной работы программные модули были внедрены в пакет 
АСЕГАМ-СОМРОЗ. Оценена их производительность для модельных задач с квазирегулярными 
конечноэлементными сетками. С учетом особенностей структур матриц жесткости и масс, получаемых при 
решении обобщенной задачи на собственные значения для электроупругого тела, определены предпочтительные 
методы для их обработки. 


Ключевые слова: пьезоматериалы, метод конечных элементов, разреженные матрицы, обобщенная задача на 
собственные значения, метод Ланцоша, подпространство Крылова, модуль препроцессинга, модуль 
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Введение. Устройства из пьезоматериалов широко используются, давно активно изучаются и 
совершенствуются. Отдельно следует отметить медицинские ультразвуковые приборы (оборудование для 
диагностики, ультразвуковые скальпели) [1-4] и мобильные генераторы энергии [5]. Исследование [6] описывает 
комбинирование фото- и пьезоэлектрического эффектов для создания действенных компактных источников 
энергии. В науке и промышленности изучаются новые материалы, рассчитанные на эксплуатацию в 
специфических условиях. В [7] рассматривается создание бессвинцового пьезоактивного состава, подходящего 
для эксплуатации при различных температурах. 

В исследовании пьезоэлементов важную роль играет этап модального анализа, позволяющий установить 
частоты резонанса и антирезонанса устройства. Эти данные: 

— необходимы для выяснения рабочей частоты устройства; 

— позволяют определить коэффициент электромеханической связи — важный индикатор эффективности 
устройства; 

— являются входной информацией в численных экспериментах для задач на вынужденные колебания. 

Цель исследования — создание численных методов решения задачи определения частот резонанса для 
системы упругих тел. Достижение заявленной цели требует решения двух задач. Первая: разработать методы 
дискретизации задачи на основе метода конечных элементов (МКЭ). Вторая: провести программную реализацию 
выбранного метода на языке С# на платформе „пе. Все известные программы учитывают контекст библиотеки 
классов комплекса АСЕГАМ-СОМРОЗ [8]. При решении обобщенной задачи на собственные значения широко 
применяются методы, основанные на обращении матриц. Однако они неприменимы к матрицам большой 
размерности. В представленной научной работе это ограничение преодолевается следующим образом: 

— дополнительно реализована логика построения матриц масс; 

— созданы программные интерфейсы для обмена данными о задачах на собственные значения с модулями пре- 
и постпроцессинга. 

Материалы и методы. В первую очередь предлагаемый подход призван решать статические задачи 
электроупругости при реализации метода осреднения [9], который задействуют для расчета эффективных 
свойств пьезокомпозитов. В связи с этим на этапе построения глобальных матриц МКЭ представлены только 
матрицы жесткости. В данном исследовании дополнительно реализовали логику построения матриц масс и 
разработали программные интерфейсы (аррИсаНоп ргоэтатите пиегасе, АР1, англ. — интерфейс прикладного 
программирования) для обмена данными о задачах на собственные значения с модулями пре- и постпроцессинга. 
Разработанные подпрограммы оценивались с точки зрения производительности и применимости для задач 
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большой размерности. Проводились численные эксперименты с целью валидации созданных алгоритмов для 
таких задач малой размерности, которые позволяют получить решение общими методами в вычислительном 
пакете МАТГАВ. Далее выполнялось тестирование на задачах с большим числом неизвестных и с учетом 
распараллеливания отдельных операций. 

Математическая модель решаемой задачи состоит из определяющих соотношений [9]: 


Р‚кб?и + азроюи -У:с0=ф, У:В=0, (1) 
с=с®.(Е+Вё)-е7.Е, О+е,О=е,: (&+6,8)+57-Е, (2) 
= = (Уи+Уи!)/2, Е =-УФ. (3) 


Здесь со — тензор напряжений; р; — плотность тела; = — тензор деформаций; и — вектор перемещений; р — 
вектор электрической индукции; Е — вектор напряженности электрического поля; /, — вектор массовых сил; 


ф — электрический потенциал; о, В.;, с, — коэффициенты демпфирования; с®, е?т, э$ — тензоры упругих 


констант, пьезомодулей и диэлектрических проницаемостей; индекс ] —Й номер тела в модели. 


Дискретизация выполняется заменой: 
и(х,г) = № (х)- 0, Фр = М(®-ФО. 


Здесь №, — матрица функций формы для поля перемещений; № — вектор функций формы для электрического 
потенциала; И (1), Ф(!) — глобальные векторы соответствующих узловых степеней свободы. 


В таком случае исходная задача (1-3) приобретает вид: 
М.-а+К-а=Е. (4) 
Здесь матрицы М и Кявляются глобальными матрицами масс и жесткости соответственно, а вектор а 
представляет собой общий вектор неизвестных: 
а= [(, Ф]. 
В задаче теории электроупругости: 


= Ми 0 К = Кии Ка 
— № о) (# -Кы/ 0 


Матрицы Мш› Кши Къ — симметричные. В случае гармонических колебаний на собственной частоте в; 


можно записать: 
а=у, 91(©и), 
обозначив через у, соответствующий собственный вектор. 
Рассмотрим свободные колебания если Г =0. В этом случае задача (4) представляется в виде: 
—0?Му +К-м, =0. (6) 

Таким образом, исходная задача сводится к обобщенной задаче на собственные значения (6). Для ненулевого 
у; неравенство (6) решается нахождением матрицы, обратной К. Однако при этом разреженная матрица 
становится заполненной, то есть метод непригоден для матриц больших размеров. Поэтому нужно использовать 
другие методы, не требующие нахождения обратной матрицы. Для решения этой задачи в данной работе 
программно реализован модифицированный метод Ланцоша [10]. Автор этой модификации — Т.С. Мартынова. 
Описание разработки в данной статье не приводится. Из используемых в методе операций наиболее затратной с 
точки зрения вычислительных ресурсов оказалось решение системы линейных алгебраических уравнений 
(СЛАУ), необходимое для выполнения спектральной трансформации. 

Матрицы М и К — разреженные, с небольшим числом ненулевых элементов. Для хранения таких матриц 
в оперативной памяти используются несколько форматов: 

— триплеты или координатный формат; 

— С5В (англ. сотргеззе4 зрагзе го\у — сжатый разреженный ряд); 

— С$С (англ. сотргеззе4 зрагзе сои — сжатый разреженный столбец); 

— формат хранилища ЗКуПпе (англ. ЗК $ — метод). 

Координатный формат предполагает хранение троек (триплетов) значений (1, }, К), представляющих собой 
координаты (1, ]) и значения (К) ненулевых элементов. СЗК иногда называют СВ$ или Йельским форматом. Он 
предполагает хранение разреженной матрицы в виде трех массивов. Рассмотрим матрицу размера № с № 
ненулевыми элементами. Опишем возможную организацию ее хранения. Все ненулевые элементы нужно 
разместить в одном массиве размера № . Позиции этих элементов в столбцах разместить в другом массиве 
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размера №7 ‚ а третий массив размера № использовать для хранения индексов первых элементов строк. 
Аналогично реализуется хранение в формате С$С. 

Формат $К$ предполагает хранение ленты матрицы переменной ширины, включающей в себя все ненулевые 
элементы. В этом случае допускается хранение нулей. Эффективность этого формата зависит от перенумерации 
строк матрицы. Методы сокращения размера ленты описаны в [11], однако требует отдельного исследования их 
применимость к матрице жесткости, получаемой при решении трехмерной задачи с использованием МКЭ. 

Для решения СЛАУ задействовали итерационный симметричный метод ГО (Зуттешле ГО Меоч, 
ЗУММГО [12]), адаптированный к перечисленным выше форматам хранения. 

Результаты исследования. В начале исследования выбрали оптимальный формат хранения для разреженных 
матриц. Координатный формат позволяет быстро добавлять и изменять элемент матрицы. Эти операции 
необходимы на этапе сборки глобальной матрицы и при учете граничных условий. Кроме того, для плохо 
обусловленных матриц, к которым относится К, часто применяют предварительное преобразование для 
нормирования. Его также удобно выполнять в координатном формате. Однако такой формат неэффективен, если 
речь идет об алгебраических операциях. 

СУК плохо приспособлен для изменения структуры матрицы: добавляя ненулевой элемент, нужно выполнять 
вставку в два массива. При этом матрица умножается на вектор очень просто и эффективно. 

ЗК$ имеет аналогичные проблемы с добавлением ненулевых элементов и сильно зависит от перенумерации 
неизвестных в задаче. Остановимся на примере квазирегулярной сетки, которая используется в пакете АСЕГАМ- 
СОМРО$З для работы с представительными объемами композитов. Ширина ленты, содержащей все ненулевые 
элементы, может быть определена заранее и зависит от числа узлов и типа конечного элемента. В общем случае 
произвольной конечноэлементной сетки сложно заранее оценить размер ленты. 

В численных экспериментах использовали четыре способа нумерации неизвестных. На рис. 1 представлена 
структура матриц жесткости, построенных для одной и той же задачи, но с различной перенумерацией узлов 
конечноэлементной сетки. 


г) 


Рис. 1. Структура матриц жесткости при различных методах нумерации узлов: а — неизвестные упорядочены по узлам; 


б — сначала упорядочены по узлам перемещения, затем потенциалы; в — узлы отсортированы по слоям КЭ-сетки, а 
неизвестные — по узлам; 6, г — узлы отсортированы по слоям КЭ-сетки, а неизвестные — как в примере 
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Итак, сетка представляла собой куб с регулярным разбиением восьмиузловыми конечными элементами. Для 
иллюстраций использовалась модельная матрица из 500 строк. Матрицы, представленные на рис. 1 а и1б, не 
подвергались дополнительной перенумерации узлов и отличаются только нумерацией степеней свободы. В Та: 

а = [Ми М.›Ф!»---›Инх>Ику >И, Фи]. 

В | бнеизвестные, отвечающие за распределение потенциала, собраны в конце вектора: 

а=[и м, И, ,,-.> Им, Иму Им ›Фрь-- Фи]. 

Неизвестные в матрицах !1в и 12 занумерованы аналогично, но узлы конечноэлементной сетки 
предварительно перенумерованы согласно их координатам путем поочередной сортировки всех узлов по каждой 
из координат. Такой прием широко используется для построения более эффективных модулей решения СЛАУ, 
так как позволяет работать с матрицей в специальном ленточном формате, удобном для распараллеливания. Для 
комплекса АСЕГАМ-СОМРО$ реализованы подобные внешние модули [13—15], однако в данной работе 
использовались только форматы хранения разреженных матриц общего вида. 

В таблице 1 сведены данные о времени, необходимом на выполнение базовых операций с матрицами в 
различных форматах. 


Таблица 1 
Время для выполнения базовых операций с матрицей жесткости в задаче теории электроупругости. 
19 652 строки 
Формат Затраченное время, мс 
Операция 

хранения Та 16 1 в 12 

С5В Конвертация из координатного формата 123 132 97 117 

С5В Умножение на вектор, 100 операций 260 260 260 260 

5К5 Конвертация из координатного формата 690 703 124 268 

5К$ Умножение на вектор, 100 операций 60 558 | 61450 7 616 22 113 


Результаты экспериментов показали, что операция преобразования из координатного формата хранения в 
компактные занимает мало времени. При этом перенумерация узлов сетки для формирования блочно-ленточной 
матрицы позволяет получить заметный прирост производительности даже без изменения внутренней структуры 
матрицы в памяти. Формат СЗВ оказался оптимальным с точки зрения эффективности операции умножения 
матрицы на вектор. При обращении матрицы формат 5К$ более эффективен, однако для задач с числом 
неизвестных порядка 103 заметно быстрее работают итерационные методы решения обобщенной задачи на 
собственные значения. 

Далее экспериментально оценили производительность программной реализации метода Ланцоша. Измерили 
вклад всех операций в общее время решения. В результате установили, что операция решения СЛАУ занимает 
до 95 % от общего времени работы алгоритма. В ходе работы алгоритма строится подпространство Крылова, и в 
зависимости от его размерности меняется число СЛАУ, которые необходимо решать. Отметим, что размерность 
подпространства Крылова выбиралась на основе эвристик относительно числа искомых собственных значений. 
При этом СЛАУ отличаются только правой частью, благодаря чему требования к выделяемой памяти остаются 
невысокими. Среди базовых операций, применяемых в ходе решения СЛАУ методом ЗУММГО, наибольшие 
вычислительные затраты нужны для умножения матрицы на вектор. 

Для увеличения производительности алгоритма реализовано простейшее распараллеливание с общей 
памятью. Для формата СВ$ выделялись блоки строк. Они передавались в отдельные потоки, которые вычисляли 
соответствующие компоненты результирующего вектора. Прирост производительности составил 40-50 % при 
использовании 8 потоков. При этом для матриц порядка 103 элементов прирост был около 40 %, а для матриц 
порядка 10% — около 50 %. 

Обсуждение и заключение. В рамках данной работы реализован метод решения обобщенной задачи на 
собственные значения для матриц, получаемых при моделировании электроупругих тел. Созданы программные 
модули на языке С# для построения матриц масс методом конечных элементов и выполнения вспомогательных 
операций в рамках метода Ланцоша (работа с векторами подпространства Крылова, переортогонализация, 
нахождение собственных векторов). Вычислительная сложность обусловлена в основном операциями 
умножения разреженных матриц на вектор. В связи с этим проводились численные эксперименты по 
определению оптимальных форматов хранения матриц, оптимальной структуры матрицы, получаемой в 
результате перенумерации узлов КЭ-сетки и степеней свободы в узлах. Разработана версия итерационного 
алгоритма ЗУММГО с использованием параллельных вычислений. Итоговая схема работы включает три пункта. 


Оганесян П.А. и др. Реализация базовых операций для разреженных матриц 


Во-первых, строятся глобальные матрицы в координатном формате с алгоритмом перенумерации (рис. 1 в). Во- 
вторых, данные преобразуются в формат СЁ$. В-третьих, данные обрабатываются методом Ланцоша, который 
включает метод ЗУММГО для решения СЛАУ. Результаты работы включили в программный пакет АСЕГАМ- 
СОМРОЗ. 
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Аннотация 

Введение. Формирование параметров качества поверхностного слоя и эксплуатационных свойств деталей 
происходит на протяжении всех этапов их изготовления. Однако решающее влияние чаще всего оказывают этапы 
финишной обработки. Поэтому в современном цифровом машиностроении задача технологического 
обеспечения высокого качества поверхностного слоя детали является одной из важнейших при решении 
проблемы повышения качества, надежности и увеличения жизненного цикла производимых машин. Ведущую 
роль в повышении эксплуатационных характеристик деталей машин играет обработка поверхностным 
пластическим деформированием, сущность которой заключается в том, что требуемые параметры качества 
деталей достигаются не удалением слоя материала, а его пластическим деформированием. В процессе обработки 
производится изменение как размеров деталей, так и физико-механических характеристик поверхностных слоев, 
управляя которыми технолог имеет возможность значительно увеличивать жизненный цикл производимой 
продукции. Целью настоящих исследований является обеспечение необходимых параметров качества 
поверхностного слоя при обработке эксцентриковым упрочнителем. 

Материалы и методы. В статье представлены результаты исследований нового метода обработки 
поверхностным пластическим деформированием — осциллирующим эксцентриковым упрочнителем. 
Рассматриваемый метод обработки позволяет получать высокое качество обработанной поверхности, 
осуществлять обработку крупногабаритных деталей в местах, являющихся концентраторами напряжений, 
обрабатывать сварные швы, небольшие участки поверхностей, упрочнение которых необходимо для выполнения 
деталью своего служебного назначения. Выполнен комплекс теоретических исследований, по результатам 
которых определены параметры единичного взаимодействия индентора с поверхностью детали, диаметр 
пластического отпечатка и его глубина. 

Результаты исследования. Получены зависимости для определения шероховатости поверхности, глубины 
упрочненного слоя и степени деформации. Полученные формулы прошли проверку адекватности 
экспериментальными исследованиями. 

Обсуждение и заключение. Полученные результаты исследований могут быть использованы при 
технологическом проектировании процессов обработки поверхностным пластическим деформированием. 
Определены дальнейшие задачи по исследованию рассматриваемого метода обработки. 


Ключевые слова: осциллирующий инструмент, эксцентриковый упрочнитель, шероховатость поверхности, 
глубина упрочненного слоя, степень деформации 


Благодарности: авторы выражают благодарность редакции журнала и рецензентам за внимательное отношение 
к статье. 


Для цитирования. Тамаркин М.А., Тамаркин Э.Э., Хашаш О.С.А. Формирование качества поверхностного слоя 
при отделочно-упрочняющей обработке деталей эксцентриковым упрочнителем. А4уапсеа Епотеептив Кезеагсй 
(Коя1оу-оп-Ооп). 2023;23(2):130-139. № рз://Ч01.0го/10.23947/2687-1653-2023-23-2-130-139 


© Тамаркин М.А., Тищенко Э.Э., Хашаш О.С.А., 2023 


'А4уапсеа Епятеениз Кезеагсй (Козоу-оп-Ооп). 2023;23(2):130-139. е155М 2687—1653 


Опвта[ агие 
бигасе ОпаШу Еогишо опдег Раг5 Епи$В шо ап 59`епоепто ТгеабпепЕ ув ап Ессепылс Нагдепег 


МИКВай А. Ташанст, ЕНпа Е. Твспепко Сл, Отаг СА НазвазВ © 


оп 5{ае Тесбтиса! Ошуегзну, Возюу-оп-Ооп, ВКяз$1ап Недеганоп 
4 Ппа_изпепКо @ тай.га 


АЪ5етгасЕ 

Гитодисйоп. Те Ююглайоп оЁ Фе диаШу рагатеег$ оЁ Ше за! Расе 1ауег ап Фе орегайопа! ргорегиез оЁ Ше рап5 
оссиг$ ШгоиевоиЕ а| ${асез оЁР Феи тапагасвте. Но\жеуег, пе 4еслууе ппрасе 15 105 ойеп ехецеа Бу Ше з(асез оЁ 
Визы ло. ТБегеРоге, 11 подег 41а] епопеенто, Ше {а5К оф ргосезз заррогЕ оЁ ШеВ адиаПеу оЁ Фе зи!Ёасе 1ауег оЁ Ше 
раке 15 опе оЁ Фе сваПепзез ш з9]уше фе ргоет оЁ паргоушс Фе диаШу апа генабиу ап4 шсгеазшо Фе Не сусе 
оЁ тапиРасвиге4 тасЫшез. ЗигЁасе р!азис ЧеогтаНоп иеацтепе 13 шзиитета! ш пиргоуше Фе регоппапсе 
свпагасвег1$Ис$ оЁ тасНше раг. [5 ез5епсе 15 фаЕ Ше гедие4 даа у рагатеег$ оЁ ра$ аге оаше4 по! Бу гетоуше 
а |ауег о таепа1, Ба Бу р1азис аеогтайоп. Рийпе Ше ргосеззте, Бой Ше 4птеп$101$ оЁ Ше раг$ ап4 Фе рБуз1са1 
ап тесБашса| ргорегиез оЁ Фе зигЁасе 1ауегз аге сВапзе4. ш 11$ сазе, Ше {есБпо|0515Е Ваз Ше оррогипйу ю 
яешЙсапНу шсгеазе Ше Пе сусе оЁ Ше тапа{асвгед ргодис($ @гои?В сопгоШпе е ргосез$. ТВезе за е$ аге апте4 
а ргоу14т» Ше гедите4 диаШу рагатеегз оЁ Пе за! Расе 1ауег ипдег ргосеззтя УИ ап ессепелс Баг4епег. 

Маепай апа Мейоа$. Те агасе ргезепз Фе гези!5$ оЁ гезеагсЬ оп а пем тефо4 оЁ зигРасе р1азНс аеГогтаноп 
(цеантепе — \1 ап озсШаНи? ессепёс Багдепег. ТВе сопз14еге4 ргосеззтх тео епа]ез ю оБашт 6126 ацаШу оЁ 
Фе пезже4 зигЁасе, ю ргосез$ 1агое-$17е4 раг$ 11 р1асез ШФа{ аге зтезз сопсешгаюгз, 10 ргосез$ \е!4$, зта| агеа$ оЁ 
зиГасез, \уВозе Ваг4ептас 15 пееде4 Гог Ше рагё ю РП 1$ пцепдеа зегутсе. А зеё оЁ Пеотенса| зв 1ез \аз сагле4 оит; 
Феш геза5 ргоу1Аед д&егииите Фе рагатеегз оЁР а зше ицегасноп оЁ ше ш4дещег ап4 Фе зигЁасе оЁ Фе ра, Ше 
Фате(ег оЁ Фе разйс пари: апа 15 аер®. 

Кези5. Переп4епсез ог 4еегийите Ше зигРасе гоиеВпез$, фе Чер оЁ Фе Пагдепе4 1ауег ап@ Фе Честее оЁ 
еРогтаНоп \еге оашеч. ТВе гезШипз Рогии1аз \еге (е$е4 Гог адедиасу Бу ехрегипепа1 зи 1ез. 

Оёсизяюпт апа Сопсшяоп. ТБе ощате4 гезеагсВ гезиз сап Бе изе4 шт фе {есВпо]оз1са| 4ез1еп оЁ за асе разйс 
еРогтаНоп ШеайтепЕ ргосеззез. Еиг ег {азКз ог Ше зш@у оЁ Ше соп$14еге4 ргосеззтз плефо4 аге деегилтеч. 


Кеумога$: озсШайпе 1001, ессепёлс Ваг4епег, загРасе гоп пезз, Ваг4епе4 Тауег дер, 4еотее оЁ 4еогтаНоп 


Аскпо\еазетет 6: Фе апфог$ уошША ПКе 10 ФапК Фе е4йопа1 Боага оЁ Фе оигпа! ап Фе геуле\уег$ ог Фет айепнуе 
аиши4е 0 Ше агисе. 


Ког сНайбоп. Татагкт МА, ТизЬсБепко ЕЁ, Назразь ОСА. Зи асе ОпаШу Еогпуше ип4ег Райз ЕиизЫт® апа 
биепоетше ‘Тгеайптепё \Ш ап Ессепилс Наг4епег. Айуапсеа Епвтеетив Кезеатсй (Ко5оу-оп-Ооп). 
2023;23(2):130-139. Врз://401.0г$/10.23947/2687-1653-2023-23-2-130-139 


Введение. Надежность и долговечность деталей машин в значительной степени зависит от качества их 
поверхностного слоя. Из многочисленных работ по технологии машиностроения известно, что формирование 
параметров качества поверхностного слоя происходит на всех стадиях их изготовления. Однако решающее 
влияние чаще всего оказывают этапы финишной обработки. Поэтому в современном цифровом машиностроении 
уделяется повышенное внимание технологическому проектированию высокоэффективных операций финишной 
обработки деталей, что позволяет решить сложнейшую задачу увеличения их жизненного цикла. Ведущую роль 
в повышении эксплуатационных характеристик деталей машин, осуществляемом на финишных операциях, 
играет обработка поверхностным пластическим деформированием (ППД). В отличие от традиционных методов 
обработки резанием, параметры качества поверхностного слоя при ПИД достигаются путем осуществления 
пластической деформации специальными инструментами или рабочими средами. В процессе обработки 
одновременно с изменением размера обрабатываемых деталей производится изменение физико-механических 
характеристик поверхностных слоев, управляя которым технолог имеет возможность значительно увеличивать 
жизненный цикл производимой продукции. 

Следует отметить, что широкому распространению многих методов ППД в промышленности препятствует 
слабая изученность их основных закономерностей, трудности, возникающие в процессе проектирования 
оптимальных сочетаний режимов обработки и конструктивных параметров средств технологического 
оснащения. Во многих случаях режимы обработки назначаются исходя из результатов частных 
экспериментальных исследований, что приводит к низкой эффективности обработки [1-6]. 

Целью настоящих исследований является обеспечение необходимых параметров качества поверхностного 
слоя при обработке эксцентриковым упрочнителем. 
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Материалы и методы. Необходимость применения ППД в условиях современных машиностроительных 
производство приводит к созданию новых методов обработки. Одним из таких методов является обработка ППД 
осциллирующим инструментом — эксцентриковым упрочнителем. 

На рис. 1 представлена кинематическая схема эксцентрикового упрочнителя, состоящего из вибрирующего 
корпуса 1, подвешенного на плоских пружинах 2. Колебания вибрирующего корпуса 1, действующие нормально 
к обрабатываемой поверхности, возбуждаются вращением эксцентриковой массы (дебаланса) 3 вокруг 
вертикальной оси. Ось вращения эксцентриковой массы ограничена от перемещений относительно 
вибрирующего корпуса 1. Вращательное движение эксцентрику передаётся от электродвигателя 5 через гибкий 
вал 6. К корпусу 1 прикреплена инструментальная головка 4 с инструментом соответствующей геометрической 
формы. Движение инструмента 4 ограничивается ограничителем 7 (обрабатываемая деталь). При этом 
инструмент представляет собой индентор со сферической заточкой или может выполняться в виде ролика или 
шарика. Вибрационную систему в эксцентриковых упрочнителях можно представить как одномассовую систему 
с двумя степенями свободы, находящуюся под действием силы, изменяющейся по гармоническому закону. Для 
исследования динамики системы рассмотрим особенности её свободного колебания под действием 
центробежного вибровозбуждения и характер движения системы, ударяющейся об ограничитель (деталь). 

При свободном колебании вибрирующая система, закрепленная на конце плоских пружин 2 (рис. 1), 
совершает гармонические колебания, которые возбуждаются вращением эксцентрика 3 с постоянной 
угловой скоростью. 

Предложенное устройство может быть эффективным при обработке фасонных деталей не самого сложного 
профиля, а в некоторых случаях и при обработке простых поверхностей, типа плоскостей или тел вращения. 

Следует прежде всего проверить возможность обеспечения широкого диапазона энергии воздействия 
инструментальной головки на поверхность обрабатываемой детали в сочетании с сравнительно низкими 
высотными характеристиками шероховатости поверхности. 

Ввиду меньшей жёсткости плоской пружины (в нашем случае двух пружин) в направлении Х, в сравнении с 
жёсткостью в направлении У, система описывает траекторию, близкую к эллипсу с большей полуосью в 
направлении оси Х. Для анализа закона движения системы разложим траектории по осям Х и У. 

Уравнение движения центра тяжести есть не что иное, как математическое выражение второго закона 
Ньютона. 


Рис. 1. Схема эксцентрикового упрочнителя: | — корпус; 2 — плоская пружина; 3 — эксцентриковая масса; 


4 — инструментальная головка; 5 — электродвигатель; 6 — гибкий вал; 7 — ограничитель (обрабатываемая деталь) 
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Рис. 2. Повороты корпуса упрочнителя в процессе колебаний 


Результаты исследования. Чтобы установить основные закономерности влияния технологических 
параметров процесса на взаимодействие осциллирующего индентора с обрабатываемой поверхностью 
необходимо учитывать кинетическую энергию индентора, количество инденторов, радиус индентора, 
коэффициент полезного действия устройства, физико-механические свойства обрабатываемой детали [1—6]. 

Учитывая все силы, действующие на движущийся индентор, можно записать уравнение движения 
инструментальной головки 

в направлении У: 


4?у ау 
= — о. 2 РВ 
В су-и рт +7, „Гб? С0$ Е — т, 8, (1) 
в направлении Х: 
2 Ах 
в т„иГ@ $1 Е, (2) 
т ГА 
4?у 4х 
где т, ть т, — проекции сил инерции системы на оси У и Х соответственно; ст, у, сх — проекции сил 
г” [2 
ах ау 
сопротивления пружины на оси У и Х; в —щ проекции сил сопротивления ро, РАДЫ на оси У их; 
ГА т 
ТтГ6? С08 СЕ, Т„тГо? т 6 — проекции возмущающейся силы на оси У и Х; т.е — сила тяжести (вес) 


вибрирующей системы; т. — масса вибрирующей системы; ст, с — жёсткость пружин в направлении У иХ; и — 
сопротивление среды; те» — масса эксцентрика; г— расстояние от оси вращения эксцентрика до его центра 
тяжести; & — угловая скорость; { — текущее значение времени; у, х — текущее значение координат. 

Ввиду значительно большей жёсткости пружин в направлении У по отношению с жёсткостью в направлении 
Х, амплитуда перемещения индентора в направлении У значительно меньше амплитуды в направлении Х. 
Поэтому примем допущение, что система совершает гармонические колебания только в направлении Х, т.е. 
будем рассматривать индентор только с одной степенью свободы. 

Если пренебречь затухающим колебанием, то уравнение движения будет иметь вид: 

х=ьЬзш( ат + В), (3) 

где Ь — амплитуда колебаний, В — разность фаз между возбуждающей силой и перемещениями центра 
тяжести индентора. 

Подставив это выражение в уравнение (2), находим: 


в т ига’ , (4) 
\/(<-— @?т,)* + аи’ 
@и 
пир (5) 
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Значение и определяется из выражений: 


@т,.б 
ый (6) 
2л 
где @, — частота собственных колебаний; 9— логарифмический декремент затухания. 


Продифференцировав уравнение (4) по времени и проведя исследование функции на экстремум, получим 
выражение для определения максимальной скорости индентора: 


з 
т.г 


= (т) 


У, : 
(с -?т, )? + 22 


Наибольшая кинетическая энергия индентора определяется из уравнения: 


— ту? р т.т? тг? 08 (8) 
2 2[(с — в?т, )? + 622] 


Анализ процесса взаимодействия сферических инденторов с деформируемым полупространством 
(поверхностным слоем обрабатываемой детали) описан в классических работах И.В. Кудрявцева [1, 2, 4, 5]. 
Диаметр пластического отпечатка можно определить по зависимости: 


(9) 


(10) 


где Т — кинетическая энергия инструментальной головки; Нр — динамическая твердость материала детали 
(отношение энергии удара сферического индентора к объему вытесненного материала при ударе); О; — диаметр 
индентора; и — коэффициент полезного действия устройства; М — число инденторов. 

При обработке эксцентриковым упрочнителем параметры шероховатости обрабатываемой поверхности 
могут получать постоянное (установившееся) значение, которое воспроизводится при дальнейшей обработке 
поверхности детали. Рельеф образующейся поверхности может быть как изотропным, так и анизотропным и 
формируется путем многократного наложения следов единичного взаимодействия. 

При взаимодействии осциллирующего индентора с исходными выступами микронеровностей происходит их 
уменьшение по высоте с одновременным уменьшением глубины впадин микронеровностей. С увеличением 
времени обработки исходный профиль шероховатости поверхности полностью передеформируется. В результате 
формируется новый микрорельеф, причем для каждого метода ППД он имеет специфический характер [7-19]. 

Окончательно сформированная шероховатость обрабатываемой поверхности получила название 
«установившейся». Как правило, еб высотные параметры не зависят от исходной. Они формируются в 
конкретных условиях каждого метода обработки и зависят от его технологических параметров. На основании 
методики работы [3, 4] получена зависимость для определения установившейся шероховатости поверхности при 
обработке эксцентриковым упрочнителем: 


Ка = 0,0075 (р 


Параметры упрочнения поверхностного слоя, к которым относятся глубина упрочненного слоя и степень 
деформации, оказывают важнейшее влияние на увеличение жизненного Цикла обработанных деталей. В 
результате теоретических исследований получены аналитические зависимости для их расчета при обработке 
эксцентриковым упрочнителем: 


(12) 


(13) 


Приведенные зависимости соответствуют физическому смыслу явлений, происходящих при обработке, и 
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прошли проверку в ходе проведения комплексных экспериментальных исследований. 
При проведении экспериментальных исследований использовались образцы из различных материалов, часто 
применяемые для изготовления деталей машин: качественные и легированные стали (сталь 45, ХВГ, сталь 30, 
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сталь ЗОХГСА и др.), алюминиевые сплавы (АЛ1, АВТ, Д16 и др.). Плоские образцы подвергались обработке 
эксцентриковым упрочнителем с различными режимами. Использовались шариковые и роликовые инденторы. 

По теоретическим зависимостям построены графики зависимостей шероховатости обработанной 
поверхности, глубины упрочненного слоя и степени деформации от режимов обработки, характеристик рабочих 
сред и обрабатываемых материалов. 

На графиках (рис. 3—6) сплошной линией показаны кривые, построенные по теоретическим формулам, а 
точками — результаты экспериментальных исследований. Выполнено построение доверительных интервалов с 
доверительной вероятностью 95 %. 

Наблюдается высокая сходимость результатов, что свидетельствует о том, что теоретическая зависимость 
правильно отражает явления, происходящие при обработке ПИД осциллирующим инструментом — 
эксцентриковым упрочнителем. 


Ва, пис 


0,48 2 1 
0,46 
0,44 
0,42 
0,40 
0,38 


| | |— 
2,0 ов 3.0 3,5 4.0 О, шт 


Рис. 3. Зависимость шероховатости поверхности от радиуса индентора: 
1 — материал образца сталь 45, 2 — материал образца ХВГ 
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Рис. 4. Зависимость шероховатости поверхности от твердости детали по Бринеллю 
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Ь, шт А 5 № А 
0,25 т 9 
0,20 + 81 
0,15 + 71 
0,10 - 6 + 
0,05 | | ре то 5 
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а) 6) 
Рис. 5. Зависимость параметров упрочнения от твердости по Бринеллю для различных материалов: 
а — от глубины упрочненного слоя; б — от степени деформации 


В, шт А 
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| 
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а) 6) 
Рис. 6. Зависимость параметров упрочнения от диаметра индентора: а — глубина упрочненного слоя; 
б — степень деформации; 1 — материал образца — сталь 45; 2 — материал образца — сталь ХВГ 


Обсуждение и заключение. По результатам проведенных исследований можно сделать следующие выводы: 

1. Получена теоретическая зависимость, позволяющая определить кинетическую энергию индентора при 
обработке осциллирующим инструментом — эксцентриковым упрочнителем. 

2. Получены зависимости для определения диаметра и глубины пластического отпечатка, а также 
шероховатости поверхности по параметру Ва, глубины упрочненного слоя и степени деформации, позволяющие 
прогнозировать качество обработанной поверхности. 

3. Выполнено сравнение результатов теоретических и экспериментальных исследований процесса обработки 
эксцентриковым упрочнителем. Расхождение результатов не превышает 15 %. 

4. Исследуемый процесс подлежит дальнейшему изучению с целью определения других параметров качества 
обработанной поверхности, например, величины остаточных напряжений в поверхностном слое, а также 
расширения диапазона режимов и конструктивных параметров обработки с целью определения их оптимальной 
области [5, 6, 20]. 

5. Полученные зависимости для определения основных параметров качества поверхностного слоя дают 
возможность прогнозирования результатов обработки и могут быть использованы для проектирования 
технологических процессов эксцентриковым упрочнителем. 
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Аспекты теплозащиты машиностроительного и энергетического 
оборудования: применение стойких к оксидированию 
комбинированных покрытий на основе никеля 


В.Н. Варавка© ‚ О.В. Кудряков, В.И. Г рищенко@ 
Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 
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Аннотация 

Введение. В тех областях энергетического машиностроения, где используется тепловая энергия перегретого 
пара, важным аспектом обеспечения надежности и безопасности оборудования является теплостойкость 
используемых материалов. При изготовлении индукционных пароперегревателей оптимальным материалом для 
паропровода (змеевика) является медь. Однако её предельная стойкость к оксидированию не превышает 400 °С, 
что существенно ограничивает эффективность работы парогенераторов. Поэтому целью работы было 
исследование кинетики окисления комбинированного гальванического покрытия системы Мо-№-Сг, 
нанесенного на медные трубчатые образцы и предназначенного для теплозащиты змеевиков парогенераторов. 
Материалы и методы. На опытных медных трубчатых образцах было сформировано комбинированное 
гальваническое покрытие системы Мо-М№1-Сг с общей толщиной 12-35 мкм. Подслой Мо толщиной около 
1,5 мкм на поверхности медной трубки был сформирован для предотвращения диффузии Си в №-покрытие. 
Слой хрома толщиной 1,5 мкм на поверхности покрытия выполнял роль индикатора процесса окисления. 
Сравнительный анализ процессов окисления поверхности меди и комбинированного покрытия системы Мо-М№- 
Сг на медной подложке выполнен с использованием методик оптической и электронной микроскопии, 
энергодисперсионного анализа, а также прецизионного определения параметров роста оксидных пленок. 
Результаты исследования. Экспериментально определены интервалы термической устойчивости медной 
подложки и никелевого покрытия. Полученные экспериментальные зависимости характеризуют 
параболический закон окисления меди с образованием однофазной диффузионной зоны СчО при температурах 
выше 350 °С и никеля при температурах выше 750 °С, когда начинается переход монооксида №0 и в оксид 
№03. Рост оксидных пленок по квадратичным законам приводит к быстрому увеличению толщины пленок, 
накоплению в них напряжений, растрескиванию и скалыванию. 

Обсуждение и заключение. Показано, что гальваническое покрытие Мо-М№1-Сг устойчиво к нагреву при 
длительной эксплуатации вплоть до температур 750-800 °С. Описаны функциональные роли Мо и Сг в 
архитектуре покрытия. Работа акцентирована на прикладном аспекте использования исследуемого покрытия 
для повышения термической устойчивости  змеевика-паропровода промышленных индукционных 
пароперегревателей малой и средней мощности. 


Ключевые слова: парогенераторы, теплостойкость, окислительный процесс, гальванические покрытия, 
микроструктура, электронная микроскопия, гравиметрический анализ 
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АБЯгас 

Гитодисйоп. ш Фе агеаз оЁ ро\ег епзшееппе \Веге фе Шегта| епегоу оЁ зарефезе4 з{еат 15 изе4, ап ипрогапе 
азресЕ оЁ ргоу1Атх ше тепа Шу ап4 заеу оЁ едшртепе 1$ Ше Беа{ гез1$апсе оЁ Фе тшаена!$ етроуе4. ш Ше 
тапиасваге оЁ шдисНоп зиареезег$, ше орйта! таепа| Рог Ше $еат р1ре (соП) 1$ соррег. Но\меуег, из иШитае 
тез15ёапсе 1ю ох14аНоп 40е$ по{ ехсее4 400 °С, ушсЬ э1етиЙсапИу Ши фе еРйс1епсу оЁ ${еат сепегаогз. ТБегеоте, 
Фе оБ]фесйуе оЁ Ше \огК \аз 1№ю зшАу Фе Кшейс$ оЁ ох1Чавоп оЁ Ше сотЫтей са]уапс соаНпе оЁ Фе Мо-М№-Сг 
зузет аррЦе4 то соррег аб \Щаг затр[ез ап4 ищеп4е4 юг Фегта| ргоесНоп оЁ %еат зепегаюг сой$. 

Маептаб ап Мефо45. А сотЫпе4 еескор!айпе оЁ фе Мо-М1-Сг зубет \/ИБ а ю{а1 ЧусКпез$ оЁ 12—35 т \/а$ 
Юге оп Фе ехрегитетйа! соррег таб \аг затр[ез. А Мо зи ауег \Ий а @1сКпез$ оРаБоци 1.5 ит оп Фе загсе оЁ 
Фе соррег бе \уаз Югте4 {0 ргеуеп ве А а51оп оЁ Си шю Ше № соайпе. А 1.5 им ЧК сБготпит 1ауег оп Фе 
соаНпс заГасе асе4 аз ап саг оЁ Ше ох14аНоп ргосезз. А сотрагануе апа[у51$ оЁ Ше ох!ЧаНоп ргосеззез оЁ ше 
соррег зи!{Расе ап Фе сот ше соайпе оЁ Фе Мо-М№-Сг зумет оп а соррег забзгае \а$ сагмей оиё изше Ше 
тефо4$ оЁ орИса| ап4 @есйоп писгозсору, епегсу Чзрегяуе апа[уз1$, ап ргес151оп 4еегимтаноп оЁ фе это\ив 
рагатеег$ оЁ ох1ае Н|пз. 

Кези5. Тре пцегуа|з оЁ Фегта! за Шу оЁ Фе соррег забзгае ап4 п1сКе! соанпе уеге ехрегитещаПу деегилтеч. 
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Введение. Электрохимическое осаждение металлов широко распространено в промышленности, являясь 
основой гальванотехники. К одной из особенностей развития этой отрасли науки принято относить то, что её 
развитие проходило «почти исключительно эмпирически» [1], отталкиваясь в первую очередь от потребностей 
различных отраслей индустрии. Несмотря на то, что в настоящее время теоретические основы электрохимии 
проработаны достаточно глубоко [2—6], прикладной аспект здесь по-прежнему имеет приоритетное значение и 
определяет большинство решаемых научных задач, связанных, в основном, с особыми условиями 
использования гальванических покрытий. В таких особых условиях эксплуатируются, например, 
пароперегреватели, работа которых связана со значительным изменением состава и температуры пара по длине 
змеевика-паропровода [7]. В зависимости от мощности парогенератора температура змеевика по его длине 
может изменяться от 150 до 650°С, а для суперсверхкритических параметров пара в особо мощных 
современных паровых турбинах и выше [8-10]. Настоящая работа посвящена исследованию возможности 
применения гальванических покрытий для защиты паропровода индукционного пароперегревателя от 
оксидирования при высоких температурах. Интерес к теплозащитным покрытиям вызван тем, что в отношении 
бытовых и малой мощности промышленных парогенераторов по совокупности физико-технологических 
свойств (электропроводность, теплопроводность, способность к пластической деформации, обрабатываемость 
резанием и т.д.) медь является на данный момент незаменимым конструкционным материалом для 
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изготовления змеевиков. Однако окислительная способность меди тоже высока и её сопротивление 
оксидированию не превышает 400 °С. Исходя из условий работы рассматриваемых пароперегревателей, это 
обстоятельство ставит задачу применения покрытий, у которых уровень стойкости к окислению находится выше 
теплового барьера 600 °С. 

Материалы и методы. Учитывая сложную конфигурацию змеевика, наличие большой протяженности 
криволинейных поверхностей и его немалые габаритные размеры, в качестве наиболее технологичного способа 
нанесения теплозащитного покрытия было выбрано электрохимическое осаждение. 

Для выбора состава такого покрытия необходимы точные данные о режимах его эксплуатации. С этой 
целью был проведен тепловизионный анализ термических условий работы экспериментального индукционного 
трехкатушечного 6-ти виткового парогенератора мощностью 10 кВт (рис. 1) [11]. Паропровод был выполнен из 
профилированной медной трубки ©25х1,5 мм из технической меди марки М2 по ГОСТ 617-2006. 


Рис. 1. Экспериментальный индукционный пароперегреватель с медным паропроводом: 


а— общий вид; б — наиболее сильно окисленные участки паропровода (змеевика) на выходе пара 


Количественный термический анализ условий эксплуатации паропровода, включая форсированные режимы 
работы, был выполнен с помощью бесконтактного тепловизора модели ЕшКе Т1401 РКО (производство НаКе 
Сотр., ЗА) [12] с основными техническими характеристиками: 

— инфракрасный спектральный диапазон: от 7,5 до 14 мкм (длинноволновый); 

— тепловая чувствительность: <0,075 °С при температуре объекта 30 °С (75 мК); 


— погрешность: =2 °С (при низких температурах) или 2 %; 

— степень защиты: по ГОСТ 14254-96 (МЭК 60529): 1Р54. 

Теплоёмкость водяного пара (с, = 33,6 Дж/моль-К) при нормальных условиях) примерно вдвое ниже 
теплоёмкости воды (с, = 75,3 Дж/мольК), что существенно меняет условия теплоотвода в змеевике и 
способствует интенсификации окисления поверхности паропровода, поэтому по его длине наблюдается 
значительная неоднородность степени окисленности (рис. 1 6). Результаты термического анализа (рис. 2) 
показали, что максимальные температуры нагрева фиксируются на выходном полувитке змеевика (рис. 16). 
Интервал их значений составил 530-540 °С с абсолютным максимумом при 541,38 °С. Средние значения 
температуры большинства перегретых (окисленных) полувитков находились на уровне 420-460 °С. 

Исходя из полученных результатов, в дальнейших исследованиях по получению теплостойких 
покрытий для защиты от оксидирования поверхности змеевика необходимо ориентироваться на 
предельную температурную нагрузку 600 °С. В связи с этим целесообразно использование покрытий на 
основе никеля. № составляет основу большинства современных жаропрочных суперсплавов, используемых 
в теплоэнергетике [13—15], а технология нанесения гальванического № достаточно хорошо отработана [15—18]. 

При нанесении экспериментальных гальванических покрытий на образцы медных трубок технической 
меди М2 были использованы стандартные режимы осаждения и составы электролитов, содержащих Ми 
Ст, рекомендуемые ГОСТ 9.305 и 9.306. При эксплуатации парогенератора покрытие практически не 
подвергается механическому воздействию, поэтому при работе оно не должно демонстрировать выдающиеся 
механические свойства. В то же время, при нанесении на криволинейную выпуклую поверхность в покрытии 
формируются внутренние растягивающие напряжения. По мере увеличения количества рабочих теплосмен их 
уровень возрастает. В связи с этим толщина исследуемого покрытия на паропроводе не должна быть слишком 
большой. Она принималась средней из рекомендуемых в литературе интервалов значений для никелевых 
гальванических покрытий, выполняющих защитно-декоративные функции, и составляла ориентировочный 
уровень 20 мкм. 
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Рис. 2. ИК-изображение общего вида тепловых полей работающего пароперегревателя (окно программы «Е Ке Соппесф», 
расположение змеевика аналогично изображенному на рис. 1 а (все значения — в градусах Цельсия) 


Для изучения микроструктуры покрытий и топографии их поверхности использовался двулучевой 
(электронный/ионный) сканирующий электронный микроскоп /Е]$$ СгоззВеат 340 (ЗЕМ), позволяющий с 
помощью ионного пучка (ЕВ) проводить травление и выполнять кросс-секции (сечения заданной 
конфигурации) образцов непосредственно в вакуумной камере микроскопа с высокой точностью 
позиционирования. Элементный состав исследуемых поверхностей контролировался с помощью 
энергодисперсионного рентгеновского детектора (ЕРАХ) модели Х-Мах 50М (ОхЕЮга шзилитеп $), встроенного 
в электронный микроскоп. 

Исследование кинетики оксидирования покрытий было выполнено путем измерения массового показателя 
высокотемпературной газовой коррозии (прирост массы образца в результате нагрева и окисления). Для 
определения степени окисления медных образцов как с покрытием, так и без покрытий, все образцы 
взвешивались до и после опытов на разных этапах взаимодействия. Гравиметрические исследования 
проводились на аналитических весах марки «ВЛР-20» с точностью взвешивания 10- г. 

Результаты исследования 

1. Качественный анализ кинетики окисления. Особенности химического взаимодействия никеля с 
кислородом проявляются в том, что № образует две модификации монооксида: о-№Ю с гексагональной 
решеткой (ниже 252 °С) и В-\№О с гранецентрированной кубической решеткой. Переход происходит при 
непрерывном нагреве в интервале 250-300 °С. Экспериментально было установлено, что при нагреве до 630 °С 
протекает диффузионный процесс через тонкую пленку монооксида №0, выше 640°С устанавливается 
химический процесс образования №0, который при нагреве выше температур 800 °С может приводить к 
формированию оксида №20: [19]. 

Прикладной аспект использования теплостойкого никелевого покрытия на медной подложке осложняется двумя 
обстоятельствами: неограниченной растворимостью компонентов системы Си-М (рис. 3) и возможностью эффекта 
взаимодиффузии Киркендалла [21, 22] на границе «№-покрытие — Си-подложка». Эти явления, ввиду растворения 
меди в покрытии, работают в направлении снижении стойкости покрытия к оксидированию. 
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Рис. 3. Диаграмма фазовых состояний системы «медь-никель» [20] 
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С целью исключения негативного влияния указанных обстоятельств на устойчивость покрытия к 
оксидированию на опытных медных трубчатых образцах было сформировано комбинированное гальваническое 
покрытие системы Мо-М№-Сг (рис. 4). Подслой Мо толщиной около 1,5 мкм на поверхности медной 
трубки (рис. 46) был сформирован для предотвращения диффузии Си в №-покрытие при длительной 
эксплуатации паропровода ввиду практически нерастворимой системы компонентов Си-Мо и ограниченной 
растворимости №-Мо [20]. Слой хрома толщиной 1,5 мкм на поверхности покрытия (рис. 4 6) выполнял роль 
индикатора процесса окисления (подробнее об этом см. далее). Общая толщина покрытия на 
экспериментальных образцах с покрытиями составляла 12-35 мкм. Элементное распределение в поперечном 
сечении исходного покрытия Мо-К№-Сг (до начала эксперимента с нагревом образцов) представлено на рис. 5. 


а) 6) 


Рис. 4. Исходное покрытие Мо-№-С в поперечном сечении, ЗЕМ: а — покрытие с маркерами толщины; 


6— однородная микроструктура покрытия и водородная пористость на границе с подложкой 
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Рис. 5. Цветовые карты распределения химических элементов по глубине покрытия Мо-М1-Сг, ЕРАХ: 
а — общий вид покрытия в поперечном сечении (ЗЕМ); 
6-0 — распределение элементов на снимке общего вида Сг (6), Си (в), Мо (г), № (д) 


Для экспериментального изучения кинетики окисления в условиях, максимально приближенных к условиям 
эксплуатации паропровода, образцы для исследования были изготовлены из медных трубок с 
соответствующими толщиной стенки и диаметром (рис. 6). 

После нанесения покрытия по технологии, включающей элементы термической обработки [11], покрытие 
приобретает зеленоватый оттенок, характерный для монооксида никеля №0. При эксплуатации паропровода 
покрытие наносится только на наружную (выпуклую) поверхность образцов. Однако для корректного 
определения прироста массы покрытия на экспериментальные образцы покрытие наносилось с обеих сторон. 

Одновременный нагрев образцов из чистой меди и образцов с покрытием проводился при фиксированной 
температуре в интервале 350-1000 °С в печи $МОГ, 6,7/1300 (2,4 КУ!) в воздушной атмосфере. Выдержка при 
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заданной температуре составляла 30 мин. Для статистической картины эксперимента нагрев при каждой 
заданной температуре осуществлялся для 5—7 образцов при их раздельной загрузке в печь. Избирательные 
результаты эксперимента визуализированы на рис. 7. 

Из полученных экспериментальных данных следует, что в условиях проведенного нагрева медь 
относительно теплостойка до температуры 300-350 °С. При этих температурах на медной поверхности, 
независимо от еб кривизны, формируется плотная тонкая оксидная пленка коричневого цвета (рис. 7). 


Рис. 6. Образцы для изучения кинетики окисления; показаны внешние поверхности эталонных образцов чистой меди 
(справа) и меди с покрытием Мо-М№-Сг (слева), подготовленные для проведения экспериментов 


их 2% =- 
ШИ ... 
а #/ ° 


Рис. 7. Сравнение наружной поверхности образцов после нагрева до указанных температур: 
слева — образцы с покрытием, справа — образцы из чистой меди 
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Рис. 8. Медная поверхность образца, оксидированного при температуре 650 °С, $ЕМ: 


2нт 
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а — расположение кристаллов СчО на поверхности меди; б — морфология кристаллитов СиО 


В соответствии с литературными данными [23, 24], она, по-видимому, является двухслойной: на поверхности 
образца расположен тонкий подслой Си20, а снаружи находится слой СиО. В силу малой толщины последнего, 
внутренние напряжения в пленке невелики. Она имеет хорошую адгезию с подложкой, низкую шероховатость, 
не разрыхляется и не скалывается с поверхности. При повышении температуры нагрева ускоряется рост 
внешнего оксидного слоя СчО. Уже при температуре 450 °С он становится весьма рыхлым и осыпается с 
поверхности (рис. 1 6). При этом под ним обнаруживается подслой Си›О характерного красноватого 
оттенка (рис. 7). В процессе дальнейшего нагрева шелушение и осыпание слоя оксида СиО прогрессирует — он 
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практически не держится на поверхности вплоть до температур 650—700 °С. Это, по-видимому, обусловлено 
характером кристаллизации монооксида меди: его кристаллиты имеют строгую огранку, близкую к 
кубической (рис. 86), слабое сопряжение друг с другом и, главное, высокую неоднородность мест 
зарождения (рис. 8 а). Начиная от температуры нагрева =750 °С, оксидная пленка меди уплотняется, 
прочность её сцепления с поверхностью повышается. На вогнутой поверхности образцов за счет сжимающих 
конфигурационных напряжений пленка СчО держится достаточно прочно или может отслаиваться целиком со 
всей поверхности образца не рассыпаясь. При температурах 800-900 °С оксидная пленка СиО ведет себя 
аналогично и на внешней (выпуклой) поверхности образцов. 

Поверхность образцов с покрытием Мо-М№-Сг практически не меняется вплоть до температуры нагрева 
750 °С. При дальнейшем нагреве выше 800 °С покрытие сначала покрывается пленкой оксида Сг›Оз, имеющим 
характерный ярко-зеленый цвет (рис. 7). Затем, при температурах нагрева выше 900 °С, окисляются и более 
глубокие слои покрытия. Амфотерный оксид хрома Сг›Оз (рис. 96) по своей морфологии и характеру 
кристаллизации принципиально отличается от монооксида меди СчО (рис. 8). Кристаллы оксида Сг›Оз имеют 
характерную форму многогранника с преобладанием кристаллитов призматической формы. Благодаря 
большому разбросу кристаллитов по размерам, они, в отличие от монооксида меди, формируют на 
поверхности слой с высокой плотностью укладки кристаллитов. 
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Рис. 9. Поверхность образцов с покрытием Мо-М1-Сг, ЗЕМ: 
а — после нагрева до температуры 650 °С; б — после нагрева до температуры 850 °С 
Если до начала окисления слоя хрома (ниже 800 °С) поверхность покрытия имеет весьма слабовыраженный 
характер кристалличности (рис. 9 а), то появление оксида Сг›Оз придает поверхности хрестоматийный 
поликристаллический вид (рис. 9 6). В целом по сечению никелевое покрытие свой состав и строение не меняет 
вплоть до температуры 850 °С (рис. 10). Химический состав поверхности покрытия при этой температуре 
свидетельствует о начальной степени окисления слоя хрома. Наличие тонкого слоя окиси хрома на поверхности 


подтверждается данными энергодисперсионного анализа (ЕРАХ) как по глубине покрытия, так и по 
поверхности. Из рис.11 видно, что кислород сосредоточен в значительно более узком поверхностном 
слое (-1 мкм), чем слой хрома (-3 мкм), что характеризует начальную стадию окисления. 
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Рис. 10. Толщина и структура покрытия Мо-М№-Сг в поперечном сечении после нагрева до 850 °С, ЗЕМ 
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В целях сокращения объема статьи данные ЕБАХ приводятся не полностью. Его результаты 
свидетельствуют о том, что количество кислорода на поверхности возрастает с 30 до 50 ат. % за счет 
аналогичного снижения концентрации хрома, что указывает на окисление именно хрома (т.к. 
концентрация № при нагреве от 650 °С до 850 °С не меняется, а присутствие № в результатах 
детектирования обусловлено проникновением рентгеновского излучения при ЕРАХ-анализе сквозь тонкий 
слой хрома в никелевую основу покрытия). Состав оксидов соответствует соединению Сг2Оз. 


8 
ИЕ 


Сг, вес. % 


30 


20 


10 


> 
— 
>) 
>) 
> 
95 
(<>) 
> 
> 
бл 
> 


Ь, мкм 


№, вес. % 
100 


ма эх 
> <> <> <) 


> 


0 10 20 30 40 50 Ь, мкм 


г) 
Рис. 11. Распределение основных химических элементов в покрытии Мо-М1-Сг по глубине В после нагрева до 850 °С, 
ЕЛРАХ: а — покрытие в сечении, показано направление сканирования (ЗЕМ); б-г — содержание элементов по 
направлению сканирования: б — кислорода; в — хрома; г — никеля 


Машиностроение и машиноведение 


Все стадии процесса окисления покрытия представлены на рис. 12. При температуре выше 800 °С 
поверхность покрытия начинает окисляться, на что указывает изменение его цвета — покрытие приобретает 
цвет «кенсингтонской зелени». 
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а) 
Рис. 12. Последовательные стадии окисления никелевого покрытия, оптическая микроскопия, х100: 
а — начальная стадия: прорастание оксида никеля (темно-серые выделения) на поверхность оксида хрома (зеленое поле) 


при 850 °С; 6 — заключительная стадия: собирательная фазовая картина в области скола покрытия при 1000 °С, где: 
1 — оксид хрома Сг2Оз; 2 — оксид никеля №2Оз; 3 — поверхность меди (скол покрытия); 4 — край образца 


Этот цвет соответствует оксиду хрома Сг>Оз. Уже при температуре 850 °С на поверхности покрытия можно 
встретить редкие единичные образования оксида никеля М№2Оз (на рис. 12 а серого цвета), которые при 
дальнейшем нагреве постепенно увеличивают свою площадь, занимаемую на поверхности. По мере их роста, 
который преимущественно распространяется тангенциально поверхности, под действием внутренних 
напряжений оксидная пленка №2Оз трескается и затем скалывается, обнажая поверхность подложки — чистой 
меди красноватого оттенка (рис. 12 6). 

Таким образом, основным результатом этой части исследований следует считать экспериментально 
установленные температурные интервалы допустимого применения материалов для изготовления паропровода 
парогенератора. Так, медный паропровод без покрытий работоспособен до температуры 300 °С и может быть 
использован только для генерации влажного (не перегретого) пара. Начиная с температуры -400 °С и до 
—7Т00 °С формирующаяся оксидная пленка СиО очень рыхлая и легко осыпается с поверхности меди. При более 
высоких температурах пленка уплотняется, становится толще и адгезия её к медной подложке возрастает. 
Однако она по-прежнему склонна к скалыванию в процессе теплосмен. Её наличие на поверхности паропровода 
существенно замедляет теплоотвод, а продолжающиеся при нагреве химические реакции высокотемпературной 
газовой коррозии работают в направлении уменьшения толщины трубопровода. Ввиду неоднородности 
протекающих процессов эксплуатация паропровода в таких условиях становится непредсказуемой с точки 
зрения аварийных ситуаций. Использование комбинированного гальванического покрытия системы Мо-М№1-Сг 
повышает работоспособность паропровода до температуры 750-800 °С. При достижении температуры -850 °С 
покрытие начинает окисляться наравне с медью. При 950 °С и выше окисленное покрытие проявляет 
склонность к скалыванию и его эксплуатация подвержена тем же рискам, что и медный трубопровод. 
Отличительной особенностью исследованного покрытия является самотестирование: если в процессе 
эксплуатации температура нагрева превышает 800 °С, то поверхностный слой хрома окрашивает покрытие в 
ярко-зеленый цвет и сигнализирует об опасности перегрева. Индикаторный слой гальванического хрома после 
окисления может быть легко восстановлен и эксплуатация паропровода после этого продолжится. 

Количественный анализ кинетики окисления 

В рамках выполненных исследований был проведен также количественный анализ кинетики оксидирования 
образцов из чистой меди и образцов с покрытием Мо-М№-Сг. В качестве измеряемого параметра использован 
наблюдавшийся в процессе нагрева удельный прирост массы М=Аии$, где Ат — увеличение массы образца, г; 
$ — площадь окисленной поверхности образца, см?. В ходе нагрева трубчатых образцов из чистой меди и 
образцов с покрытием, проведенного по описанной выше методике качественного анализа, были получены 
экспериментальные данные значений М, представленные в таблицах | и 2. 

В таблицах приведены интервалы разброса полученных данных для фиксированных значений 
температур нагрева (таблица 1) или времени выдержки в печи (таблица 2), а также среднее значение М из 
каждого интервала. 

Статистическая обработка данных, приведенных в таблицах 1 и 2, выполненная с помощью пакета 
прикладных программ МаСАР, включавших процедуры интерполяции, позволила получить кинетические 
зависимости, показанные на рис.13 и 14. Прямолинейные графики полученных зависимостей, 
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представленные на рис. 13 6 в «аррениусовых координатах» (-м М — 1000/Т), характеризуют параболический 


закон окисления меди при температурах выште 350 °С и никеля при температурах выше 750 °С [19, 25—27]. 


Удельный прирост массы М, г/см? 


Экспериментальные данные удельного прироста массы образцов 
при печном нагреве в воздушной атмосфере в течение 30 мин 


Таблица 1 


№ 


Удельный прирост массы М, 10-5 г/см? 


Температура 
эксперимента р Медная трубка : 
в печи, °С Медная трубка с №-покрытием* 
п/п без покрытия 
1 350 1,35 + 0,31 — 
2 450 6,23 + 1,37 _ 
3 550 24,73 + 3,60 в 
4 650 58,51 + 9,02 1,52 + 0,24 
5 750 138,88 + 17,85 4,25 + 0,58 
6 850 241,03 + 25,25 13,12 + 1,08 
7 1000 564,70 + 49,76 35,71 2,78 


*за критерий отсутствия окисления (прочерк в таблице) принимался минимальный, измеряемый на 


используемых аналитических весах, прирост массы, равный 105 Г 


Экспериментальные данные удельного прироста массы образцов 
при различном времени выдержки в воздушной атмосфере печи 


Таблица 2 


Удельный прирост массы М, 10-$ г/см? 
№ экспери- Время выдержки м тит м кит 
ыы и едная трубка без едная трубка с №1-покрытием при 
покрытия при 600 °С 800 °С 
1 5 13,68 = 2,02 3,03 + 0,41 
2 15 29,43 + 4,15 5,93 = 0,54 
3 30 35,76 = 4,84 7,11 0,67 
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Рис. 13. Температурные зависимости прироста массы М образцов чистой меди (1) и образцов 
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с покрытием Мо-М1-Сг (2): а — в абсолютных единицах; б — в относительной координатной системе 
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Рис. 14. Кинетика изменения во времени прироста массы М образцов чистой меди: 
1 — при 600 °С; 2 — образцов с покрытием Мо-М№-Сг при 800 °С 


Рост оксидных пленок по квадратичным законам происходит с образованием однофазных диффузионных 
зон, в рассматриваемом случае состоящих из оксидов СО и №0 соответственно. Он приводит к быстрому 
увеличению толщины пленок, накоплению в них напряжений, растрескиванию и скалыванию. Дополнительный 
вклад в ускорение этого процесса вносит криволинейная наружная поверхность медной трубки [28, 29]. 

Обсуждение и заключение 

1. Выполненный комплекс исследований показал, что гальваническое комбинированное покрытие системы 
Мо-М№1-Сг является достаточно эффективной защитой медного паропровода от окисления. Покрытие способно 
обеспечить длительную работоспособность парогенератора вплоть до температур нагрева 750-800 °С. 

2. Длительная теплостойкость покрытия обеспечивается внутренним слоем № с рекомендуемой толщиной 
20-30 мкм. Исследование кинетики оксидирования покрытия, выполненное методами оптической и 
электронной микроскопии, энергодисперсионного анализа, а также с использованием методик прецизионного 
определения параметров роста оксидных пленок, показало, что никелевое покрытие индифферентно к нагреву 
вплоть до температур 600-650 °С. В интервале температур 700-900 °С окисление покрытия происходит с 
образованием монооксида №0 по параболическому закону. При более высоких температурах окисление 
прогрессирует за счет формирования пленки оксида №2Оз, что достаточно быстро приводит к её росту, 
растрескиванию и скалыванию. 

3. Комбинированная архитектура исследованного никелевого покрытия включает два тонких слоя Мо и Сг. 
Подслой Мо толщиной около 1,5 мкм расположен на поверхности медной трубки (подложки). Его функция — в 
предотвращении взаимной диффузии № и Си на границе «покрытие-подложка» при длительной эксплуатации 
парогенератора, поскольку растворение меди снижает термостойкость никеля и нарушает работоспособность 
покрытия. Наружный слой хрома толщиной 2-3 мкм выполняет роль индикатора степени окисления покрытия. 
Первым признаком избыточной окисленности покрытия является появление ярко зеленого оттенка на 
поверхности покрытия, что связано с образованием оксида хрома Сг>Оз при температурах >800 °С. Индикатор 
перегрева — слой хрома — легко обновляем и способствует лонгированию жизненного цикла парогенератора. 
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Аннотация 

Введение. Оценка производственной технологичности изготавливаемых изделий — этап технологической 
подготовки и ключевой аспект развития современных производственных механообрабатывающих систем — 
является актуальной задачей современного машиностроения. В этой связи теоретические и практические 
исследования по разработке методических подходов к определению весовой значимости количественных 
показателей при оценке производственной технологичности деталей являются весьма актуальными. Целью 
представленной работы явилась разработка метода оценки, нацеленного на повышение качества обработки 
деталей и эффективности функционирования многономенклатурных производственных систем на основе 
разработки дополнительных количественных показателей оценки производственной технологичности. 
Материалы и методы. Для оценки влияния количественных производственных показателей, связанных с 
затратами времени при простое оборудования, создана модель прогнозирования потоков событий поступления 
партий деталей на изготовление на определенную операцию и потоков обработанных деталей с использованием 
аппарата теории массового обслуживания. Такой подход позволяет учесть, как конструкторско-технологические 
характеристики деталей, особенности конкретной производственной системы, так и складывающуюся 
производственную ситуацию. 

Результаты исследования. Посредством оценки возможного влияния на составляющие при расчете штучно- 
калькуляционного времени (Тик) на уровне технологической операции была определена степень влияния 
показателей технологичности. Установлены взаимосвязи между показателями технологичности и затратами по 
всем статьям технологической себестоимости обработки заготовки (Соп), а также затратами, связанными с 
организационными простоями оборудования (Сир.о.). С помощью применения аппарата парных сравнений при 
принятии решений применительно ко всем структурным элементам производственных затрат определена 
степень влияния показателей технологичности относительно других показателей. 

Обсуждение и заключение. Описан подход к выполнению данной проектной процедуры, позволяющий 
учитывать состав и возможности технологического оборудования конкретного производства и реально 
складывающуюся производственную ситуацию. Разработанные формализованные модели позволяют 
комплексно спрогнозировать влияние показателей технологичности деталей на эффективность 
функционирования механообрабатывающих систем при их изготовлении. 


Ключевые слова: технологическая подготовка производства, оценка технологичности изделий, количественные 
показатели производственной технологичности, механообрабатывающие производственные системы, 
эффективность функционирования производства 
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Введение. Развитие машиностроительного производства в современных условиях невозможно без серьезного 
наращивания научных исследований, связанных с разработкой теории и методических принципов формализации 
всех этапов производства изделий, являющихся основой будущих интеллектуальных систем поддержки создания 
и изготовления технических объектов. Наиболее сложными в этом аспекте выступают задачи конструкторской и 
технологической подготовки производства [1, 2]. Несмотря на многочисленные работы в данном направлении, 
до настоящего времени не созданы автоматизированные системы, позволяющие взять на себя даже 
незначительные функциональные действия конструкторов и технологов, связанные с выполнением творческих 
проектных решений. 

Обязательным условием технологической подготовки эффективного функционирования 
механообрабатывающих систем является мониторинг и анализ складывающейся производственной ситуации, а 
также информации о состоянии оборудования и средствах технологического обеспечения. Только обладая в 
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полной мере этими сведениями можно принимать рациональные производственные решения. Даже опытный 
технолог не в состоянии собрать и проанализировать такой большой объем информации. Поэтому решения 
принимаются зачастую субъективно и необоснованно, проектирование технологических процессов и их 
реализация разнесены во времени, а использование вычислительных систем тормозится отсутствием моделей, 
описывающих процесс технологической подготовки производства. 

Решению поставленных задач посвящены научно-исследовательские работы по созданию системы 
планирования многономенклатурных технологических процессов [3], в основе которых лежит концептуальный 
подход к формализации всех проектных процедур обеспечения технологической подготовки 
механообрабатывающих производств, позволяющих учесть конкретные особенности, технологические 
возможности оборудования и оснастки. Одной из таких проектных процедур является оценка обеспечения 
технологичности изделий, чему традиционно уделяется недостаточно внимания. Роль данного этапа существенно 
недооценена. 

Все технико-экономические показатели функционирования производственной системы определяются 
сложностью выпускаемой продукции и степенью оснащенностью производственных мощностей. Часто имеет 
место несогласованность этих двух показателей, что приводит к невозможности обеспечения требований к 
качеству выпускаемой продукции, простоям и нерациональному использованию оборудования. Именно при 
оценке технологичности изделий, наряду с известными задачами и методами их решения, должны быть получены 
объективные данные о целесообразности изготовления изделий в условиях конкретной производственной 
системы. 

Материалы и методы. Большое значение имеют научные исследования в области создания 
формализованных моделей установления связей между технологическими и конструкторскими задачами 
подготовки производственных машиностроительных систем. В связи с увеличивающимся мировым 
соперничеством в производственной сфере, первоочередной задачей является повышение эффективной работы 
оборудования при реализации технологических процессов с учетом соответствия заданных требований к 
качеству деталей, которые, в свою очередь, устанавливаются в процессе конструирования. 

Вопросы, поставленные современными условиями функционирования производственных комплексов, 
связаны с обеспечением производственной технологичности изготавливаемых изделий. В настоящее время 
методы оценки технологичности изготавливаемых изделий, с учетом необходимости соответствия требованиям 
стандартов, напрямую зависят от квалификации технолога (конструктора) и их знаний. Такой подход не 
гарантирует рационального принятия технических решений и сдерживает автоматизацию проектных процедур. 

Оценка технологичности, как этап подготовки производства, выполняется для установления взаимосвязи 
между затратами на изготовление изделия и его конструктивными особенностями. Результаты такой оценки 
зачастую имеют противоречивый характер, отсутствует полноценное математическое описание процедуры 
ее проведения. 

Для разрешения сложившейся ситуации целесообразным представляется реализация на практике следующих 
шагов [4—6]: 

— установление относительных весовых характеристик показателей технологичности на основе параметров 
изготовляемых изделий. Труднореализуемым, но первостепенным является решение этой задачи на стадиях 
разработки рабочей конструкторской документации, когда отсутствуют технологические решения по 
изготовлению; 

— развитие существующей номенклатуры количественных показателей для выполнения процедуры оценки 
производственной технологичности, которые позволяют учесть специфические подходы к подготовке 
производства для конкретных производственных комплексов. 

При создании методического обеспечения проектной процедуры оценки технологичности изделий 
необходимо учитывать, что оно должно основываться на обширной конструкторско-технологической базе 
данных с учетом ее структуры и взаимосвязей между элементами моделей, используемых при проектировании и 
реализации технологических процессов. Система планирования многономенклатурных технологических 
процессов отвечает этим требованиям и позволяет, наряду с возможностью проведения оценки известных и 
применяемых в производстве количественных показателей, предложить новые [7]. 

В соответствии с принципами, заложенными при создании системы планирования технологических 
процессов, основным критерием эффективности является время работы производственной системы по 
изготовлению всего комплекта изделий. Оно включает все затраты производственного цикла и напрямую связано 
с себестоимостью производства деталей. Базируясь на этом положении, авторами предложен подход, 
позволяющий на основе установления взаимосвязей между элементами, входящими в оценку себестоимости 
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изготовления деталей, и количественными показателями производственной технологичности дать заключение о 
значимости данных показателей для конкретных производственных условий. 

Результаты исследования. Последовательность реализации разработанного подхода включает в себя 
несколько проектных процедур, реализация которых позволяет учесть, как конструктивные особенности 
обрабатываемых деталей, так и организационно-технологические особенности производственной системы, 
учитывающие состав и технологические возможности оборудования, а также специфику программы 
выпускаемых изделий. 

Первоначально на уровне технологической операции была установлена степень влияния показателей 
технологичности посредством оценки возможного влияния на составляющие при расчете штучно- 
калькуляционного времени (Тиш,к.). На рис. 1 представлена детализация структуры (Тиж.х.) для токарной операции, 
используя которую был проведен анализ и установлены возможности влияния производственных показателей 
технологичности на каждое отдельное значение в расчетах (Ти»..). Аналогичные исследования, позволившие 
установить аналитические зависимости между количественными показателями производственной 
технологичности и структурными элементами технологических операций, были выполнены для других групп 
технологического оборудования. 
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пер. раб Ба уст.сн 


Рис. 1. Структурная схема Тит.к. при точении 


Установлены взаимосвязи между показателями технологичности и затратами по всем статьям 
технологической себестоимости обработки заготовки Сон (рис. 2), а также затратами, связанными с 
организационными простоями оборудования Ср... (рис. 3). 
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Рис. 2. Составные части Соп — затраты по всем статьям технологической себестоимости обработки заготовки 
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Рис. 4. Диаграмма Ганта реализации технологических процессов 
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Для оценки влияния количественных производственных показателей, связанных с затратами времени 
оборудования при простое, создана модель прогнозирования потоков событий поступления партий деталей на 
изготовление на определенную операцию и потоков обработанных деталей с использованием аппарата теории 
массового обслуживания. Данный методический подход использован в системе планирования технологических 
процессов [8, 9]. В качестве примера на рис. 4 представлены результаты в виде диаграммы Ганта. Такой подход 
позволяет учесть конструкторско-технологические характеристики деталей, особенности конкретной 
производственной системы и складывающуюся производственную ситуацию. 

На рис. 5 представлена укрупненная схема структуры для определения затрат производства при изготовлении 
деталей, использованная для оценки удельного веса количественных показателей оценки технологичности в 
системе планирования технологических процессов. Выполнен анализ возможности влияния каждого показателя 
на эффективность функционирования всей производственной системы при изготовлении партии заданных 
деталей для конкретных производственных условий. 

Результаты представленного анализа и установленные взаимосвязи между показателями технологичности и 
эффективностью работы механообрабатывающих систем позволили перейти к решению вопроса об 
установлении значимости количественных показателей производственной технологичности. Представленный 
фрагмент (рис. 6) содержит в реляционной форме информацию о вышеперечисленных связях и дополнен 
сведениями об удельном весе элементов затрат (в процентах). Данные получены на основе статистической 
обработки результатов функционирования реальной производственной системы. В случае отсутствия такого рода 
сведений имеется возможность использовать общемашиностроительные или отраслевые. 
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Рис. 5. Укрупненная схема для определения затрат производства 
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Рис. 6. Фрагмент влияния количественных показателей производственной технологичности на составляющие затрат 


Для установления весовой значимости показателей производственной технологичности остается только 
решение задачи, заключающейся в определении их влияния непосредственно на каждый элемент в 
представленной структуре производственных затрат. Решение данной задачи выполнено с привлечением 
аппарата парных сравнений при принятии решений (итеративный процесс Бержа [10]). Данный метод позволил 
определить степень влияния показателей технологичности относительно других показателей применительно ко 
всем структурным элементам производственных затрат. В таблице 1 представлено в качестве примера попарное 


сравнение показателей технологичности относительно основного времени (То). 


Таблица 1 
Попарное сравнение показателей технологичности относительно основного времени (То) 
То (14 %) 
Оцениваемы Весовые 
11| 12 113 | 14 115 | 16 |17 119 | 20 21 22 23 Итого 

е показатели показатели 
И —- 2 10010 1 2 1 2 2 13 0,098485 
12 0 — 110002000 1 1 2 2 9 0,068182 
13 0 1 1011 1010 2 Р 0 0 7 0,05303 
14 1 2 21-11 1 1 1 1 1 1 1 13 0,098485 
15 2 2 2111-11 1 1 1 1 2 2. 16 0,121212 
16 2 2, 1 110-01 1 1 0 0 8 0,060606 
17 2 0 1 1 1 1 - 1 1 1 0 0 9 0,068182 
19 0 2 211 11212 - 1 1 0 0 12 0,090909 
20 0 2 011 1 1 1 1 — 1 1 1 10 0,075758 
21 1 1 011 1 1 1 1 1 — 1 1 10 0,075758 
22 0 0 В | 1 1 — 1 13 0,098485 
23 0 0 21110121212 1 1 1 - 12 0,090909 

Сумма | 132 1 
Промышленная проверка и апробация разработанных моделей выполнялась в условиях 


ООО «ГАЗПРОММАШУ, специализирующегося на серийном выпуске подогревателей газа прямого действия с 
промежуточным теплоносителем и модифицированного ряда станций, регуляторов, блоков фильтров и клапанов, 
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кранов высокого давление. Исходными данными для проведения экспериментов являлись: сформированная и 
заполненная база данных, содержащая информацию по технологическим возможностям оборудования 
(производственное подразделение № 1), информация по реальному состоянию и технико-экономические 
характеристики функционирования участка (таблицы 2—6), программа изготавливаемых деталей (чертежи 


отдельных деталей представлены на рис. 7, 8). 


Таблица 2 
Соотношение составляющих производственных затрат 
(участок №1 ООО «ГАЗПРОММАШ») 
Зобщ. = 5682 (100 %) 
Соп = 4662 руб. (78 %) (Спр.о. = 1020 руб. (22 %) 
< в < < < ее < 
< : . > > 
а © в > 2 о 2“ о = 2 = о 8 2х 
с > я: и > хх 5х р% а] (® 
со м -+ 2 се д © + > 60 са об Ш а а — — 
Е Хх я 5 ^^ 5 Е и = ^^ и з 
$ © © я © 8 © $ 
Таблица 3 
Затраты, связанные с заработной платой рабочего за выполнение операции 
Зоп = 1872 руб. (32 %) 
Тит.к (29 %) Т, (3 %) 
Тит (25 %) Тиз. (4 %) _ 
То (22%) Тобе (2%) _ Тр.з.1 (1,5 %) Тр.э.2 (1 %) Тир (0,5 %) -- 
Тлех.обс 
См. табл.9 _ _ См. табл.8 См. табл. 9 _ _ 
(1,0%) 
Таблица 4 
Затраты, связанные с изготовлением изделия 
Топ (22 %) 
То т, (8 %) Тлех.обс Торг.об 
(14 %) (1%) (2%) 
Тизм 
— Туст (4 %) Тур (3 %) (1 %) = = 
< м 5 > 8. ^^ —^ > з — 
аа че, 1,8 18| 8 | ее 
= Е Е т № Е 58 а 
в | ее | Ш неее 
Таблица 5 
Затраты, связанные с организационной подготовкой 
Тиз. 1 (1,5 %) 
> © = © © > > > ве . > > < Е < 
5 © : 3 = 5 > 8 55 г 3; $5 = ЕС 2.“ 
= а = а = — 5 = = = — =“ 8 — ре 8 — 5 — $5 &— 
о м ео 
Таблица 6 
Затраты, связанные с наладкой станка 
Те. з.2 (1,0%) 
нана || Е | Еа | де |5 | Е | 9 | Ем | Е9 
ив с ие = 
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Результаты расчетов степени относительного влияния количественных показателей производственной 


технологичности деталей на отдельные структурные элементы производственных — затрат 
(участок №1 ООО «ГАЗПРОММАТ») представлены в таблицах 7-10. 
Таблица 7 
Влияние количественных показателей производственной технологичности 
То (основное (технологическое) время на изготовление или обработку единицы продукции) 
Оцениваемые а оо ка ао В Весовые 
показатели Я В О О О И И показатели 
1.1 1212 1 0 0 0 1 2 1 2 2 13 0,098485 
1.2 0-1 0 0 2 0 0 1 1 2 2 9 0,068182 
1.3 01| - 0 1 1 0 0 2 2 0 0 7 0,05303 
1.4 11212 _ 1 1 1 1 1 1 1 1 13 0,098485 
1.5 21212 1 _ 1 1 1 1 1 2 2 16 0,121212 
1.6 21211 1 0 = 0 0 1 1 0 0 8 0,060606 
1.7 21011 1 1 1 — 1 1 1 0 0 9 0,068182 
1.9 01212 1 1 2 2 в 1 1 0 0 12 0,090909 
1.10 01210 1 1 1 1 1 — 1 1 1 10 0,075758 
1.11 1110 1 1 1 1 1 1 = 1 1 10 0,075758 
1.12 01012 1 2 2 2 1 1 1 = 1 13 0,098485 
1.13 01012 1 0 2 2 2 1 1 1 — | 12 0,090909 
Таблица 8 
Влияние количественных показателей производственной технологичности 
Тупр (время на подвод инструмента к заготовке) 
Оцениваемые — < яч мона я © Е Весовые 
показатели м — ттт -- - И. = показатели 
1.1 _ 2 2 0112212 2 2 2 18 0,136364 
1.2 0 — 0 01011 111 1 1 1 1 7 0,05303 
1.3 0 2 - 21111 111 1 2 2 2 15 0,113636 
1.4 2 2 0 10| 1 11 1 1 0 0 9 0,068182 
1.5 1 2 1 21-21211 12 2 2 2 19 0,143939 
1.6 1 1 1 У О ее И ОО 1 1 1 1 10 0,075758 
1.7 0 1 1 10-11 1 1 1 1 9 0,068182 
1.9 0 1 1 1 111 2 2 2 2 14 0,106061 
1.10 0 1 1 11011 110 - 1 2 2 10 0,075758 
1.11 0 1 0 11011 11011 = 2 2 9 0,068182 
1.12 0 1 0 21011 1100 0 = 1 6 0,045455 
1.13 0 1 0 21011 11010 0 1 —- 6 0,045455 
Таблица 9 
Влияние количественных показателей производственной технологичности 
Ти.з.1 (норма времени на организационную подготовку) 
Оцениваемые аа чмо о ааа Е Весовые 
показатели О О В о о № показатели 
1.1 1212 1 1 1121212121212 |2 21 0,135484 
1.2 СО О ОО А ИО О ВО А ОГО А МО А АО ОВО В АО В АО ВОВ ИВО) 1 0,006452 
1.3 012-11 1 1 1 1 1 1 1100 10 0,064516 
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Оцениваемые а чи оо а = ао В Весовые 
показатели ре = показатели 
1.4 1 1 1 00000 3 0,019355 
1.5 1 2. 1 0 - 1212121212121212 20 0,129032 
1.6 1 2 1 210-022121212 1 1 18 0,116129 
1.7 012 1 2001 - 1 1 1 1 1 1 и 0,070968 
1.8 012 1 2010 1 - 1 1 1 1 1 И 0,070968 
1.9 012 1 2010 1 1 _ 1 1 1 1 И 0,070968 
1.10 012 1 2010 1 1 1 - 1 1 1 И 0,070968 
1.11 012 1 2010 1 1 1 1 - 1 1 И 0,070968 
1.12 0220200 1 1 1 1 1 _ 1 13 0,083871 
1.13 0121212 1 1 1 1 1 1 1 1 - 14 0,090323 
Таблица 10 
Влияние количественных показателей производственной технологичности 
Аст (затраты, связанные с использованием технологического оборудования) 
о Весовые 
Оцениваемые показатели 5 ый 2 = 7 7 Е показатели 

1.8 - 1 1 2 1 1 6 0,2 

1.9 1 - 2 1 1 1 6 0,2 
1.10 1 0 0 0 0 1 0,033333 
1.11 0 1 2 _ 1 1 5 0,166667 

1.12 1 1 2 1 - 1 6 0,2 

1.13 1 1 2 1 1 _ 6 0,2 


1.1 — показатель обрабатываемости материала; 

1.2 — показатель сложности конструкции детали; 

1.3 — коэффициент точности и шероховатости поверхностей детали; 

1.4 — показатель унификации конструктивных элементов; 

1.5 — показатель использования материала; 

1.6 — показатель возможности изготовления заданной номенклатуры деталей; 

1.7 — показатель использования технологических возможностей производственных систем; 

1.8 — показатель технологичности детали по однородности технологического оборудования; 

1.9 — показатель прогнозирования уровня загрузки технологического оборудования при обработке заданной 
номенклатуры деталей; 

1.10 — показатель многовариантности принятия решений при проектировании технологического процесса; 

1.11 — показатель многовариантности принятия решений при реализации технологических процессов; 

1.12 — показатель производственной технологичности детали, отражающий возможность соблюдения 
принципа единства баз при разработке технологического процесса в части поверхности детали, являющейся 
основной конструкторской базой; 

1.13 — показатель производственной технологичности детали, отражающий возможность соблюдения 
принципа единства баз при разработке технологического процесса в части поверхностей деталей, являющихся 
вспомогательными конструкторскими базами. 


На основе представленных моделей и ранее известных зависимостей расчета количественных показателей 
была проведена оценка производственной технологичности деталей. При этом было использовано разработанное 
и зарегистрированное авторами программное обеспечение. Таким образом, с учетом информации о реальном 
СОСТОЯНИИ производственной системы проведена корректировка конструкторской документации, номенклатуры, 
последовательности реализации изготовления отдельных групп деталей и планирование технологической 
подготовки производства. Сравнительный анализ результатов расчетов представлен в таблице 11. 
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Таблица 11 
Сравнительный анализ результатов 
№ 1 вариант 2 вариант Эффективность 
Общее время на изготовление изделий 20 
1 наименований с годовой программой 61193,42 ч 53073,35 ч 15 % 
выпуска 16 600 шт. 
2 Расчетное количество оборудования, 23 17 35 % 
участвующее в ТП 
3 Коэффициент загрузки оборудования 0,67 0,72 7% 


Обсуждение и заключение. Результаты представленных теоретических исследований и их апробация в 
условиях реального производства позволили предложить метод оценки производственной технологичности 
деталей. Он позволяет провести комплексную оценку на основе разработанных аналитических зависимостей 
определения весовых коэффициентов, определяющих значимость каждого показателя технологичности с 
позиции эффективности функционирования механообрабатывающей системы. Отличительной особенностью и 
научной новизной работы является учет реально складывающейся производственной ситуации при оценке 
технологичности. Это позволяет использовать данную проектную процедуру не только традиционно на 
начальных стадиях технологической подготовки производства, но и на этапах реализации технологических 
процессов с целью рациональной организации производства. 

Созданные формализованные модели создают основу для полной последовательной автоматизации 
проектных действий при выполнении оценки технологичности изделий и обеспечивают предпосылки для 
создания перспективной интеллектуальной системы прогнозирования рациональности изготовления деталей на 
конкретном производстве и принятия обоснованных организационных и технологических решений. 
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Прогнозирование поведения участников дорожного движения 
в условиях проселочных дорог для беспилотных автомобилей 


С.А. Иванов® ‚Б. Рашид® 
Центр беспилотных технологий университета Иннополис, г. Иннополис, Российская Федерация 
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Аннотация 

Введение. Благодаря модулю прогнозирования траекторий движения динамических объектов беспилотный 
автомобиль способен безопасно двигаться по дорогам общего пользования. Однако все современные методы 
прогнозирования оценивают производительность только в городских условиях и не рассматривают свою 
применимость к домену проселочных дорог. Цель данного исследования заключается в анализе адаптивности 
существующих методов прогнозирования и разработке подхода, который будет демонстрировать лучшую 
производительность при работе в новых условиях. 

Материалы и методы. В качестве решения предлагается использовать нейронную сеть, включающую в себя 
следующие подмодули: графовый кодировщик сцены, мультимодальный декодировщик траекторий, модуль 
фильтрации траекторий. Также предлагается применить адаптированную функцию потерь, которая штрафует 
сеть за генерацию траекторий, выходящих за границы дорожного полотна. Данные элементы задействуют 
распространённые практики решения задачи прогнозирования, а также адаптируют её для домена 
проселочных дорог. 

Результаты исследования. Проанализированы основные отличия и условия работы модуля прогнозирования в 
условиях проселочных дорог. Выполнена симуляция нового домена путем модификации существующих 
наборов данных. Проведено сравнение популярных методов прогнозирования и оценена их применимость к 
новым условиям. Представлен новый, более адаптивный к новому домену, подход. 

Обсуждение и заключение. Проведенное сравнение с другими популярными методами показывает, что 
предложенное авторами решение обеспечивает более точные результаты прогнозирования. Также были 
выявлены недостатки предложенного подхода и описаны возможные пути их устранения. 


Ключевые слова: прогнозирование траекторий, прогнозирование поведения, нейронные сети, беспилотные 
автомобили, искусственный интеллект, автономные автомобили 
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Введение. Последние достижения в области искусственного интеллекта (ИИ) активно внедряются в 
различные сферы деятельности. Одно из таких достижений — беспилотные транспортные средства (БТС). 
Целью текущих исследований является создание алгоритмов, позволяющих БТС безопасно двигаться по 
дорогам общего пользования. Это позволит значительно снизить количество дорожно-транспортных 
происшествий [1]. 

Научное сообщество уже выделило основные модули беспилотного автомобиля. Один из них — система 
прогнозирования будущего поведения участников дорожного движения (агентов) [2]. Четкое понимание того, 
как будет развиваться окружение и в какую сторону будут двигаться динамические объекты (пешеходы, 
автомобили, велосипедисты), крайне необходимо БТС для поиска и использования безопасной и эффективной 
траектории движения. 

Множество научных статей посвящено задаче прогнозирования таких траекторий [3-12]. Однако в 
настоящее время не проводится активных исследований по применению существующих методов вне городских 
условий. А делать это крайне важно, поскольку автономные автомобили будут использоваться и на 
проселочных дорогах тоже [13]. Городские условия сильно структурированы: автомобили в основном следуют 
по полосам движения, а пешеходы движутся по специальным зонам. В этом смысле область проселочных дорог 
является полной противоположностью, а значит, она будет иметь дополнительные сложности при разработке. В 
данной работе сосредоточено внимание именно на этих сложностях: рассматриваются существующие методы 
прогнозирования и их применимость к новым, менее структурированным условиям. 

Цель исследования заключается в следующем: 

— анализ основных отличий и условий работы модуля прогнозирования в условиях проселочных дорог; 

— симуляция менее структурированного домена проселочных дорог путем модификации существующих 
наборов данных; 
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— сравнение современных методов прогнозирования, в том числе на предмет их применимости к новым условиям; 

— описание нового подхода и доказательство его более высокой точности работы в сравнении с другими 
методами прогнозирования. 

Материалы и методы. На первый взгляд, представленный домен может показаться более простой версией 
городских условий из-за того, что проселочные дороги характеризуются меньшим дорожным трафиком. 
Однако отсутствие сложных дорожных развязок, специальных пешеходных зон, большого количества знаков, 
разметки ит. д. делает домен проселочных дорог менее структурированным, то есть меньшее количество 
правил и специфических шаблонов движения увеличивает хаотичность и снижает предсказуемость поведения 
машин и пешеходов. 

Следующие особенности домена проселочных дорог будут оказывать сильное влияние на выбор 
архитектуры модуля прогнозирования: 

— наличие перекрестков. Несомненно, они на проселочных дорогах более простые в сравнении с городскими, 
но в то же время этот факт простоты означает, что модель должна учитывать мультимодальность и оценивать 
вероятность выбора каждого возможного направления движения на перекрестке при подъезде агента к нему; 

— проселочные дороги не имеют разметки полос движения, пешеходных переходов, велосипедных дорожек 
ит.д. Вместо нее НО-карта будет содержать только информацию о границах дорожного полотна. 
Следовательно, этап кодирования сцены должен учитывать эту особенность, чтобы эффективнее описывать 
окружающий контекст; 

— пешеходы и велосипедисты будут двигаться по одной дороге с обычными и беспилотными автомобилями. 
Следовательно, модель должна быть адаптивной для предсказания будущих траекторий движения как 
автомобилей, так и пешеходов/велосипедистов. 

Модуль прогнозирования подразумевает наличие систем распознавания, трекинга и локализации БТС и их 
точную работу. Авторы статьи задействуют датасет Агооуегзе, который хранит необходимые записи работы 
всех систем в удобном виде [14]. 

Датасет состоит из записей дорожных сцен, зафиксированных на улицах Майами и Питтсбурга США). 
Каждая из записей содержит локальную часть карты местности (границы полос, дорог, пешеходные переходы) 
и список всех распознанных агентов, включающий текущую позицию и историю передвижений каждого из 
них. Каждая из записей делится на две части: две секунды истории наблюдений и три последующих секунды, 
для которых делается прогноз (горизонт прогнозирования). Данные будущего движения объектов также 
доступны и используются для вычисления точности методов прогнозирования и тренировки моделей. 

Информация об агентах представлена в дискретном формате. Временной интервал между измерениями 
фиксирован, в этой работе он равен 0,1 секунды (10 Гц). 


Для каждого момента времени! модуль получает историю наблюдений 5? для каждого обнаруженного 
агента 1. История наблюдений состоит из текущего и прошлых состояний агента, где каждое из состояний 5, — 


это 20 позиция в глобальной системе координат. Авторы делают допущение, что информация о высоте 
является избыточной. 

Датасет также предоставляет доступ к НО-карте, которая содержит информацию о границах полос и 
дорожном полотне, пешеходных переходах. Для симуляции домена проселочных дорог выполнена 
модификация датасета таким образом, чтобы исключить из дорожных карт всю информацию, кроме границ 
дорожного полотна О. Тем самым уменьшено количество информации о дорожном контексте и усложнена 
задача прогнозирования. 

Следовательно, контекст сцены представляется как: 


в о НВ. (1) 

где к— это общее количество отслеживаемых агентов на сцене. 
Такой подход подразумевает прогнозирование траектории только для одного агента за одно выполнение, 
поэтому далее 5? трактуется как 5” для упрощения. Чтобы обобщить модель для всех распознанных 


агентов, требуется повторить предложенный подход для всех К агентов на сцене. Агент, для которого в текущий 
момент делается прогноз, считается целевым агентом. 
Для оценки точности работы методов прогнозирования датасет содержит записанные будущие траектории 


5/7 для каждого целевого агента: 
нь (2) 
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где Н обозначает количество следующих временных шагов. В данном случае параметр Н будет равен 30, так 
как горизонт планирования равен трем секундам с частотой дискретизации 10 Гц. 

Домен работы модуля прогнозирования мультимодален, то есть в абсолютно идентичных дорожных 
ситуациях будущее поведение агентов может существенно различаться. Например, автомобиль, 
приближающийся к перекрестку, может продолжить свое движение прямо или выполнить поворот. Чтобы 
принять это во внимание, на выходе модели необходимо генерировать М возможных будущих траекторий и 
М вероятностей выполнения каждой из них. 


Следовательно, цель модуля прогнозирования — создать такую функцию }, которая принимает на вход 
контекст сцены с и генерирует М пар возможных будущих траекторий и их вероятностей: 
Сел ©/ у 
1(с)=[5 5 о оьВынЫОь, (3) 


При этом хотя бы одна генерируемая траектория 5” должна быть максимально приближена к реальной 


траектории $7 › а вероятность её выполнения р должна быть приближена к единице. 

Архитектура модели. Предлагаемый подход подразумевает использование нейронной сети, состоящей ИЗ 
подмодулей кодирования сцены, декодирования и фильтрации траекторий. Архитектура системы 
представлена на рис. 1. 


НР карта с агентами ! ° Кодировщик сцены \ / Декодер траекторий `; ] 
8 } } : } } ОЧЕН 
| ; оо ; ; Выход: : : Фильтрация траекторий _ 
} | : Кодировщик | | м*он+! : А К го М траекторий 
| } | и ь | | ] НИК. у * » | 
, ' р : полилиний © ®@ ® ; : р ] * у р 


| Граф взаимо- оф 


дейстивий } ; К иоьенния Е 


} ; | 

| | . 

} | УесогМе! ! : } > . 
4 : . : М вероятностей 


Рис. 1. Архитектура системы 


За кодирование сцены отвечает адаптированная к новым условиям нейронная сеть, основанная на векторном 
представлении. Данный выбор обусловлен тем, что на проселочных дорогах НО-карта будет содержать 
ограниченное количество информации (только границы дорожного полотна и историю наблюдений 
динамических объектов). Популярные методы представляют контекст дорожной сценыс в формате 
изображения и обрабатывают его с помощью сверточных нейронных сетей, однако векторное кодирование 
позволяет избежать накладных расходов, которые связанны с генерацией изображения [4—5]. 

Представленный кодировщик основан на модели УестюгМе, однако его входной формат данных был 
модифицирован, чтобы получать информацию только о границах дорожного полотна и состояниях агентов [3]. 
Данный кодировщик представляет границы дороги и состояния агентов с помощью ломанных линий, которые 
далее обрабатываются графовой нейронной сетью. Это позволяет кодировать взаимодействие между ломаными 
линиями. Подробности реализации описаны в статье УесогМе! [3]. 

Декодировщик траекторий представляет собой задачу регрессии нескольких возможных траекторий и 
генерирования набора вероятностей. Для решения данной задачи задействована модель многослойного 
персептрона. Реализация декодировщика вдохновлена моделью МТР [4], однако авторы статьи предлагают 
другую формулу вычисления наилучшей траектории т* из множества М траекторий. Также предлагается 
использовать дополнительный механизм, который штрафует модель за прогнозы, выходящие за пределы 
области движения. 

Авторы оригинальной модели МТР предлагают тренировать многослойный персептрон с помощью функции 


потерь, которая представляет сумму р и Та, ‚› где: 
т, =1.(57,5%,} (4) 
М 
ТР — Ра 105 Р» (5) 
т=1 
В данном случае а — это среднеквадратичная ошибка между реальной траекторией $7 И лучшей 


траекторий т* из М генерируемых. 


1. (57.57) 1>(5 -5', } (6) 
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х дж 
где 5: — реальная будущая позиция агента в момент времени 1, а 5; —_ прогнозируемое будущее состояние 


наилучшей траектории 7"*. 
Г 


выполнения лучшей из прогнозируемых траекторий т* до 1 и уменьшает вероятность других траекторий до 0. 
Г является бинарным индикатором, равным 1, если условие с верно, и 0 — в противном случае. 
В оригинальной статье лучшая из прогнозируемых траекторий т* определяется как та, которая имеет 
минимальное значение среднеквадратичной ошибки в сравнении с реальной траекторией: 
т = агетт (57,5; (7) 


т 


функция потерь, основанная на перекрёстной энтропии, которая увеличивает вероятность 
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Авторы статьи предлагают использовать следующую модификацию: 


т = атетт т. (57,5: } (8) 
теЕД 

где А — подмножество генерируемых траекторий, которое имеет схожее финальное направление с реальной 

траекторией 57. 

Идея заключается в том, чтобы при расчете наилучшей траектории т* убрать из рассмотрения траектории, в 
которых финальное направление агента значительно отличается от направления в реальной траектории. Если 
разница направлений отличается меньше, чем некоторый порог), то генерируемая траектория считается 
корректной, то есть тсл. В рассматриваемом случае у=30°. Следовательно, лучшая траектория т* должна 
иметь схожее финальное направление и наименьшее значение функции потерь. 

В данной работе также задействованы предварительные знания о домене (рйог Кпо\е4эе) для достижения 
большей сходимости модели [15]. Поскольку при движении в домене проселочных дорог с НО-карты доступна 
только информация о границах дорожного полотна, вводится дополнительная переменная Даа в функцию 
потерь. Благодаря ей модель будет штрафовать прогнозируемые траектории, которые выходят за пределы 
дороги в случаях, если хотя бы одно состояние $; & О. Модель штрафует только наилучшую траекторию, 
поскольку только для неё можно определить направление уменьшения ошибки путем приближения наилучшей 
из генерируемых траекторий т* к реальной траектории 5”. 


Следовательно, Г, определяется как: 


а 


11а #\" 
И = г" (5 5 | › (9) 


где Г, равен 1, если 5, р, и0 — в обратном случае. 
Финальная функция потерь определяется как: 


В =Г, 
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+: +В-Г., 


и? (10) 
где аи Вр — гиперпараметры нейронной сети, использующиеся для обучения. В данном случае оба этих 
параметра равны 0,5. 

Для фильтрации похожих и дублирующихся траекторий предлагаемый подход использует фильтрацию 
конечного набора траекторий М на финальном этапе. Этот модуль необходим, поскольку в некоторых случаях 
количество возможных траекторий агента может быть меньше М, например, когда автомобиль движется по 
прямой дороге с постоянной скоростью, модель может сгенерировать только одну траекторию: автомобиль 
продолжает двигаться прямо. Однако необходимость генерировать именно М траекторий приведет к тому, что 
все прогнозы будут схожи между собой. 

Предлагаемая фильтрация основана на финальном направлении и позициях состояний 5:: если направление и 
сумма отклонений между состояниями $; реальной и генерируемой траекторий меньше порогового значения од, 
то траектории считаются схожими. Авторы усредняют каждое состояние траекторий и суммируют вероятности 
траекторий р. 

Такой подход был реализован на языке программирования Рушоп на фреймворке глубокого обучения 
РуюгсН. Модель обучалась на видеокарте СеРогсе ВТХ 2080 Т1 в течение 40 эпох, тренировка заняла четыре часа. 

Результаты исследования. Для оценки точности моделей прогнозирования в этом разделе используются 
широко распространённые метрики для задачи прогнозирования траекторий: среднее значение 
отклонения (ауегазе @1зр]асетеп етог, АШЕ), финальное отклонение (йпа! 1зр!асетепё еггог, ЕОЕ) [6], 
М155Кае (МВ) и ОНЙгоа4 гае (ОБ). 

Для мультимодальных случаев с генерацией нескольких траекторий АШЕ и ЕОЕ принимается как 
минимальное АПЕ и ЕПЕ среди М траекторий (траектория с наименьшим значением метрик) [5]. 
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Метрика М133Кае рассчитывается как процент «промазанных» прогнозов. Прогноз считается 
«промазанным», если метрика АРЕ сгенерированной траектории более двух метров. Метрика ОК 
рассчитывается как процент траекторий, в которой хотя бы одно состояние $; выходит за пределы области 


движений Ш. 
Для визуализации контекста сцены с, а также реальной и прогнозируемых будущих траекторий 5/и57 был 


реализован скрипт на языке программирования Ру®оп с использованием библиотеки Маро ШБ. 

В этом разделе сравнивается работа нескольких разных методов в случае неструктурированного домена. В 
сравнении задействованы следующие методы: 

— фильтр Калмана; 

— предложенный кодировщик сцены с генерацией одной единственной траектории (Зе пгаесюгу ори); 

— предложенный кодировщик сцены со сведением задачи к классификации среди предопределенных 
траекторий: по наборам из 64 и 415 предопределенных траекторий (Е1хе4 зе! с1аз;1НсаНоп); 

— предлагаемый подход (Ргорозе4 арргоасВ). 

В таблице 1 представлено сравнение точности методов при работе в неструктурированных условиях. 
Сравнению подвергается несколько методов, включая предлагаемый подход. 


Таблица 1 
Сравнение моделей в неструктурированном домене работы 

Мешоа Модез | АРЕ, | ЕОЕ, | АБЕз | ЕОЕб | МВ2, | МВ2 | ОВ 
Ка|пап Яег 1 3,78 8,05 3,78 8,05 0,89 0,89 5,89 
5шае наесюгу ошри 1 3,12 6,75 3,12 6,75 0,89 0,89 3,26 
Ехед зе с]1аз;1Исайоп 415 3,27 7,00 1,74 3,57 0,84 0,52, 3,61 
Ехед зе с]1азз1Ясаноп 64 2,6 5,63 1,52. 2,91 0,82 0,49 2,58 
Ргорозе4 арргоасв 6 2,36 5,29 1,32 2,55 0,78 0,38 1,84 


Фильтр Калмана. Простейший способ прогнозирования поведения заключается в получении текущего 
состояния объекта (текущая полоса движения, скорость, направление ит. д.) и распространении этого 
состояния на будущие шаги на основе некоторых допущений, например, что автомобиль будет продолжать 
следовать по своей полосе или будет иметь постоянную скорость и/или ускорение. Другим популярным 
методом для таких задач является использование фильтра Калмана [12]. 

Согласно данным таблицы 1, фильтр Калмана работает хуже, чем все представленные методы, основанные 
на нейронных сетях. 

На рис. 2 представлены два случая. В первом случае фильтр Калмана успешно выполняет прогнозирование, 
поскольку транспортное средство движется прямо, без каких-либо поворотов или изменений скорости. Во втором 
случае фильтр Калмана приводит ошибочный прогноз из-за отсутствия знаний о контексте дорожной ситуации. 


У У 
30 
—\ 
м 30 
20 = 
10 = 10|. 3 р 
5 › . 
о ры 
Е о : О 
[72 
5 
3 -10 
[- 0 | 
Е | А ИВ | -20 ООО ВО И 
— -20 10 0 0 2х -20 10 0 10 20 3х 
= 
& Рис. 2. Пример предсказаний с помощью фильтра Калмана. Пунктирные линии — граница дорожного полотна, 
= линии красного цвета — целевой агент с историей наблюдений, голубого — другие агенты, 


зеленого — реальная траектория, желтого — прогнозируемая траектория, красные кресты указывают предсказанные 
состояния за пределами дорожного полотна 
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5тае щесюту ошри. Данный метод подразумевает использование графового кодировщика сцены, 
который идентичен используемому в предложенном подходе. Вывод сети подразумевает генерацию только 
одной траектории. Данная модель тренируется с использованием среднеквадратичной функции потерь. 

Как показано в таблице 1, нейронная сеть даже с генерированием единственной траектории демонстрирует 
лучшие результаты в сравнении с фильтром Калмана. 

На рис.3 представлена визуализация работы данного метода прогнозирования. Изображение слева 
показывает, что модель может успешно предсказывать поворот агента. Изображение справа показывает, что 
генерации одной траектории недостаточно. Нейронная сеть пытается представить оба возможных исхода: 
движение прямо и поворот направо. В результате модель выводит среднее значение двух исходов. 


ь И | а 
30 ю. 30 | 
20 я 20. 
10 101 ®. 
0 о. 
-10 -10 | 
-20 -20 ‚в 


-30 -20 10 0 10 02 3х 350 220 4-0 0 Ю 020 эх 


Рис. 3. Пример генерации одной траектории. Линией красного цвета показан целевой агент с историей наблюдений, 
зеленого — реальная траектория движения, желтого — прогнозируемая траектория 


Ерхе4 зе! Чазяйсаноп. Реализация была вдохновлена методом прогнозирования СоуетМее [5]. Данная модель 
состоит из предлагаемого векторного кодировщика сцены, за которым следует иной декодер траекторий. 
Декодер представляет собой задачу классификации по предопределенному набору траекторий, состоящему из 
физически реализуемых траекторий транспортного средства с достаточным покрытием. Для экспериментов 
были созданы два набора: из 415 и 64 возможных траекторий. Второй набор имеет такое же покрытие, как и 
первый, но обеспечивает меньшую плотность траекторий. Подробная информация о наборах траекторий 
содержится в статье Соуег Ме [5]. 

Визуализация работы представлена на рис.4. Модель классификации успешно справляется с 
мультимодальностью на перекрестках, но в некоторых случаях отсутствие достаточного покрытия набором 
траекторий негативно влияет на результаты. 

ОИ ЗЕ У НИ Р 


20 | { | мы — 


10 10 
0 0 . Н 0,4 
-10 30 . 0,2 
-20 0,0 


| | р | | | -20 11 | д | | | 
-30 -20 -10 0 10 20 30Х -30 -20 -10 0 10 20 30х 
Рис. 4. Пример прогнозирования с помощью классификационной модели. Линии красного цвета — целевой агент с 


историей наблюдений, зеленого — реальная траектория. М прогнозируемых траекторий с разной вероятностью выполнения 
р: представлены с использованием красно-желтых оттенков 


Как показано в таблице 1, данный метод работает точнее, чем генерация одной траектории, но хуже, чем 
предлагаемый подход. Кроме того, увеличение плотности набора траекторий за счет использования набора из 
415 траекторий не улучшило результаты. Авторы связывают это с наличием шума в наборе данных, который 
исходит от системы отслеживания, используемой при сборе данных. 
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Рторозе@ арргоасй. Предлагаемый подход устраняет недостатки всех вышеописанных методов. Это 
мультимодальный способ прогнозирования, который не страдает от ограничений предопределенного набора 
траекторий. 

Кроме того, согласно таблице 1, предлагаемый подход превосходит все остальные методы по всем 
показателям. Как показано на рис. 5, метод успешно захватывает два возможных исхода на перекрестке: 
движение прямо или выполнение поворота. 


У — — —_ — __ И _ У — — _ _ —__ . _ Р1 
ЕТ. т — 1,0 
30 % ы 

30 ъ * 0,8 

20 
0,6 

10. 
0,4 

0 4 
0,2 

-10 1 
-20 ы -20 0,0 


-30 -20 -0 о 10 20 30Х -30 -20 -0 0 10 20 30 Хх 


Рис. 5. Пример прогнозирования с помощью классификационной модели. Линии красного цвета — целевой агент с 
историей наблюдений, зеленого — реальная траектория. М прогнозируемых траекторий с разной вероятностью выполнения 
Р:. представлены с использованием красно-желтых оттенков 


На рис. 6 показан пример фильтрации схожих траекторий в случае единственного возможного исхода. 
Вероятность того, что агент завершит начатый поворот, близка к единице, поскольку он уже находится в 
процессе поворота. Следовательно, в этом случае вероятность выполнения других исходов близка к 0. 
Предлагаемый модуль успешно фильтрует схожие траектории. 
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о. 0 
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-30 -20 -10 о 10 20 30Х -30 -20 -10 о 10 20 30х 
Рис. 6. Эффект фильтрации. Слева представлен весь набор прогнозов, справа — только отфильтрованный набор 


Ограничения. Хотя авторы оригинальной статьи о модели МТР [4] указывают, что их метод решает 
проблему коллапса режима (то4е соПарзе), проведенные авторами данной статьи эксперименты этого не 
подтверждают. Такая проблема все еще возникает в некоторых случаях. Предполагается, что это связанно со 
следующими особенностями: функция потерь не штрафует нейронную сеть за генерацию всех возможных 
траекторий, которые может исполнить целевой агент, пока лучшая из них максимально близко к реальной 
траектории. Но также модель никак не поощряет сеть к прогнозированию разнообразных возможных путей. 
Следовательно, сети выгодно делать несколько похожих прогнозов в одном направлении, в котором она более 
уверена, чем делать по одному прогнозу для каждого возможного пути. 

Одним из возможных решений данной проблемы может быть использование декодера траекторий, 
представленного в моделях ТиТ, ОепзеТитТ [10-11], которые подразумевают генерацию финальных целей на 
первых этапах работы. В данных моделях сначала генерируются все возможные финальные цели для агента, а 
далее генерируются траектории, которые описывают движение от стартовой позиции до каждой из целей. Это 
позволяет отфильтровывать похожие финальные цели на ранних стадиях и предотвращать коллапс режима. 

Обсуждение и заключение. В выполненной работе исследованы современные методы решения задачи 
прогнозирования траектории. Рассмотрена адаптивность методов к неструктурированным дорожным 
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условиям — проселочным дорогам. Выявлена недостаточная точность методов, и предложен новый подход к 
прогнозированию. 

Предлагаемый подход основан на моделях УесюгМе( и МТР, но был адаптирован для домена проселочных 
дорог. Кроме того, был предложен модуль фильтрации траекторий и дополнительный механизм для функции 
потерь, который штрафует траектории за выход за пределы зоны движения. 

Представленное сравнение показывает, что предлагаемый подход превосходит другие популярные методы. 

Были выявлены ограничения подхода МТР: выходные данные по-прежнему имеют тенденцию к 
коллапсу режима. Предложение для дальнейших модификаций заключается в использовании методов, 
которые генерируют финальную цель на ранних этапах прогнозирования и тем самым меньше подвержены 
коллапсу режима. 
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Визуальная когерентность в дополненной реальности 
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Аннотация 

Введение. 2020-е годы ознаменовались появлением нового поколения компьютерных тренажеров с 
применением технологии дополненной реальности. Одно из преимуществ данной технологии — возможность 
безопасного моделирования опасных ситуаций в реальном мире. Необходимым условием использования этого 
преимущества является обеспечение визуальной когерентности сцен дополненной реальности: виртуальные 
объекты должны быть неотличимы от реальных. Все мировые ГТ-лидеры рассматривают дополненную 
реальность как следующую волну радикальных изменений в цифровой среде, поэтому визуальная 
когерентность становится ключевым вопросом для будущего ГТ, а в аэрокосмических приложениях визуальная 
когерентность уже приобрела практическое значение. Примером может служить разработка корпорацией 
Боинг пилотского тренажера с дополненной реальностью (2022). Визуальная когерентность — сложная 
комплексная проблема, одним из аспектов которой является обеспечение корректной колористической 
тонировки виртуальных объектов в сцене дополненной реальности. Цель работы — разработка нового 
метода такой тонировки. 

Материалы и методы. В разработанном методе (названном спектральной трансплантацией) используются 
двумерные спектральные преобразования изображений. 

Результаты исследования. Предложена технология спектральной трансплантации, обеспечивающая прямую 
передачу характеристик цвета, яркости и контраста от реального фона к виртуальным объектам. Разработан 
алгоритм автоматического выбора оптимального вида спектрального преобразования. 

Обсуждение и заключение. Будучи полностью автоматическим процессом без регистрации условий 
освещенности, спектральная трансплантация решает ряд сложных проблем визуальной когерентности. 
Спектральная трансплантация может стать ценным дополнением к другим методам обеспечения визуальной 
когерентности. 


Ключевые слова: компьютерные тренажеры, дополненная реальность, визуальная когерентность 
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Введение. Современные тренажеры фактически по умолчанию подразумевают использование виртуальной 
реальности (ушиа| геашу — УВ). Преимущества такого подхода хорошо известны, поэтому не будем на них 
останавливаться, но отметим ряд существенных и, что важнее, непреодолимых недостатков, обусловленных 
самой природой технологии виртуальной реальности. УК — это цифровая, дискретная технология, в то время 
как реальный мир непрерывен. Поэтому моделирование реального мира в УВ неизбежно связано с ошибками, 
что снижает эффективность обучения. Однако для обучающих систем еще более серьезным негативным 
аспектом является то, что решения человека во многом основаны на подсознательном учете множества деталей 
реальной картины мира. Этот процесс принципиально невозможно воспроизвести с помощью чисто 
компьютерных технологий (например, УК) по двум причинам: мы до сих пор не знаем (и вряд ли когда-либо 
узнаем) каков механизм работы человеческого мозга, последние спекуляции на тему искусственного 
интеллекта только подтверждают это; детали реального мира, учитываемые при принятии решений, 
практически бесконечны по количеству, возникают случайным образом и имеют самую разную природу 
(визуальную, акустическую, тактильную...). 

Появление в 2020-х годах обучающих систем дополненной реальности (аистеше4 геашу — АК) снижает 
остроту этой проблемной ситуации. Примерами могут служить разработка компанией Воешех пилотского 
тренажера дополненной реальности на основе известного проекта Кб АТАК$, которая началась осенью 
2022 года, а также параллельно стартовавший похожий проект британской ВАЕ Зу\етз или тренажер 
авиадиспетчеров из данной статьи. Все информационное богатство окружающего нас мира в АВ представлено в 
явном виде и не требует моделирования. Но для реализации преимуществ АК, связанных с параллельным 
присутствием реальных и виртуальных объектов в сценах, необходимо решить проблему визуальной 
когерентности (у15йа1 собегепсе — УС): виртуальные объекты должны быть неотличимы от реальных. В данной 
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статье предлагается метод решения проблемы визуальной когерентности в рамках проекта по разработке 
системы обучения авиадиспетчеров. 

АК является производной формой УК. АВ сохраняет все возможности УК, но, помимо этого, как гибридная 
технология обладает значительными преимуществами, вытекающими из параллельного сосуществования 
виртуальных и реальных объектов, что и привлекает внимание разработчиков к УС. Более того, 
исследования [1] показывают, что среди негативных психофизиологических последствий использования 
устройств дополненной реальности доминирует оптический дискомфорт, который возникает из-за разницы в 
восприятии реальных и виртуальных объектов в одной и той же сцене из-за отсутствия УС. Лидеры ГТ-отрасли 
рассматривают АК как следующую «большую волну» революционных изменений в цифровой электронике. 
Поэтому проблема УС становится ключевой для [Т в целом, и эти лидеры проявляют растущий интерес к 
методам ее решения [2]. Однако в аэрокосмических приложениях проблема визуальной когерентности уже 
приобрела практическое значение. Авторы столкнулись с проблемой УС при разработке системы обучения 
авиадиспетчеров: стремительный рост интенсивности воздушного движения в аэропортах привел к увеличению 
частоты столкновений самолетов с другими воздушными судами и аэродромным транспортом при наземном 
маневрировании (>50 случаев по всему миру в 2018 году до начала пандемии). Авиадиспетчеры, работающие 
на вышках аэропортов, не всегда готовы адекватно реагировать на подобные нештатные ситуации, что требует 
дополнительной подготовки. Наиболее эффективная форма такого обучения предполагает представление 
диспетчеру ситуации опасного сближения объектов на летном поле, что невозможно с реальными объектами, 
но может быть абсолютно безопасно реализовано в сценах дополненной реальности. В нашем приложении 
аварийные ситуации были безопасно смоделированы с помощью АК на реальном аэродроме, при этом 
используемые виртуальные самолеты должны быть неотличимы от реальных. 

Исчерпывающий обзор известных методов УС можно найти в статье [3]. Согласно авторской 
классификации, все методы УС можно разделить на два основных класса: с измерением параметров освещения 
и с оценкой условий освещения. В первом случае обязательной процедурой является предварительное 
измерение условий освещенности, проводимое с помощью специального оборудования. Эта процедура — 
длительный и трудоемкий процесс и не представляется возможной, если используется заранее полученное 
изображение или видео реального мира. Во втором случае сложность восстановления картины освещения по 
изображениям приводит к допущениям и ограничениям, что делает результаты неоднозначными. Поэтому, 
несмотря на впечатляющие результаты, полученные исследователями с помощью методов упомянутых в 
обзоре [3], уровень УС все еще часто недостаточен, особенно в сценах АК с реальными естественными 
ландшафтами в условиях рассеянного освещения, которые типичны для авиационных приложений. Как 
показывает приведенный ниже обзор публикаций, существует дефицит исследований такого рода. 

Цель данной работы — разработка универсального и автоматического метода, обеспечивающего прямой 
перенос характеристик цвета, яркости и контраста с реального фона на виртуальные объекты без цифрового 
ЗР-моделирования, которое необходимо в существующих подходах УС. Метод основан на математическом 
аппарате двумерных спектральных преобразований, мы назвали его «спектральной трансплантацией». 

Основными результатами данной работы являются: 

— базовая схема для метода спектральной трансплантации, который обеспечивает прямой перенос 
характеристик цвета, яркости и контраста с реального фона на виртуальные объекты. Метод предполагает 
замену части спектра изображения виртуального мира на такую же часть спектра изображения реального мира 
с последующим обратным преобразованием спектра с пересаженной частью; 

— алгоритм автоматического выбора оптимального вида спектрального преобразования для использования в 
спектральной трансплантации. 

Важно отметить, что УС зависит от многих факторов: освещения, теней, цветового тона, взаимных 
отражений, текстуры поверхности, оптических аберраций, конвергенции, аккомодации и пр. Соответственно 
используются различные техники АВ-визуализации. В нашем случае УС обеспечивается только для факторов 
общего освещения и окраски виртуальных объектов в АК. Это является одним из наиболее сложных вопросов 
УС, особенно для сцен на открытом воздухе. Поэтому для достижения полной УС спектральная 
трансплантация должна использоваться в сочетании с другими методами УС. 

Список источников в [3] включает 175 позиций; этот обзор включает практически все подходы к 
достижениям в УС (за исключением последних, основанных на нейронных сетях, рассматриваемых ниже). 
Поэтому здесь кратко опишем некоторые характерные примеры, которые соответствуют упомянутым 
основным классам. 
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Измерение условий освещения 

Используя зонд освещенности с диффузными полосами между зеркальными сферическими квадрантами, 
Дебевек и др. [4] продемонстрировали, как полный динамический цветовой диапазон сцены может быть 
восстановлен из одной экспозиции. На основе изображения, полученного с помощью зонда, интенсивность 
нескольких источников света может быть оценена путем решения простой линейной системы уравнений. 
Результаты были использованы для рендеринга виртуальной диффузной сферы. 

Альхаками и Тасерян [5] оценивают направление падающего света (прямое освещение) реальной сцены 
используя методы компьютерного зрения с помощью 360°- камеры, прикрепленной к АВ-устройству. Система 
моделирует отраженный свет от поверхностей при рендеринге виртуальных объектов. Затем определяются 
параметры теней для каждого виртуального объекта. 

Оценка условий освещения 

Кнорр и Курц [6] предложили схему оценки условий освещения в реальном мире на основе фото 
человеческого лица. Метод основан на обучении модели вида лица на основе базы данных по лицам при 
известном освещении. Авторы затем восстанавливают наиболее правдоподобные условия освещения в 
реальном мире в базисе сферических гармоник для запечатленного лица. 

Стоит упомянуть работу [7], где описана комбинация измерения и оценки освещенности. Авторы измеряют 
отражательные свойства реальных объектов, используя карты глубины и цветные изображения вращающегося 
объекта на поворотном столе с помощью В СВ-О камеры. Форма объекта реконструируется путем интеграции 
изображений глубины объекта, полученных с разных точек обзора. Отражательная способность объекта 
определяется путем оценки параметров модели отражения по восстановленным изображениям формы и цвета. 

Наиболее близкими аналогами предлагаемого метода являются подходы, которые, как и спектральная 
трансплантация не предполагают предварительных измерений освещения и моделирования условий освещения, 
геометрии сцены, отражения поверхности, а также предусматривают автоматическую обработку. 

Среди таких аналогов — методы передачи цвета от изображения к изображению. В работе [8] представлен 
способ автоматического переноса цветовой статистики (средних значений и стандартных отклонений) с 
эталонного изображения на целевое. Дополнительные параметры используются для того, чтобы избежать 
ручной обработки, которая необходима для определения особенностей передачи цвета в случаях, когда 
изображения имеют сильное различие в цветовой палитре. Эти дополнительные параметры объединяют 
дисперсии эталонного и целевого изображений. Авторы статьи утверждают, что, хотя ручное изменение этих 
параметров требуется крайне редко, тем не менее иногда это необходимо. Кроме того, статистическая природа 
метода вызывает вопросы относительно вида и объема статистики. Также неочевидна способность метода 
обрабатывать некоторые типы изображений (содержащие блестящие объекты, тени). 

Сяо и Ма [9] представили алгоритм для решения проблемы достоверности передачи цвета с точки зрения 
деталей сцены и цветов. Авторы рассматривают сохранение цветового градиента как необходимое условие 
достоверности сцены. Они формулируют проблему передачи цвета как задачу оптимизации и решают ее в два 
этапа — согласование гистограмм и оптимизация с сохранением градиента. Предложена метрика для 
объективной оценки эффективности алгоритмов передачи цвета на основе примеров. 

Преимуществами разработанного метода, по сравнению с [8, 9] и их многочисленными аналогами, являются 
его универсальность, полностью автоматический характер и возможность передачи не только цвета, но и всех 
основных характеристик изображения с помощью одной простой процедуры. 

Предложенный метод использует двумерные спектральные преобразования. Различные виды изображений 
оптимально описываются различными типами спектральных преобразований («оптимально» — в смысле 
соответствия визуального восприятия для реальных и виртуальных объектов). Среди активно используемых в 
цифровой обработке изображений со времен появления цифрового телевидения — дискретное преобразование 
Фурье, дискретное косинусное преобразование, преобразование Хадамара, 5-преобразование и преобразование 
Карунена-Лоэва. 

Материалы и методы. Схема метода спектральной трансплантации (версия, использующая преобразование 
Фурье [10]) показана на рис. 1. В качестве входных данных используются кадры реального мира (\оП4 Нате — 
\УЕ) и виртуального мира (ушаа! Нате — УЕ) (рис. 2). Это естественно для «видео» АК (когда реальный мир 
наблюдается через видеокамеру). Для «прозрачной» АК, когда реальный мир наблюдается через прозрачные 
очки, реальные изображения захватываются с помощью камер, расположенных на АЁ-очках. 
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Рис. 1. Схема метода спектральной трансплантации. 
Версия, использующая преобразование Фурье 


6) в) 
Рис. 2. Видеокадры реального мира (\/Р) и виртуального мира (УЕ): 
а — УЕ, аэропорт, облачная погода; 6 — УЕ, аэропорт, солнечная погода; в — УЕ, виртуальный самолет. \Е представляют 
собой небольшие фрагменты (<25 %) изображений, опубликованных на сайтах зудпеуапрой.сот.аи и 65ай.сот 


Целью данного метода является перенос основных характеристик изображения из \\Е в УЕ. Схема работы 
метода очень проста, хотя операции имеют большой вычислительный объем. Метод реализуется в пять 
этапов (рис. 1): 

1) Выделение цветовых (КСВ) каналов для \Е — \\Ег, \!Ее, УМЕЬ и для УЕ — УЕг, УЕю, УЕЬБ. Модель КСВ 
используется из-за ее общности и корреляции между каналами, которая специфична для этой модели. 

2) Вычисление двумерного прямого преобразования Фурье (тес Еоипег бапзФогт — ОЕТ): РЕТСОМЕ), 
РЕТС\Е»), РЕТСМЕЬ), РЕТ(УЕг), РЕТ(УЕ»), РЕТ(УЕБ). Формула РЕТ приведена ниже: 


М-—М-1 | 


1 т Ш 
Х (К.О т >, >, х,.(т‚п)ехр( 2 + м (1) 


гдес = К, С, В — индекс для красного, зеленого и синего каналов цветного изображения; М, М — номера строк 
и столбцов матрицы пикселей преобразованного изображения; К, [ — аргументы пространственной частоты; 
хе(т,п) — значение пикселя с пространственными координатами (тп) в канале с; Х‹(К,1) — комплексные числа. 

3) Этот этап является ключевым. Трансплантация низкочастотной части (1ю\-Недиепсу раг-ГЕР) 
осуществляется между парами спектров УЕ и УЕ для каждого из красного, зеленого и синего каналов. Это 
означает, что ГЕР УЕ заменяется на соответствующую ГЕР УЕ. Идея спектральной трансплантации 
основана на следующем свойстве ОЕТ: общий характер изображения (т.е. цветовой оттенок, яркость, 
контраст) зависит от пространственных частот, содержащихся в Г.ЕР (включая постоянную составляющую) 
его двумерного спектра. 
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Таким образом, пересаживая Г.ЕР УЕ в спектр УЕ, мы переносим основные характеристики изображения из 
УЕ в УЕ. Для этого удобнее использовать центрированную форму двумерного спектра, где постоянная 
составляющая находится в центре матрицы спектральных коэффициентов, а низкочастотные компоненты 
симметрично расположены вокруг постоянной составляющей. В центрированном спектре ГЕР является 
центральной частью матрицы ОЕТ, и ГЕР имеет размер М! х М (М <М, М <М№). Если М! =\М, то обозначение 
для квадратной матрицы Г.ЕР может быть Г.ЕР(012..Е), где 0 — постоянная компонента, Е — номер наибольшей 
пространственной частоты в матрице Г.Р. 

Размер ГЕР для трансплантации зависит от размера преобразованного изображения (этот размер определяет 
спектральное разрешение) и от объема характеристик изображения, которые должны быть заимствованы из 
УЕ. На данном этапе исследований размер ГЕР определяется эмпирически. Например, наилучшие визуальные 
результаты для изображений размером 512х512 пикселей были получены при использовании ГЕР (012345). 

4) Восстановление каналов КОВ для УЕ с помощью двумерного обратного преобразования Фурье (геуегзе 
Еопиег гапзотт — ВЕТ). Одновременно на этом этапе происходит смешение характеристик УЕ и УЕ. В 
результате получаются КСВ-каналы изображения УЕ с основными цветовыми, яркостными и контрастными 
характеристиками \\Е, а также с характеристиками, унаследованными от исходного УЕ. 

5) Восстановление откорректированного цвета УЕ путем слияния ВОВ-каналов, полученных на 
предыдущем этапе, вырезания виртуальных объектов и построения АК-сцены путем наложения вырезанных 
виртуальных объектов на \Е. 

Очевидно, что если этот метод используется для обработки видеопотока УЕ, то нет необходимости в 
вычислениях ПОЕТ и, соответственно, ГЕР для каждого кадра реального мира, поскольку основные 
характеристики изображения меняются только при радикальном изменении записываемой сцены. Такие 
изменения можно легко обнаружить по скачкам среднего значения пикселей. В эти моменты необходимо 
пересчитать спектральное преобразование для ГЕР. 

Поскольку различные виды изображений оптимально описываются различными видами спектральных 
преобразований упомянутыми выше, целесообразно разработать автоматический алгоритм выбора 
оптимального вида преобразования для использования в спектральной трансплантации. 

Мы предлагаем оценивать разницу между зрительным восприятием УЕ и УЕ по среднеквадратичному 


расстоянию // между ГЕР спектров мощности их изображений (для всех цветовых каналов): 
1 М, М 2 
А. = >> [В, &,1)- Ри (&,1] ‚ с=В,С,В, 
М, М, к=01=0 (2) 


где Руи Ри — двумерные спектры мощности соответственно УЕ и УЕ. Например, в случае преобразования 
Фурье, формула для Р имеет вид: 
Му то ШГ 
Р (К, = х (т‚п)ехр(-72щ(-— +—))!. 
(К, = ХУ, ( АИ м) 


т=0 п=0 


(3) 


Мы предлагаем определять оптимальный вид спектрального преобразования близостью векторов Ди 
среднего вектора, вычисленного по критерию минимальной суммы квадратов расстояний между средним 
вектором и векторами / для всех рассматриваемых преобразований. 

Пусть ДКА, АС, АВ) — нормированный вектор расстояния между спектром ГЕР5 УЕ и УЕ для 
преобразования /. Пусть а(ДаК, ДаС, АаВ) — средний вектор, а О; — расстояние между Д} и Да. Тогда сумма 5 
квадратов расстояний от векторов /; всех рассматриваемых преобразований до среднего вектора равна: 


5=Ур =УУ(А,-А,), с= В, С, В. 
о“ 1с р 


(4) 
Координаты Лав, Лас, Лав среднего вектора вычисляются как решение системы дифференциальных 


уравнении с частными производными: 
—=0 =А, С, В 5 
ОЛ ыы о. ( ) 


Выбор оптимального вида спектрального преобразования определяется условием близости: 


шшр.. 
т (6) 

Другим очевидным критерием выбора оптимального вида преобразования является длина векторов Д. 
Однако экстремумы такого критерия могут быть связаны со способностью или неспособностью определенных 
преобразований правильно обнаружить разницу между определенными видами УЕ и УЕ. Поэтому мы 
рассматриваем использование среднего вектора как более надежный способ отбора. 
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Аналогично расчету ОЕТ для УЕ, выбор оптимального вида преобразования производится только один раз 
в начале спектральной трансплантации, если только \Е не меняется резко. 

Результаты исследования. Предложенный метод был протестирован с использованием преобразования 
Фурье без выбора оптимального вида преобразования. УЕ (реальная сцена аэропорта) и УЕ (виртуальная 
модель самолета) имели размер 512х12пикселей и 24-битные цвета. Исследовали два различных условия: 

1) \Е — фото аэропорта в облачную погоду, рис. 2 а; 

2) \Е — фото аэропорта в солнечную погоду, рис. 2 6. 

В обоих случаях УЕ содержал З)-модель самолета, показанную на рис. 2 в. Протестировали трансплантации 
ГЕР(0), ГЕР(ОТ), [ЕР(012), [.ЕР(0123), ГЕР(01234), [.ЕР(012345). Некоторые результаты тестирования показаны 
на рис. 3. Наилучшие визуальные результаты были получены при использовании Г.ЕР(012345). На рис. 3 
изображения после спектральной трансплантации намеренно показаны без других эффектов УС (тени, 
освещение и т.д.), чтобы продемонстрировать чистые результаты этого метода. 


ЕР(0123) Е ЕР(012345) 
.. —_— с 


8 СЕР 


< : я . Р 7 х-я Р- Г с” . ". . ы #.. > - 
С а, 
(0123) #-— ЕР (012345) = 
6) в) 
Рис. 3. АВ-сцены без трансплантации и с трансплантацией Г.ЕР разных размеров. 


а — АВ-сцена, состоящая из УР и УЕ без трансплантации Г.ЕР; 6 — АР-сцена, после трансплантации 
Г.ЕР(0123); в — АВ-сцена после трансплантации ГЕР(012345) 


а) 


Верхний и нижний ряды на рис. 3 соответствуют противоположным условиям для УЕ: светлые и темные 
\УЕ с разными оттенками. Эксперименты с любыми другими УЕ не добавят существенно новой информации, 
поскольку будут иметь условия между теми, что уже представлены на рис. 3. 

Для демонстрации механизма спектральной трансплантации было проведено численное моделирование. На 
рис.4 показана Фурье-трансплантация с использованием небольшой (8х8) пиксельной матрицы, 
представляющей один из цветовых каналов \\Е и УЕ. Такой маленький размер матрицы позволяет наглядно 
проиллюстрировать процедуру трансплантации. В данном примере матрица \!Е может быть ассоциирована с 
изображением с вертикальной градиентной заливкой, а матрица УЕ — с изображением с горизонтальной 
градиентной заливкой. Другим отличием между УЕ и УЕ является диапазон значений пикселей: 8—15 для \Е 
(«более светлое изображение», 8 — постоянная составляющая) и 0-7 для УЕ («более темное изображение»). 
Показана трансплантация ГЕР(01"), где 1' означает часть первой пространственной частотной компоненты 
(используется из-за очень низкого разрешения матрицы 8х8). 3) вид матрицы ВФ после трансплантации 
свидетельствует о передаче свойств от УЕ к УЕ: край поверхности приподнялся, первый пиксель получил 
значение постоянной составляющей УЕ. Этот пример демонстрирует, как в результате спектральной 
трансплантации УЕ начинает приобретать вертикальный градиент и постоянную составляющую. 
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Зр мем оЁ \Е рха ЗБ яле\ оЕУЕ рахе| тах ЗБ ме\м оЁУЕ ташх и ГРЕ (01°) 
ы У» 20 
10 8 => 15 
6 РР 
ь 4 т 7 И 
о 2 И И 
о А рае ВИ та 
Туо-аппепзюпа! \\Е зресат Туо-аппепзюпа! УЕ зресблит УЕ зресилит ув ГЕР (01°) 


ГЕР() ГЕР(01”) гапзр!атаноп 


Рис. 4. Численное моделирование спектральной трансплантации для матриц 8х8 пикселей 


Спектральная трансплантация предоставляет несколько вариантов изменения параметров этой процедуры: 
изменение размера Г.ЕР; выбор отдельных компонентов спектра для трансплантации; использование различных 
коэффициентов трансплантации для разных трансплантируемых компонентов. 

На рис. 5 показан эффект трансплантации с различными параметрами для различных типов виртуальных 
объектов — виртуальных моделей самолетов, отличающихся текстурой поверхности, маркировкой и блеском. 
На рис. 5а и б изображен виртуальный самолет со сложными текстурами, текстовыми обозначениями и 
отражениями виртуальных источников света. На рис.5г‚,д и е — виртуальный самолет с простыми 
контрастными цветами. Части а и г содержат виртуальные объекты без трансплантации; б и д содержат 
виртуальные объекты после трансплантации Г.ЕР(0123); в и е содержат виртуальные объекты после 
трансплантации Г.ЕР(012345). Виртуальные объекты намеренно показаны без других эффектов УС (тени, 
освещение и т.д.), чтобы продемонстрировать чистые результаты метода. 


\ъхххь. 


Е — 


Рис. 5. Сцены с пасмурной \\Е: а, г — АК-сцены, состоящие из \!Е и УЕ без трансплантации Г.ЕР; 6, д — АР-сцены, 
составленные после трансплантации Г.ЕР(0123); в, е — АВ-сцены, составленные после трансплантации Г.ЕР(012345) 


Информатика, вычислительная техника и управление 


187 


БИр://уезыиК-допз.ги 


188 


'А4уапсей Епятеетиз Везеагсй (Коз1оу-оп-Ооп). 2023;23(2):180-190. е155М 2687—1653 


Важно подчеркнуть, что представленные рисунки иллюстрируют возможности настройки 
предложенного метода, а не конечный результат, поскольку он требует настройки на конкретный \Е. 
Демонстрация хорошо выполненных, но неполных результатов, как это часто практикуется в работах по 
УС, не кажется нам корректной. 

Обсуждение и заключение. Основным осложняющим фактором для описанного метода, представленного 
на рис. 1, являются высокие вычислительные затраты. Наиболее перспективным способом решения этой 
проблемы является прямое преобразование параметров ГЕР УЕ в параметры рендеринга УЕ. Это позволяет 
отказаться от громоздких процедур трех ОЕТ- и трех ВЕТ-расчетов на втором и четвертом шагах обработки и 
требует только трех расчетов \!Е ПОЕТ, по одному разу для каждого участка потока \!Е без существенного 
изменения базовых характеристик \Е. Такой подход позволяет обрабатывать поток УЕ в режиме реального 
времени. 

Другой проблемой является выбор оптимального размера ГЕР. По мере увеличения объема используемых 
пространственных частот, они начинают содержать информацию о содержимом УЕ, поэтому ограничение 
размера ГЕР необходимо для устранения эффекта гибридного изображения [11]. Сложность оптимального 
выбора обусловлена тем, что он связан как с размером Г.ЕР, так и с характером изображения. Последние 
достижения в области глубокого обучения позволяют предположить, что новым подходом, связанным с 
визуальной когерентностью посредством спектральной трансплантации, может стать использование 
генеративно-состязательных сетей (сепегануе аЧуегзайа! пебуок — САМ) для передачи реалистичной 
информации об освещении с исходного изображения на целевое изображение тем же образом, как это делают 
САМ для передачи стиля изображения В частности, было бы интересно сравнить работу САМ в случае наборов 
данных, состоящих либо из КОВ-изображений, либо из изображений, представленных в частотной области с 
помощью ПОЕТ. Мы считаем, что использование последнего подхода поможет выбрать оптимальный ГЕР. САМ 
уже широко используются в исследованиях УС [2], как и нейронные сети в целом [12]. 

В дальнейших исследованиях, связанных с темой данной работы, будут рассмотрены следующие вопросы: 

— автоматическое определение оптимального размера ГЕР для трансплантации с заданным объемом 
характеристик, заимствуемых из \!Е; 

— автоматическое определение точных моментов, когда необходимо рассчитать новую ГЕР УР для 
трансплантации при обработке видео \Е в реальном времени (как указывалось выше, это необходимо делать, 
если базовые характеристики \Е радикально меняются); 

— использование того же подхода в обратном направлении (от виртуального к реальному) для применения 
виртуального освещения к реальным сценам (как виртуальное освещение влияет на окружающую среду). 

Как полностью автоматический процесс без измерения освещения, предлагаемый метод спектральной 
трансплантации решает ряд сложных проблем УС. Например, как наилучшим образом выровнять 
характеристики цвета, яркости и контраста между реальными и виртуальными компонентами в сценах АК. Все 
эти задачи решаются посредством одной простой процедуры без моделирования условий освещения, геометрии 
АБ-сцены или ВКПЕ, что исключает неизбежные ошибки моделирования. Предложенный метод может стать 
ценным дополнением к другим инструментам УС. 
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Аннотация 

Введение. Очевидно, что в ближайшее время сохранят актуальность вопросы оснащения движущихся 
робототехнических комплексов элементами автономного управления. Это требует развития моделей 
группового преследования. Отметим, что оптимизация в задачах преследования сводится к построению 
оптимальных траекторий (кратчайшие траектории, траектории с дифференциальными ограничениями, 
показатели расхода топлива). При этом не рассматриваются аспекты автоматизированного распределения по 
целям при групповом преследовании. Для восполнения этого пробела выполнена представленная научная 
работа. Ее результатом должно стать построение модели автоматизированного распределения преследователей 
по целям в групповом преследовании. 

Материалы и методы. Для изучения группового преследования множества целей сформирована матрица. 
Управляющие параметры движения преследователей модифицированы по минимальной кривизне траектории. 
Детально рассмотрены методы погони и сближения. Показаны возможности модификации метода 
параллельного сближения. Матричное моделирование задействовали для построения схемы группового 
преследования множества целей. Перечисленные процессы проиллюстрированы функциями в заданных 
системах координат и анимацией. Как база функций построены блок-схемы фазовых координат преследователя 
на следующем шаге, времени и расстояния достижения преследователем цели. В ряде случаев расположение 
целей и преследователей определено как точки на окружности Аполлония. Матрица сформирована по 
выборкам, соответствующим распределению преследователей по целям. 

Результаты исследования. Рассмотрены девять вариантов погони, параллельного, пропорционального и 
трехточечного сближения на плоскости и в пространстве. Рассчитано максимальное значение времен 
достижений целей. Отмечены случаи, когда вектор скорости преследователя направлен произвольно и в точку 
на окружности Аполлония. Отмечено, что трехточечный метод сближения удобен, если цель движется по 
баллистической траектории. Для модификации метода параллельного сближения на плоскости строится сеть 
параллельных линий. При этом учтены длина дуги линии (которая может быть произвольной формы) и массив 
опорных точек траектории цели. С данными элементами составлено и решено уравнение. На массиве выборок с 
соответственными значениями времен найдено минимальное время, то есть определено оптимальное время 
одновременного группового достижения множества целей. Для унифицированного обращения к библиотеке 
выражен управляющий вектор через однопараметрическое семейство параллельных плоскостей. Сформирована 
библиотека расчетов управляющих векторов. Показан пример применения матричного моделирования к 
групповому преследованию. Представлена схема группового преследования множества целей. Для двух целей и 
трех преследователей рассмотрены шесть выборок, соответствующих распределению преследователей по 
целям. Данные представлены в виде матрицы. По итогам научных изысканий создана и зарегистрирована 
программа для ЭВМ «Модель параллельного сближения на плоскости группы преследователей с 
одновременным достижением цели». 

Обсуждение и заключение. Исследованы методы использования матриц при моделировании группового 
преследования. Показана возможность модификации метода параллельного сближения. Матричное 
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моделирование группового преследования позволяет выстроить его схему для множества целей. Матрица 
распределения преследователей по целям будет генерироваться в каждый момент времени. Методы 
формирования матриц распределения преследователей и целей представляют интерес при проектировании систем 
виртуальной реальности, для задач с моделированием процесса группового преследования, убегания, уклонения. 
Метод динамического программирования при формировании матрицы распределения преследователей по целям 
открывает возможность автоматизации распределения с оптимизацией по заданным параметрам. 


Ключевые слова: алгоритм группового преследования, оптимизация в задачах преследования, 
автоматизированное распределение по целям, матрица достижения преследователями целей, 
автоматизированное принятие решений, автономное управление, параллельное сближение, пропорциональное 
сближение, трехточечный метод сближения, библиотека расчетов управляющих векторов 
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Введение. Алгоритмы преследования изучаются с точки зрения их классической и оптимальной реализации. 
Исследуется их роль в дифференциальных играх преследования. Прикладная сфера готовых решений весьма 
широка, т. к. результаты таких научных изысканий применимы в различных информационных технологиях и 
системах, в частности в поисковых. Безусловно также, что продолжительное время будут актуальны вопросы 
оснащения движущихся робототехнических комплексов элементами автономного управления, что также 
требует качественной реализации рассматриваемых алгоритмов. 

В работах [1-4] исследовалось согласованное поведение группы преследователей и целей. Для общих 
теоретических и практических вопросов в задачах преследования рассматривались работы [5-9]. Наведение 
преследователя на цель анализировалось в свете информации, представленной в [10-13]. 

При всем теоретическом и практическом интересе к данной теме оптимизация в задачах преследования 
ограничивается построением оптимальных траекторий. Предлагаются, например, кратчайшие траектории, 
траектории с дифференциальными ограничениями, показатели расхода топлива. Но не рассматриваются 
аспекты автоматизированного распределения по целям при групповом преследовании. Для восполнения этого 
пробела проводилась данная научная работа. Ее главным результатом представляется построение модели 
автоматизированного распределения преследователей по целям в групповом преследовании. Показано 
формирование матрицы достижения преследователями целей. При назначении преследователям целей 
выполняется перебор всех возможных комбинаций достижения целей и выбирается комбинация минимального 
значения критерия из формируемого множества с максимальным значением. 

Оптимизация группового преследования множества целей является перспективным направлением развития 
такой дисциплины, как оптимальное управление движением в задачах, связанных с автоматизированным 
принятием решений и автономным управлением. 

Материалы и методы. В модели описанного в статье группового преследования цели движутся по 
предопределенным траекториям. Однако эта предопределенность не имеет принципиального значения. 
Преследователи распределяются по целям автоматически, на основе минимаксного решения целевой функции. 
Затем модифицируются управляющие параметры движения преследователей. В данной работе это параметр 
минимальной кривизны траектории. Такой подход позволяет обеспечить одновременное достижения целей. 

Рассмотрим групповое преследование множества целей: М преследователей догоняют М целей. Сформируем 
матрицу распределения преследователей по целям: 

Ч, где1 =1..№,} =1..М. 

Каждая ячейка ЧФ; содержит информацию о фазовых координатах 1-го преследователя и /-й цели. В матрице 
$; содержится информация о том, каким методом 1-й преследователь преследует /-ю цель. 

Данные, хранящиеся в ячейках матрицы, определяют обращение к библиотеке расчетов управляющих 
векторов преследователя. 

В каждой ячейке матрицы ‘Ч; может рассчитываться прогнозируемое время достижения 1-м 
преследователем /-й цели: #. 

Результаты работы 


В каждой полученной выборке А, = следует найти максимальное значение времен 


ПЕК ИИ кк } 


достижений 1, = Мах{1,| ‚ например из {1, ,(,.,(›,Ё.} (таблица 1). 
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Таблица 1 
Выборки, соответствующие распределению преследователей по целям 
= Цели 
5 2 2 2 1 2 1 1 
В 1 х х х х х 
5 2 х х х х х 
5. 
= 3 х х х х х х 
Выборки А: А> Аз Аз Аз Ав 


Необходимо сформировать матрицы Фу, где {= 1...3, ]=1...2 в соответствии с возможными выборками 
АьКк = 1...‹6. Затем, после обращения находится максимальное значение 1=Мах{1,}. Расчет позволил 
установить, что наибольшее время достижения демонстрирует преследователь Р1, настигающий цель Т! из 
выборки А2. 

Итак, рассмотрим выборку Ах. Можно увеличить до значения параметра 1х все значения 1, зависящие от 
векторов скоростей преследователей и целей, а также их допустимых угловых скоростей. Этим определяется 
максимальное значение #. 

Получив массив выборок {Ах} с соответственными значениями времен {1;}, следует найти минимальное 
время {„„=Мт{1}. Так определяется оптимальное время одновременного группового достижения множества 
целей. 

Алгоритмы расчета следующего шага преследователя и оценки времени достижения преследователем 
цели. На рис. 1 представлен алгоритм функции расчета следующего шага и вектора скорости движения 


преследователя. 


Обращение к библиотеке расчетов 
Начало о 
управляющих векторов И 


Ввод текущих координат преследователя Р › вектора текущей скорости 


преследователя И „ ‚ допустимого углового вращения преследователя Ор, 


управляющего вектора И ‚ дискретного промежутка времени Л! . 


Перевод в локальную систему координат. Абсцисса — вектор И. Начало в точке 


нахождения преследователя Р. Вектор скорости преследователя , . В локальной 
системе координат вычисляется угол & между направлением скорости и, и 


управляющим вектором И (абсциссой локальной системы координат). 


[0 


пей 


Г И... Ч, : ы 
Тем в = р, | ы а 
: эт(а„. ) 


Р. =Р+УИ, № 


пем Нем 


м 


в традиционную Перевод значений Р г 


рен пем’ 


Перевод значений Р в традиционную 


пем’ 


систему координат. систему координат. 


Рис. 1. Блок-схема расчета фазовых координат преследователя на следующем шаге 
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На рис. 2 представлен алгоритм расчета времени и расстояния Достижения преследователем цели. 
Переменная & — это пороговое значение расстояния от преследователя до цели, при достижении которого цель 
считается достигнутой. 


Начало 


Вычисление расстояния между Ввод начальных фазовых координат 
преследователем и целью Д5, Г = 0. преследователя и цели. 


Обращение к функции 
расчета фазовых 


Вычисление расстояния между координат 


преследователем и целью Д5,/=1+1. преследователя. 


Ввод текущих фазовых 
координат цели. 


Рис. 2. Блок-схема функции расчета времени и расстояния достижения преследователем цели 


Если цель движется по предопределенной траектории, то алгоритм, представленный на рис. 2, может дать 
оценку времени #; достижения {1-м преследователем /-й цели. В этом случае выходным параметром функции 
может служить количество итераций процесса преследования №. Количество итераций № — выходной 
параметр функции расчета времени и расстояния достижения преследователем цели. 

Если цель предпринимает ответные шаги, чтобы избежать достижения, следует оценивать время иначе. 
Нужно строить прогнозируемые траектории как составные из сегментов прямых, дуг окружностей, квадратных 
и кубических парабол и других известных линий. Это позволит не решать краевые задачи в расчетном цикле. 

Формирование библиотеки расчетов управляющих векторов. Матрица распределения Чу, где 1= 1...М№, 
]ЕТ...М преследователей по целям будет строиться на каждом дискретном промежутке времени. В каждой 
ячейке матрицы Фу; будет храниться информация о методе преследования. На ней основывается обращение к 
библиотеке функций расчета управляющих векторов 1 (рис. 3-11). 


Т-Р 


1 


Здесь Т. — точка положения цели, 


#-| | 


Рис. 3. Метод погони на плоскости и в пространстве: и, = 


РБ. — точка положения преследователя [14] 
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К-Р 


К-Р . Т — точка положения цели, Р — точка положения 
К-Р 


Рис. 4. Метод параллельного сближения на плоскости: и = 


преследователя, К — точка на окружности Аполлония. 


Она однозначно определена точками Р, Т и вектором скорости цели У, [15] 


Рис. 5. Метод параллельного сближения в пространстве: и = ———— , Т — точка положения цели, Р — точка положения 


преследователя, К — точка на окружности Аполлония. Окружность Аполлония лежит в плоскости У, 


образованной точками Р, Т и вектором скорости цели Ит[16]. Показан случай, когда вектор скорости преследователя 
направлен произвольно. По истечении некоторого времени вектор скорости преследователя направлен в точку на 
окружности Аполлония 


‚где Р, — результат пересечения поверхности <=} (х,у), плоскости 


Ри и 
Ри 1 


Рис. 6. Метод погони на поверхности: И, = | 


9 |-> 


РРТ и сферы 5; с центром в точке Р; Радиус |У |. Л! . Точка Р’ — ортогональная проекция точки Р; на плоскость ХУ. 
и р р Ус |". : р р 
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Рис. 7. Метод параллельного сближения на поверхности: и, = ————— . Здесь Рыи — результат пересечения поверхности 


1+1 


5=/(%,У), плоскости Р.В. т и сферы 5:с центром в точке Р:. Радиус |. "ДЕ. 


Точка Р° 


и 


— ортогональная проекция точки Р.1 на плоскость ХУ. Для унифицированного обращения к библиотеке 


необходимо выразить управляющий вектор. Ф: — однопараметрическое семейство параллельных плоскостей [18] 


Х 
па ЕТ! 
Рис. 8. Метод пропорционального сближения: ЧК г. поет И о | в Е , 
"РЕ 
ЕН. 
, У, - ДЕ. с05| К: агссох 
Их. 2]. 
Д0 = К. агссо5| — - — › 
2.71. и т 
У, - ДЕ-5т| К-агссо5| — — — 
2. Т я Т. 
ре , 
у 
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Рис. 9. Трехточечный метод сближения: 
„=(1-т).О+т-Т,. |Р.-Е| 


(р 
в 


Метод удобен, если цель движется по баллистической траектории 
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Рис. 10. Модификация метода параллельного сближения на плоскости. Строится сеть параллельных линий 
1. (5) =1(5)+Т. -Т, › где 5 — длина дуги линии, Т:— массив опорных точек траектории цели. Решение уравнения 


2 2 > * 
а (5) - Р) = (У, - АЕ) относительно параметра 5 позволит найти значение 5", которое будет соответствовать 


ны а Р. -Р 
В (5 р и = РР В семействе }($) могут быть линии произвольной конфигурации[19] 


1+1 


Рис. 11. Модификация метода погоней на плоскости. Строится сеть ($), где 5 — длина дуги линии, 
Т; — массив опорных точек траектории цели. Выполняется условие, что конец линии }($) проходит через точку Ть, а точка 


2 2 
Р:инцидентна линии }($), то есть используется в качестве лекала. Решение уравнения ( рр (5) _ Р) = (У, - 1)’ 


В. эВ 


= 1+1 1 


относительно параметра 5 позволит найти значение 5", которое будет соответствовать т = Ха (5) ‚и, = | . 
Р-Р 
1+1 7 


В семействе }(5) могут быть линии произвольной конфигурации [20] 


Итак, библиотека расчетов управляющих векторов содержит методы погони на плоскости, в пространстве и 
на поверхности. Методы параллельного сближения рассчитываются на плоскости, в пространстве и на 
поверхности. Методы пропорционального сближения рассчитываются на плоскости и в пространстве. 
Трехточечные методы рассчитываются на плоскости и в пространстве. Модифицированные методы погони 
рассчитываются на плоскости и в пространстве, когда для управления преследователем задействуют 
допустимую кривизну траекторий. Модифицированные методы параллельного сближения рассчитываются на 
плоскости и в пространстве. 

Модификация методов параллельного сближения и погони позволяет строить сеть прогнозируемых 
траекторий, которые допускают различные краевые условия. Это иллюстрируют рис. 3-11. Но на них 
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представлены далеко не все методы расчета управляющих векторов. Подразумевается, что это открытая, 
пополняемая библиотека функций. 

Пример применения матричного моделирования к групповому преследованию. Рассмотрим пример 
группового преследования (рис. 12). 


Г. (5) 


Рис. 12. Схема группового преследования множества целей 


Здесь все преследователи достигают цели по модифицированному методу параллельного сближения, что 
соответствует рис. 10. В модели преследования на. рис. 10 кривизна траектории не должна быть больше 
определенного значения. Поэтому у преследователей Р› и Рз увеличивается первоначальный радиус кривизны 
траектории, что показано на рис. 12. 

Сформирована выборка Ах, в которой преследователь Р; догоняет Т;. Далее происходит первичная оценка 
времени достижения /;. Для оценки времени 1; вычисляются: 

— длина прямолинейного участка до цели, 

— длина дуги сопрягаемой окружности допустимого радиуса. 

Затем выбирается максимальное значение &=Мах{1;}. Увеличение времени 1; до кв данной модели происходит 
за счет увеличения у преследователя Р; радиуса сопрягаемой окружности со значения и; до значения и; +67. 

Рис. 13 дополнен анимированным изображением, где показан процесс группового преследования множества 
целей [21]. 


Траектории преследователей Положения преследователей 


-30 -24 #/-18 -12 -6 0 6 12 18 24 30 
Сети прогнозируемых траекторий Положения целей 
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г Положения целей и преследователей 
Траектории преследователей 


совпадают преследователей 


Сети прогнозируемых траекторий 6) Положения целей 


Рис. 13. Схемы фаз группового преследования: а — начальная фаза; 6 — конечная фаза 


По результатам исследования создана и зарегистрирована программа для ЭВМ [22], которая реализует 
алгоритм группового преследования нескольких целей. Это софтверное решение называется «Модель 
параллельного сближения на плоскости группы преследователей с одновременным достижением цели». 

Обсуждение и заключение. Описаны и визуализированы в виде функций на плоскости, на поверхности и в 
пространстве методы погони, параллельного, пропорционального и трехточечного сближения. Кроме того, 
показаны возможности модификации метода параллельного сближения на плоскости. С применением матричного 
моделирования к групповому преследованию выстраивается схема группового преследования множества целей. 
Отдельно показаны начальная и конечная фаза этого процесса. Расчет времени достижения позволил выявить 
преследователя, которому нужно наибольшее время, чтобы настичь цель из рассматриваемой выборки. 

Итак, предполагается, что матрица распределения преследователей по целям будет генерироваться в каждый 
момент времени. Цели и преследователи могут исчезать, могут появляться новые. Такая матрица может 
использоваться и стороной, представляющей цели, которая уклоняется от преследования. Итоги описанных в 
статье научных изысканий позволят сформировать принципы автоматизированного распределения 
преследователей по целям на основе выбранной целевой функции. Предлагаются алгоритмы модификации 
траекторий преследователей для достижения целей одновременно или согласно установленному графику. 
Также рассматриваются вопросы формирования библиотеки методов преследования. Метод формирования 
матрицы распределения преследователей по целям может быть востребован при проектировании систем 
виртуальной реальности для игровых задач, в которых моделируется процесс группового преследования, 


убегания и уклонения. 
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Размещение нескольких виртуальных объектов в физическом 
пространстве в приложениях дополненной реальности 


М.В. Алпатова® ‚ Ю.В. Рудяк® 
Московский политехнический университет, г. Москва, Российская Федерация 
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Аннотация 

Введение. Проблемы, связанные с размещением виртуальных объектов в реальной среде, существенно 
ограничивают возможности технологии дополненной реальности (АК). Такая ситуация выявляет пробел в 
научных знаниях, требующий дополнительного исследования. Поэтому основной задачей данного 
исследования явилась разработка метода оптимального размещения виртуальных объектов, при котором 
происходит минимизация целевой функции комфортности. Такой подход направлен на усовершенствование 
систем АК и развитие соответствующей теории. 

Материалы и методы. Проведенное исследование основывается на анализе размещения виртуальных объектов в 
АВ/УК приложениях с особым акцентом на оптимизацию. Было предложено понятие комфортности размещения, 
учитывающее размеры объекта и расстояния до границ свободного пространства по координатам Х, У, 7. 
Результаты исследования. В рамках исследования были получены формулы для оптимального размещения 
объектов с произвольной функцией комфортности. Основным критерием является минимизация разницы 
между уровнями комфортности с разных сторон объекта. Было выявлено, что успешное размещение объектов 
требует учета их размеров и зон комфортности, а также решения системы из п линейных уравнений. 
Обсуждение и заключение. Полученные результаты представляют собой важный вклад в исследование 
проблемы размещения виртуальных объектов в АВ/УВ/МК. Они открывают новые возможности для 
улучшения взаимодействия с пользователями и проведения дальнейших исследований в области 
пространственных вычислений. Возможными направлениями для дальнейшего развития являются 
динамические корректировки и интеграция полученных результатов в различные ХК-сценарии. 


Ключевые слова: дополненная реальность, виртуальные объекты, физическое пространство, рациональное 
размещение, математическая модель, уравнения 
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Введение. Быстрое развитие технологии дополненной реальности (АК) открывает новые возможности в 
различных областях — от развлечений до образования и промышленных приложений [1, 2]. Однако, несмотря 
на значительные успехи, существует множество проблем, которые необходимо преодолеть, особенно в 
контексте размещения нескольких виртуальных объектов в реальной среде. Одна из таких проблем — 
рациональное размещение виртуальных объектов в приложениях дополненной реальности для обеспечения 
оптимального и комфортного пользовательского опыта [3—5]. 

Указанная проблема возникает в связи с необходимостью понимания устройством физического 
пространства. Для эффективного размещения виртуальных объектов в реальном мире приложение должно 
иметь возможность корректно интерпретировать материальное окружение, в котором находится пользователь, 
используя сенсоры и камеры мобильных устройств [6]. 

В данной статье представлен новый подход к определению оптимального размещения виртуальных 
объектов в физическом пространстве. Эта проблема имеет некоторые сходства с другой близкой к ней темой 
генеративного контекстного дополнения сцены (сомехша| зсепе аистетаноп, СЗА), где основная цель 
заключается в создании гармоничного и удобного взаимодействия между виртуальными и физическими 
объектами [7, 8]. Однако предлагаемый авторами подход отличается от упомянутого, так как он фокусируется 
на определении оптимального расстояния между объектами с использованием монотонной функции 
комфортности, в то время как СЗА учитывает семантику сцены, контекст и смысл виртуальных объектов. 

Существующие подходы к решению проблемы оптимального размещения виртуальных объектов обычно 
ограничены предположениями о форме функции комфортности и не всегда гарантируют оптимальность 
решения. В данной статье предлагается новый подход, который значительно отличается от используемых в 
настоящее время. Это позволяет более гибко находить рациональное расположение виртуальных объектов и 
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предоставляет универсальное решение для различных сценариев и условий. Термин «рациональное 
размещение» применяется по следующим соображениям: во-первых, окончательный выбор при предлагаемых 
рекомендациях все же остаётся на усмотрение пользователя; во-вторых, несмотря на то, что в контексте данной 
работы решается задача минимизации целевой функции, она содержит элементы нечётких множеств. 

В рамках данной работы вводится понятие оптимального размещения набора виртуальных объектов в 
одномерном физическом пространстве. Разработана модель, позволяющая решать задачу оптимального 
размещения набора объектов и приводящая, независимо от выбора вида функции комфортности, к системе 
линейных уравнений для определения оптимальных расстояний между объектами. В качестве примера 
приведено решение задачи для случая с двумя виртуальными объектами. 

Цель данной статьи предложить и продемонстрировать новый подход к определению оптимального 
размещения виртуальных объектов в физическом пространстве, обосновать его применимость и 
эффективность. Представлены математические формулировки и методы решения задачи оптимального 
размещения, а также примеры практического применения полученных результатов. Это позволит выявить 
новые возможности для улучшения взаимодействия между виртуальными и физическими объектами, а также 
способствовать развитию теории и практики дополненной реальности. 

Таким образом, данная статья направлена на углубление понимания проблемы оптимального размещения 
виртуальных объектов в физическом пространстве и предложение нового подхода к еб решению. Результаты 
исследования могут быть использованы для создания более эффективных и удобных в использовании систем 
виртуальной реальности, а также для дальнейшего развития теории в данной области. 

Материалы и методы. Размещение виртуальных объектов в реальном физическом пространстве — это 
задача, которая возникает практически в каждом АВ/УК приложении. При всей своей простоте она может 
создавать большие проблемы в случае недостаточного внимания к вопросу оптимизации размещения таких 
объектов, вплоть до полного отказа большого количества потенциальных пользователей работать с 
упомянутыми приложениями. Наиболее остро стоит задача оптимизации, когда необходимо разместить сразу 
несколько виртуальных объектов в заданном физическом пространстве. При этом, даже в случае размещения 
только одного объекта, лишь совсем недавно было сформулировано понятие «комфортности» его размещения и 
предложена соответствующая модель [9], состоящая в следующем. 

Объект, встраиваемый в трёхмерную физическую область, представляется в виде прямоугольного 
параллелепипеда с характерными размерами: [ — длина; 4 — ширина; й — высота. При этом для каждой из 
координат Х, 7, У вводится следующее понятие комфортности размещения. Понятно, что размер свободного 
пространства не должен быть меньше размера объекта, но, помимо этого, для каждой координаты вводится 
понятие комфортных расстояний от объекта до границы свободного пространства. Например, для координаты Х 
введём понятие комфортного расстояния слева — ДР и справа — ШО, и, соответственно, левой и правой 


комфортности — К и К,. Обозначим расстояния слева и справа от объекта до границы свободного 


пространства Х_и Х, . Будем считать, что К_ =1, если Х_>) и уменьшается до 0 при приближении р_к 


нулю. Например, для простоты возьмём линейные зависимости К_ (х @ / р.) и К, (х, | р,) : 


Х 
—,хХ <р, 

К =: В | (1) 
х>О 


Зависимость К’, (Хх и. р,) аналогична (1). Точно таким же образом введём понятие комфортности с одной и 


другой стороны для координат Й и У. 
Если размер свободного пространства по горизонтали Г>Р +[+0)., то проблема комфортного 


размещения (К_=К, =1) не возникает, и все проблемы появляются, когда [<Ё<р +1+0,. В этом случае 


вводится понятие комфортности размещения объекта, когда комфортность с одной стороны не достигается за 
счёт комфортности с другой стороны. Введём целевую функцию комфортности: 


к. =(к-кК,). (2) 

Под оптимальным размещением будем понимать такое размещение, при котором достигается минимум К,. 
Очевидно, что это происходит, если К_ =К,. 

Как было показано в работе [9], если зависимости К_ (х_ / р) и К, (х, / р, линейны, минимум целевой 


функции (2), соответствующий оптимальному размещению объекта, достигается при следующих значениях Х_ и Хх, : 
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р 
Х_ =(2-П——, 
- } р-р. 
р 
ХХ, =Ё-Х_-1=16-——. 
: в (3) 


Формулы (3) просты и удобны для оптимального размещения одного виртуального объекта. Как отмечается 
в [9], их недостаток состоит в том, что они были выведены для случая линейных (1) функций комфортности 


К_ (Х. /Б.) и К, (х./р.) . Покажем, что они справедливы всегда, если функции К_ (х_/р.) и К, (х./р.) г 
представляют собой одну и ту же функцию (х), удовлетворяющую условию, что она монотонно возрастает для 


0<х $1, и равна 1 при х >1. Вид такой функции К_ (х. /Р.) =К (х) представлен на рис. 1. 


1,0 


0,5 


Критерий комфортности К 


0 0,5 1,0 1,5 


Отнормированное на зону комфортности расстояние Х- / О- 


1-х, 0<х<1 


Рис. 1. Зависимость функции комфортности К(х) = 
1 х>1 


Ниже на рис.2 и 3 представлены дополнительно ещё две функции — кубическая и линейная, 
демонстрирующие схожее описываемое поведение. 


1,2 


реше) — 
со [=> 


Критерий комфортности К- 
© 
^^ 


0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1 1,2 1,3 14 


Отнормированное на зону комфортности расстояние Х- / О- 


х?, 0<х<1 


Рис. 2. Зависимость функции комфортности К(х) -] | И 
‚х> 
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0,8 


0,6 


Критерий комофртности К- 


0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1 1,2 1,3 1,4 


Отнормированное на зону комфортности расстояние Х- / Р- 


х, 0<х<1 


Рис. 3. Зависимость функции комфортности К(х) = ] ] 
,х> 


Действительно, пусть нам надо встроить объект размером / с комфортными расстояниями слева П_ и справа 
р. в пространство, ограниченное размером Г, и при этом выполнены условия 1<Г<р +1[+),. Тогда, если 
односторонние комфортности представимы в виде К. (Хх _/ р.) =К [% К, (х й /р.) = иж) ‚ где 
х =х. / р, х, =Х, / р, ‚ а зависимость К(х), удовлетворяет вышеназванным условиям, то из условия 
минимума целевой функции (2) получим К_ = К, >К(Х_/Рр_)=К(Х./Р,). 

Для конкретно заданной нелинейной функции комфортности А(х) можно численно решать полученное 
уравнение одним из известных методов. Однако, поскольку односторонние комфортности описываются 
автомодельной функцией К(х), из уравнения к(Х_ /Р.) = к(Х. /р.) получаем соотношение Х_/Др_ =Х, /Р, , 
откуда с учётом равенства [.=Х_-1[+Х, вытекают равенства (3). 


Таким образом, показано, что простые и достаточно удобные равенства (3), позволяющие встраивать объект 
с оптимальной комфортностью справедливы для любой функции односторонней комфортности «(х). 

Результаты исследования. В реальной ситуации возникает необходимость размещения сразу нескольких 
объектов. В этом случае решать последовательно задачи оптимального размещения сначала первого объекта, 
потом второго, третьего и т.д. нет смысла, поскольку при размещении следующего объекта возникает 
необходимость сдвигать ранее установленные объекты, чтобы размещение было оптимальным по совокупности 


п 
всех объектов. Если выполнено условие да ‚ объекты в принципе не умещаются в свободном 
1=1 


п. 
пространстве протяжённости Г.. Если же Ё> увы +19 +00), то объекты можно разместить так, чтобы они 
1=1 


не мешали друг другу. Реально задача оптимального размещения возникает, если выполнены условия: 


У <<» 9 + +09). (4) 
11 = 
В этом случае введём целевую функцию: 
К, =К+К®+..+К®, (5) 
где КО — определяемая формулой (2) комфортность 1-го объекта. При этом, если имеются 2 


последовательных объекта с номерами Ё и #1, то комфортное их соседство будет, если расстояние между ними 
и Е мВ ыы 
не менее, чем ро +0”, что соответствует новым комфортным расстояниям Бо =6%) = ро +09, 


1=1,2,...,(и-1) поскольку комфортное расстояние до стены это одно, а до другого объекта, откуда может быть 
что-то выдвинуто, совсем другое. Рациональное расположение встраиваемых объектов определяется минимумом 
целевой функции (5). 
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Минимум целевой функции (5) с учётом (2) даёт нам систему из и уравнений: 

к =0; КО =0; ..,К®=0. (6) 
Из (6) следует: 
Вок Пера. (1) 


Система уравнений (7) может быть записана в следующем виде: 


Хх х@) 
А || 50 | 
р 5 


(у 8 
хе) Г- (х +[ ) (8) 
К —|=А ве 
Бе) Б” 
То есть мы получили п уравнений относительно п неизвестных х0,х@... хе ‚ где хо — расстояние 
первого объекта от левого края области встраивания, ХИ. 1=2,3,...п — расстояние между объектами с 


номерами [и (1-1). 
Поскольку функция К(х) монотонная, система (8) сводится к линейной системе уравнений, которая никак не 
зависит от вида самой функции комфортности ((х). 


Хх 1 х@2) 
9 Бо’ 
- + 
Хх В) хе 
50 БО’ 
[у 
зы 
И И 
Г (”) СБ (п) ь (9) 
Б Б 


Система (9) может решаться одним из известных методов. Однако, ввиду того, что матрица системы (9) 
является сильно разрежена, решение системы можно найти достаточно просто. Например, в первых (п-1) 


| ыы 
уравнениях можно в каждом {-м уравнении выразить хи через Хх, затем, подставив это в последнее 


> 1 
уравнение, получить линеиное уравнение относительно х' ) . После этого, двигаясь от первого уравнения к 


(п-[)-му, можно последовательно найти хо х® х®. 


В случае размещения двух виртуальных объектов, и=2, система (9) принимает следующий вид: 


х( х@) 
ре ро’ 
хе _ Г хо ра [о 12) 5 
БО = 56 0% 
Из системы (10) находим: 
— К 12) \ д) р@ 
ой [р 0 )5 Гр) 
- (59+) 59 +2059” 
. | ее (11) 
(2-ю 12) Боро 
х@) МЕ ы + 
ы (5' +2?) 50 +р®р@) 
где Х® — это расстояние между первым объектом и левым краем пространства; Х @ — расстояние между 


правым краем первого объекта и левым краем второго объекта. 
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Рассмотрим пример (12) когда: 


Г.=100, [® =40, [2 =30, О =10,00 =5, 08 =11,5% =13. (12) 
Поскольку выполнены условия (13), то: 
9+0 <1<Ю0 + + 00 +00 +09 +0 (13) 


Для оптимальной расстановки двух объектов можно воспользоваться выражениями (11), которые дают 
независимо от вида К(х) в этом случае Х 0) =7, 692; Х®) =12,308; ко = ко =0. При этом значение 
односторонней комфортности зависит от вида К(х). Для линейной зависимости К(х), изображённой на рис. 3, 
ко =К®=К в =К и =0,769. Если зависимость А(х) соответствует рис. 1, получим значение комфортности 


1 


С 
№ 
р 


=К@) = К@) =0,973. 


Таким образом, в работе введено понятие оптимального размещения набора виртуальных объектов в 
физическом пространстве. Разработана модель, позволяющая решать задачу оптимального размещения набора 
объектов. Показано, что решение этой задачи не зависит от вида монотонной функции комфортности. 

Обсуждение и заключение. Рассмотрен теоретический аспект важного вопроса оптимального размещения 
набора виртуальных объектов в физическом пространстве — проблемы, часто встречающейся в приложениях 
дополненной реальности. Предложив новую математическую модель и включив нечёткую логику, мы заложили 
основу для алгоритма, который потенциально может помочь пользователям найти рациональное и удобное 
расположение виртуальных объектов в их реальном окружении. 

Основы, заложенные в данном исследовании, показывают, что предложенная модель эффективно решает 
задачи, связанные с размещением множества виртуальных объектов в заданном физическом пространстве. 
Анализируя виртуальные плоскости и учитывая расстояния между виртуальными объектами и краями этих 
виртуальных плоскостей, наш метод обеспечивает рациональное размещение с учётом линейных размеров 
виртуальных объектов и заложенной комфортной зоны вокруг них. 

Полученные результаты способствуют текущему развитию приложений дополненной и смешанной 
реальности, предоставляя теоретическое решение проблемы оптимального размещения, которое, в свою 
очередь, способно улучшить взаимодействие с пользователем и общую удовлетворённость инструментами 
обсуждаемой технологии. Более того, возможные применения приведенного исследования выходят за рамки 
АВ-приложений, открывая новые пути для исследований в смежных областях, таких как виртуальная 
реальность, смешанная реальность и пространственные вычисления. 

В свете полученных результатов данного исследования, будущие разработки могут быть направлены на 
проверку алгоритма путём эмпирического тестирования, включение динамических корректировок в реальном 
времени на основе поведения пользователя и изучение интеграции нашего подхода в различные сценарии ХВ- 
приложений. Поскольку область дополненной реальности продолжает развиваться, мы ожидаем, что наше 
исследование внесёт значительный вклад в развитие технологии, способствуя её широкому распространению и 
дальнейшему обогащению пользовательского опыта. 
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Аннотация 

Введение. В области вычислительной математики известно множество способов аппроксимации модели 
механики жидкости. Учеными выработаны методы и оценки критериев качества аппроксимации, таких как 
устойчивость и сходимость. Комбинация подходов построения экономичных разностных схем, таких как 
расщепление по физическим процессам, регуляризация по Б. Н. Четверушкину, линейная комбинация 
разностной схемы «кабаре» и «крест» в совокупности ранее не реализовывалась и не оценивалась. Перед 
авторами стояла задача аппроксимировать каждую часть расщеплённой по физическим процессам модели 
гидродинамики наиболее адекватной схемой и далее исследовать корректность данного подхода. 

Материалы и методы. Математическая модель гидрофизических процессов замыкается эмпирическим 
уравнением состояния соленой воды. Выбираются значимые свойства, строится математическая модель. 
Разностные операторы аппроксимируют дифференциальные операторы. Строится алгоритм послойного 
моделирования переходных процессов. Алгоритм реализован в виде программы, которая, в основном, содержит 
поэлементные (массивно параллельные) операции. 

Результаты исследования. Получены математические модели гидродинамических процессов в водоемах, 
учитывающие три уравнения движения при наличии градиента плотности водной среды при отказе от 
гидростатического приближения. Апробирован новый способ вычисления поля давления с применением 
регуляризаторов по Б.Н. Четверушкину в уравнении неразрывности. Разработан программный модуль 
численного моделирования гидрофизических процессов движения воды с различной солёностью и плотностью. 
Это открытое программное обеспечение, допускающее не только переопределение эмпирических зависимостей 
(как алгебраических функций), но и подключение внешних моделирующих модулей для отображения 
зависимостей алгоритмически. 

Обсуждение и заключение. Разработанная модель гидрофизики, учитывающая свойства солёной воды и 
динамическую связь механического движения воды с солёностью, может применяться для изучения 
формирования неравновесного распределения параметров и идентификации наиболее стабильных параметров 
водной среды. Модель объясняет нисходящее движение кислорода, что позволит в будущем оценивать 
величины параметров водной среды, которые сложно измерить непосредственно. Она может быть использована 
в процедуре параметрической идентификации трудноизмеряемых параметров водной среды. 


Ключевые слова: математическая модель, стратификация, модель динамики морской воды, 
квазигидродинамическая модель, метод заполненности ячеек, центрально-разностная схема, метод прогонки, 
ВВЦИ, Кабаре, Крест 
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Введение. Одной из важных задач, связанных с экологией и безопасностью жизнедеятельности людей, 
живущих в прибрежных территориях, является прогнозирование и моделирование движения воды в морях и 
крупных региональных водоёмах. Кроме того, актуально изучение процесса переноса растворенных в водной 
среде веществ при учёте стратификации и зависимости плотности воды от многих переменных факторов. Такое 
прогнозное моделирование, вероятно, позволит не только оценивать качество вод, но и управлять им в 
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условиях климатических изменений и индустриальных воздействий. Более общая и главная цель — связать 
качество вод с численностью и видовым разнообразим естественных гидробиоценозов, обитающих в 
гидросфере. При решении этих задач необходимо учитывать гидродинамические характеристики водной среды, 
особенности внешних факторов воздействия, таких как неоднородность распределения температуры, соленость, 
насыщенность воды кислородом, количество растворенных в воде газов, кислотность. Эти параметры являются 
одними из ключевых параметров биологической активности водной экосистемы [1]. 

Акватория водоема также может рассматриваться как передаточная система, проводящая кислород от 
атмосферы к донным осадкам. Однако известны случаи образования химических градиентов большой 
величины на относительно небольших перепадах глубин — в пограничных слоях при относительно небольшой 
их величине за пределами этих зон. Причина появления таких участков при общей стратификации воды по 
плотности во всемирном тяготении с одной стороны и радиационном воздействии солнца на воду, приводящего 
к её нагреванию, с другой. Эти процессы могут привести к снижению скорости продукции, деструкции и 
рециклинга биогенных элементов и биоорганизмов вплоть до остановки этих процессов, а также 
предопределять биоразнообразие гидробионтов в целом и видовой состав в частности [2]. Температурная 
стратификация существенно воздействует на распределение организмов в толще воды, на перенос и осаждение 
вредных для биоорганизмов примесей. Повышение температуры поверхностных вод приводит к нарушению 
вертикального водообмена и, соответственно, к уменьшению аэрации глубоководной зоны, снижению 
растворимости и концентрации кислорода в воде. Стратификация по плотности, температуре и химическому 
составу ограничивает конвективный подъем в поверхностные слои воды биогенных элементов, углекислоты и 
продуктов неполного окисления органических веществ, поступающих в гиполимнион (холодные, солёные, 
плотные слои воды) в результате седиментации (оседания) сестона. С начала стратификации и до момента ее 
окончания поверхностный слой обедняется, а гиполимнион, напротив, обогащается этими веществами. В 
результате физико-химическая стратификация приводит к возникновению неравномерного распределения по 
глубине ряда биологически значимых веществ и является причиной самоорганизации сложной структуры 
экологических ниш [3]. Математическое моделирование механических, химических и биологических 
процессов, протекающих в водных экосистемах, чрезвычайно актуально и связано с проблемами экологии и 
безопасности жизнедеятельности населения прибрежных территорий. 

Среди выдающихся ученых, внесших значительный вклад в изучение гидрологии и океанологии, следует 
отметить В.П. Дымникова, который занимался исследованием климата и океанологии, моделированием 
атмосферы и океана. А.С. Монин и М.Ю. Белевич исследовали процессы кинематики водной среды, 
турбулентности и микроструктуры океана. Советский ученый А.С. Судольский занимался исследованием 
динамики вод и береговых процессов в различных водоемах применительно к решению задач проектирования, 
строительства и эксплуатации конкретных гидротехнических сооружений или рационального хозяйственного 
использования водоемов в целом [4]. По мнению В.И. Вернадского, одним из важнейших проявлений жизни 
является газовый обмен организмов с окружающей средой, главным образом — процессы дыхания, основанные 
на потреблении кислорода [5]. Изучением гидробиологических процессов водоемов занимались различные 
выдающиеся отечественные ученые, среди которых С.В. Бруевич, работы которого посвящены разработке 
аналитических методов исследований, формулировании основ гидро- и биогидрохимии. Г.Г. Матишов и 
В.Г. Ильичев активно изучают условия оптимальной эксплуатации водных ресурсов, занимаются разработкой 
моделей транспорта загрязняющих веществ в водоемах и исследованием оценки их воздействия на биоресурсы 
водной среды [6]. 

Для изучения влияния перечисленных процессов водной среды разрабатывается комплекс взаимосвязанных 
математических моделей, основанный на использовании точных прогностических моделей и программной 
реализации экономичных численных методов, который позволит детально исследовать кинематику процесса, 
причинно-следственные связи и состояние объекта моделирования. Существующие методы и средства 
прогнозного моделирования состояния водной среды с учетом ряда биотических и абиотических факторов, 
включая процессы распределения кислорода, углекислого газа, солей, основаны на общих научных подходах, 
упрощенных математических моделях, обладающих низкой адаптивностью, отсутствием возможности 
моделирования нелинейных динамических процессов, свойственных большинству водных экосистем, 
некорректным заданием границ расчетной области. В некоторых случаях исследования сопровождаются 
формальным определением граничных условий, дают достаточно грубые и приближенные результаты 
моделирования. 

При моделировании процесса переноса веществ, основанного на уравнениях адвекции-диффузии, 
необходима хорошая аппроксимация адвективных слагаемых, представляющих собой градиенты давления, 


Кудинов Н.В. и др. Моделирование вертикальных движений морской воды в стратифицированных водоемах 


плотности массы и полной энергии, импульса движения. Применение стандартных разностных схем при 
завышенных оценках параметров подобия приводит к потере точности вычислений из-за возрастания 
погрешности аппроксимации и усилению ограничений на шаг по времени в связи с условием устойчивости 
разностной схемы. В работах А.И. Сухинова, А.Е. Чистякова и др. [7, 8] показано, как эффективно использовать 
линейную комбинацию разностных схемам «кабаре» и «крест» с оптимальными значениями весовых 
коэффициентов для аппроксимации уравнения переноса. Эффективность этих методов достигается 
оптимизацией погрешности аппроксимации дискретной моделью сплошной среды точного решения задачи 
переноса вещества с постоянной скоростью. Исследования показали, что такой подход распространяется и на 
модели гидродинамики с переменной (знакопеременной) скоростью без эффекта сеточной вязкости. Ещё одно 
положительное качество такого рода аппроксимаций состоит в том, что с их помощью можно моделировать 
сложные структуры течения, к примеру, вихревые. В настоящее время многие исследователи используют 
подобные схемы для моделирования турбулентных течений. Сотрудниками ведущих зарубежных научно- 
исследовательских организаций, таких как З‘апЮга ОшуегзИу, Ппрепа| СоШезе Гоп4доп и др., а также 
сотрудниками Института вычислительной математики РАН Е.М. Володиным, 
А.В. Глазуновым, А.С. Грицуном, Н.Г. Яковлевым и др. [9] изданы труды, в которых математическое 
моделирование климатических изменений, гидродинамических и атмосферных процессов и явлений 
осуществлено на основе вихреразрешающих схем. Ещё больше снизить требование к шагу по времени и 
повысить пространственное разрешение модели при ограниченной компьютерной памяти позволяет 
квазигидродинамическое приближение сплошной среды. На практике в систему уравнений Навье-Стокса и 
неразрывности добавляется малое слагаемое, пропорциональное второй производной по времени от функции 
плотности. Данный подход позволяет сгладить нефизические флуктуации плотности массы и импульса, а также 
полной энергии, переносимые по пространственной сетке быстрее скорости звука. 

Существующие универсальные пакеты прикладных программ (например, пакет программ «Магз3а», 
Экоинтегратор, программный комплекс СНАВ!ЗМА, комплекс ЗАГМО, комплекс программ СНТОМ, 
САКОПМАГ, пакеты моделирования различных процессов аэрогидродинамики, программные комплексы 
РНОЕМ!С$, ЕГОЕМТ, СА РУМАМ!С$ ТООГ) не учитывают некоторые свойства моделируемых сложных 
систем, снижая таким образом точность и оперативность моделирования. К таким свойствам относятся: 
пространственная неоднородность движения водной среды, вихревые структуры течений. В математических 
моделях и в алгоритмах их численной реализации не учитывается вероятность значительного изменения 
глубины, плотность водной среды, что может привести к неустойчивости полученных численных решений. По 
этой причине такие специализированные программные пакеты могут быть использованы для моделирования 
ограниченного разнообразия гидрофизических процессов водных систем. Большая часть известного 
специализированного программного обеспечения (АРАМ, САГЗОНС, СВепз1, ТАЗСЙо\, [$С-3, РАМАСНЕ, 
ВЕМ$АО, ОАМ-ГУ, ЭКОЛОГ, ПРИЗМА, УГПТЕСОМ), предназначенная для моделирования процесса 
распространения загрязняющих веществ, взаимодействия гидробионтов, преимущественно ориентирована на 
однопроцессорные вычислительные системы, представленные, в основном, персональными компьютерами. В 
таких системах масштабируются на параллельные системы только единичные составляющие модули этих 
систем (например, ЕСОЗ1М и МАОЗ1Р). На практике в систему уравнений Навье-Стокса и неразрывности 
добавляется малое слагаемое, пропорциональное второй производной по времени от функции плотности. 
Данный подход позволяет сгладить нефизические флуктуации плотности массы, импульса и полной энергии, 
переносимые по пространственной сетке быстрее скорости звука. 

Материалы и методы. Успешность разработки математической модели гидрофизических и 
гидробиологических процессов зависит от наличия и проработанности тестовых примеров и задач для изучения 
устойчиво наблюдаемых в морях явлений, таких как вертикальное перемешивание и перераспределение 
солёности и кислорода, галоклин и теромоклин. Для изучения этих явлений в работе используется модель 
гидродинамики, учитывающая баланс массовых сил и трансграничных потоков [1, 2]: 


Р(,рь) в, 
и 1Ур+Ь. 
ро |+ | - 4» (ть) + (,рь)+(У.п)-+ра , 
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где ру — плотность потока, 0/0 — скорость изменения плотности кинетической энергии, 
08/0: — скорость изменения плотности внутренней энергии, Т — тензор напряжений Т1и=-рп, Б — 
массовая сила, й = й(1,х) — плотность потока тепла, 4 — удельный приток тепла за счёт излучения. Из-за 


того, что эти явления обычно описываются вертикальным распределением параметров по глубине, то 
целесообразно получить упрощённую модель, допускающую оперативную идентификацию ненаблюдаемых 
параметров. Выделим цилиндрическую область У с основаниями на дне и поверхности воды площадью 
поперечного сечения 5. Спроецируем скорости, потоки и силы на вертикальное направление в предположении, 
что частные производные по х, у от параметров на горизонтальное направление равны 0 за пределами цилиндра, 
те. предположим горизонтальную однородность параметров водной среды. Запишем систему (1) В 
консервативной форме так, чтобы она позволила определить плотность массы ( Р ), механический импульс ( ру ) 


и плотность полной энергии (Р(К +=) ), К=у?/2. В цилиндре выделим бесконечно малый объём и 


предположим, что на каждый такой объём, составляющий У, действуют: сила реакции опоры (дна водобма), 
равная гидростатическому давлению, и сила, аналогичная силе трения, вызываемая вязкостью жидкости и 
переносом импульса, не равная нулю при вертикальных движениях жидкости. 

Если пренебрегать горизонтальными перемещениями жидкости и предположить, что существенны только 
вертикальные движения, и тот факт, что плотность водной среды существенно зависит от солёности, принять 
изложенные ранее упрощения и соглашения, то уравнения гидромеханики [1-3] в компактной форме 
представимы системой дифференциальных уравнений в частных производных [2]: 


2(ру) (ру?) __@р 
О Ох Ох 


д д 9(ру) № вот 
5: [Р(в+у" 2) +5 [+(р(в+т?/2)) |= о 2 р &(Т-Т.,),Т = Хе) ==/с, (1) 
& 0(у5) "| э 
+ = и— |, 
0 Ох Ох\о дх 


где р = Г(р,Т,5), == Г(р,у,Г) — эмпирическое уравнение состояния морской воды и уравнение, замыкающее 
систему по внутренней энергии соответственно; р — полное гидродинамическое давление; р — локальная 
плотность водной среды; у — проекция вектор-функции скорости на вертикаль (ось 7 направлена вверх от дна 
к поверхности); & —Й объемная плотность внутренней энергии; р — давление газа, заключённого в 
элементарном объеме между соседними слоями; РЕ — объемная плотность обобщённой силы (суммы сил), 
приложенных к элементарным объёмам жидкости, помимо давления; 5 = На) — концентрация соли; 5 — 
площадь поперечного сечения цилиндрической выделенной области, в которой предполагается наиболее 


интенсивное протекание процесса апвеллинга и транспорта соли; Т — абсолютная температура воды; Т. 


епу 
температура воды внешней по отношению к выделенному объёму; К — теплопроводность воды; 

$ =9,8м/ с? — ускорение силы тяжести; ци — коэффициент, характеризующий интенсивность переноса 
импульса из-за вязкости. 


Граничные условия, характеризующие свойство непротекания жидкости через породу, составляющую дно 
водного объекта, для системы (1) могут быть записаны в виде равенств: 


д ОТ д 
(уп) =0, 2 =0, —=0, ^^ =0. (2) 
Где п — вектор нормали, направленный внутрь расчётной области. 


Международный стандарт уравнения состояния морской воды [10] определяет плотность р морской воды 
как функцию солености 5, температуры Ти гидростатического давления р, которая имеет вид: 
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(3) 


где К(з, Т‚р) — средний модуль упругости; цифра 0 соответствует одной стандартной атмосфере (101 325 Па). 


При численном моделировании на основе (1)-(3), также востребованной как (3), установим 
алгебраическую связь давления с плотностью, температурой и солёностью. Квадратное уравнение имеет два 
корня: 


_ =Ур+(-А)-р(5.Т,р)+ А: рь(5,Т,р) 


Р!2 р, ь (4) 
р, = 2В-(р(5,Т, Р)-Ро(5,Т, р), 
р = (—4ВК,(5,Т, р) + А? -2А+1-р(5,Т, р)? + (8В-Ку(5,Т,р)-2А?-+2А). ру(5,Т,р): (5) 
р(5,Т, р) + (А? —4В-К.(5,Т,р))- ро($,Т, р)?, 
где А, В, № (57, р), К (5, Т, р) — переменные параметры, связь которых с температурой, солёностью и 
плотностью определена стандартной моделью морской воды; р(5, Т, р) =р, где Т — температура воды в 


градусах Цельсия. Скорость звука в солёной воде может быть выражена формулой: с(5.Т, р) =_/бр/бр. 


Предположим, что импульс, солёность и теплота объёма воды, ограниченного цилиндрической 
поверхностью, меняются только при изменении параметров на границе области (граничные условия Дирихле), 
в частности, изменении солёности на поверхности воды. Поведение механических параметров — плотности, 
импульса — задано условиями Неймана и внешнему воздействию недоступно. Независимо от этих допущений, 
внутри границы водной среды система уравнений (1) может быть записана в векторной форме модели 
переноса-реакции с источниковой добавкой [11, 5]: 

64 + ей =Т, 


0 дх о 


где Т= (0 Е/5+5х Ру/5+Е02Т/0х? -К(Т-Т. } — вектор, характеризующий взаимодействие потока с 


епу 


окружающей его жидкостью и планетой; а=(а ру Е) —Щ_ вектор консервативных переменных состояния; 


из <“ 
Г = (ру р-ру? у(Е+ р)) — вектор потоков, выполняющих роль обратной связи и замыкающих уравнение 


баланса. Из компонент вектора Г можно вынести общий множитель }Ё=у:4. 

На практике для решения уравнения переноса вида (6) хорошо себя зарекомендовала разностная схема, 
оператор послойного перехода которой получается линейным комбинированием аналогичных операторов 
перехода схем «кабаре» и «крест» [12—15]. Учет связи потоков массы и изменения плотности оказался более 
эффективным по регуляризованн0ой по Б. Н Четверушкину — квазигидродинамичекой — системе, 
аппроксимированной по схеме ВВЦП (разности со сдвигом вперёд и центральные по пространству). 

Дискретизацию модели сплошной среды будем проводить интегро-интерполяционным методом на 


равномерной сетке 5$ =5,х5, 5, ={х, =Ш,1=0-+и,п.й=[}, 5, ={И =, ]=0+т,т-т=Т}, где В — шаг 
сетки по вертикали, { — индекс узла (контрольного/конечного объёма в терминах метода Годунова) при 
нумерации по пространству, т — шаг сетки по времени, ] — номер временного слоя. Опишем конечно- 


разностную аппроксимацию модели (1). Разностные схемы, используемые для аппроксимации уравнений 
баланса первого порядка, дают сравнительно приемлемые результаты лишь при очень малом шаге сетки [3], что 
приводит к активному потреблению ресурсов вычислительных устройств и разработчиков, создающих 
алгоритмы и программы. 

Проблема постановки граничных условий и их согласованного задания при знакопеременной скорости 
решалась методом заполненности контрольных ячеек [7]. Учет максимальной заполненности и выражение её 
функциональной зависимостью от номера узла позволяет повысить точность аппроксимации граничных условий. 

Работа по теоретическому и экспериментальному подбору методов и схем аппроксимации привела к выбору 
методов расщепления [13] и регуляризации. Уравнения неразрывности  регуляризируется по 
Б.Н. Четверушкину. Данное уравнение подвергается разностной аппроксимации по схеме ВВЦИП. Уравнения 
баланса (переноса) импульса, солёности и полной энергии заменяются явными уравнениями, полученными 
линейным комбинированием различных аппроксимаций оператора переноса («чехарда» и «кабаре») [12, 15]. В 
соответствии с методом расщепления, перенос импульса (и скорости), солёности, полной энергии 
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аппроксимирован на дробном шаге. Изменение скорости движения воды, определяющей интенсивность 
переноса массы вещества, вводится на втором дробном шаге решением волнового уравнения. Обозначим 


р=р,, р=р"1°, р=р"Н,у=у, У=у"°, у=ум, 0<о<1. В расщеплённом виде [16], после введения в 


уравнение неразрывности регуляризатора и полудискретизации на сетке 5. аппроксимации частной 


производной по времени схемой «разности вперёд по времени», первые два уравнения системы (1) могут быть 
записаны в виде: 


бу-ру 9[ д» 
РА, (рь), =-ра$ (== |, 7 
о (ру), =-р85 (х) мрт (7) 
1 р, —р, а _ Р"н —2Р» + Р"" ь _ 
=“ +(69), |= Р5 Е ‚р=Л(рТ),с=с(5Т,р), (8) 
РР р. (9) 
т 


Уравнения (7) и (9) представляют собой дискретный аналог модели переноса и изменения импульса — 
второе уравнение системы (1). Уравнение (8) прогнозирует изменение поля давления при учёте неразрывности 
потока и информации о межузловых (межъячеечных) и граничных потоках. При допущении постоянности по 
времени на каждом шаге скорости звука уравнение (8) может аппроксимироваться неявной разностной схемой 
по времени и решаться прогонкой. Именно это направление в численном моделировании и было выбрано для 
исследования. 

Уравнение (7) и четвёртое уравнение (1) аппроксимируем на эквидистантной сетке явной разностной 
схемой [13]. Слагаемые (8), определяющие волновые свойства при распространении импульса (правая часть 
уравнения), аппроксимируем центральными разностями по пространству: 


1 —1 1 1 1 1 
Код рен — 2ре + р’ х Р-Р" у Р-Р (10) 
с? д ы 2 2. 2 ь 
где К, —_ заполненность ячейки, находящейся справа от узла с индексом а К, —_ заполненность ячейки, 
находящейся слева от узла с индексом 7 - К —_ степень заполненности контрольной области Х: 1/2 < х < Х; 1/2 , 


находящейся в окрестности узла с индексом 1 (К, ‚ К, ‚ К, — переменные величины, зависящие от индекса 1) [7]. 


Оставшиеся два слагаемых (8) определяют связь скорости изменения плотности с её потоком. В (8) входит 
отношение скорости изменения плотности к отрезку времени т, выражаемой разностью со сдвигом вперёд: 


+1 _ 
р-р 


0,1 2 7 


[* 


поток массы переносится в правую часть (8) и аппроксимируется центральными разностями: 


19(58) | бор -Ф),, (ри) - Фа), 


0,1 х ах _ Вы Ъ.х К. т (11) 
где К, ,К — степень заполненности областей, находящихся в окрестности ячеек с номером 1[7]; К, 
характеризует заполненность области [х1,х..|; К — [хьхи|, № — [хх]. 
Уравнения баланса механического импульса и давления (9) аппроксимируем ВВЦП: 
к Роу Р-Р", РР 2) 


о в а 
: т ‚ 26 26 
Представим задачу разрешения разностной аппроксимации уравнения (8) в форме задачи матричной 
прогонки с переменным во времени вектором в правой части Ах=ЁР: 
Ах, +Сх+Вхи=Е. 


Аппроксимация уравнения (8), определённая на трёхточечным разностном шаблоне, в форме линейной 


системы уравнений Ар”+! = {(р", р"-", (27) .) ‚ разрешаемой относительно давления (р), имеет вид: 
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АООТ 


1 р 1 1 0 ых а не 0 
2212 И? Ь? 
1 1 4 2 1 р 
р? с2т2 |2 „62 Р4 
: . в : рии 
ы - -. ы - и в (13) 
6? с : й | , в! 
0 1 1 2 1 рим 
62 с242 12 Ь2 
1 1 
0 + 


- (= (2-7) == (р-р!) \ о) + 
А А м т 
+ (Во -рь)/2 гос Р-р, = (Ву-ры)/2) + 
+1 (55) гк зе = а, 
2тИ 1 7+1 1-1 М-1 


Подобная ленточная матрица получается в ходе аналогичной аппроксимации двух других слагаемых (8). 

В модели могут быть учтены динамические изменения потока на границах и других параметров, таких как 
температура, солёность, кислородосодержание. Поэтому моделирование движения воды нужно выполнять в 
цикле по времени послойно с удержанием оперативной информации о параметрах как минимум на двух 
соседних по времени слоях решения сеточного уравнения. Алгоритм вычисления гидродинамических 
параметров на двухиндексной сетке по пространству и времени включает в себя: 

— построение прогноза изменения импульса по первому уравнению (7); 

— приближённое вычисление функции пространственного распределения давления как функции плотности и 
температуры с предыдущего временного слоя; 

— оценку изменения плотности по второму уравнению (8); 

—вычисление градиента давления по новым значениям плотности и температуры, корректировка 
распределения импульса по третьему уравнению (9); 

— нахождение нового распределения полной энергии и температуры по аппроксимированному третьему 
уравнению системы (1). 

При алгоритмизации метода решения системы расщеплённых уравнений введены: 

— двоично-числовые маски, предопределяющие переключение шаблона разностной схемы с изменением 
знака скорости; 

— переменные сдвиги индексов соседних узлов, позволяющие записать решение уравнения (12) для 
граничных и внутренних узлов единой системой вычислительных операций. 

Такие переменные сдвиги индексов используются при вычислениях аппроксимаций градиентов как в 
объёме модели сплошной среды, так и на границах со вторым порядком точности по дискретным аналогам 
уравнений (7)—(9), аппроксимированных по ВВЦП и по линейной комбинации схем «чехарда» и «кабаре» 
(САВАВЕТ) [12—14]: 


47 -а к —а"'! + Ч + а 


п+1 — (п 4 п — п-1 
Е - [2 - му т - Е 0,у 20, 
а 4" 4 9" а 4" > 9" 97 о И В Ч! (14) 
1 7 3 +1 р ул 141 7 4 1 1 у! 141 1-Е — 0, ул < 0, 
т 3 2 й 3* ЗА 


где характеристика 4 для переноса соли (а = 5 ) и импульса д=ри. 


Параметр ш представляет собой «переключатель» потока при изменении знака скорости в разностных 
аппроксимациях градиентов давления и потока массы: 


0, у< 0, 


ттт = 50 


При использовании такого переключателя аппроксимация градиента при решении задачи переноса для 
импульса и солёности может быть записана формулой: 


04 х 2 И: — Ч и 4—4; +1 Ч Ч , 


Ох 3 й й 3 2й 


х=х, 


(15) 


Информатика, вычислительная техника и управление 


219 


Бир://уезиК-аопзв.га 


220 


Адуапсей Епзтеетпия Везеагсй (Ко51оу-оп-Ооп). 2023;23(2):212—224. е155М 2687—1653 


Применение метода заполненности ячеек позволяет корректно аппроксимировать граничные условия при 
размещении параметров состояния воды в памяти компьютера массивом числовых значений, а переменные 
сдвиги индексов позволяют уменьшить объём символьной записи подпрограммы выполнения послойного 
итерационного изменения переменных состояния, отнесённых к узлам разностной схемы [17]. 

Матричное уравнение с ленточной трёхдиагональной матрицей А (13) решается прогонкой [18]. Для 
начальной верификации программной реализации использовалась упрощённая модель состояния солёной воды 
в виде уравнения П. С. Линейкина. 

Результаты исследования. Моделирование выполнялось на основе программного комплекса, написанного 
на языке МаНаб. Отладка проводилась с использованием интерпретатора СМО Осауе. Эксплуатация также 
предполагает наличие библиотек этой системы. Программный комплекс состоит из 16-ти функциональных 
модулей. Выбор этого интерпретатора и соответствующего языка обусловлен возможностью записать и 
проверить программу, оперирующую с массивом переменных состояния, являющегося проекцией искомых 
функций р(х,), у(х,1), &(х,0) на пространственную сетку. На этом этапе исследования авторы 


абстрагировались от специфики выполнения поэлементных операций над массивами действительных чисел. 

Программная система состоит из взаимосвязанных подпрограмм: 

— выполнение одного шага по уравнению переноса по формулам (10); 

— расчет плотности морской воды (ипезсо_лгз); 

— расчет давления морской воды, находящейся в поле силы тяжести, при заданной температуре, солёности и 
плотности по эмпирическим уравнениям состояния (оТ$2Р); 

— вычисление скорости звука в зависимости от температуры, солености и давления по стандарту и 
уточнённой формуле «ЮНЕСКО» (зрее4_оЁ_зоипа); 

— оценка величины силы вязкости жидкости при движении её объёмов под действием сил давления во всех 
точках пространственной стеки (ЕогсеОЁЕйсйоп); 

— циклическое варьирование переменных состояния и времени при вертикальном движении солёной воды 
(адиа_ргосез$); 

— начальная установка константных величин, характеризующих жидкость, начальные и граничные условия 
её глобального гидрофизического равновесия (бат, зе _рагатеегз); 

— формирование ленточной матрицы, аппроксимирующей уравнение, содержащее давление (Ёапс5); 

— решение матричного уравнения методом прогонки (гип_5\еер_зНиШе); 

— преобразование температура <> полная энергия (ТЕготЕ, ЕЕготТ); 

— расчет баланса полной энергии (То{а!Ро\ег); 

— решение задачи диффузии-конвекции-реакции, в том числе при аппроксимации комбинации разностных 
схем «чехарда» и «кабаре» (АБК _зоуег); 

— оценка величин градиентов и производных по времени по центральным и направленным разностям с 
учётом изменения знака скорости и шаблона разностной схемы (41123). 

Построенная модель использована в тестовом запуске для оценки изменения плотности при повышении 
солёности поверхности водной среды на 1% на 30-Й секунде от начала моделирования и максимальной 
глубине 1км. Переходные процессы по плотности во множестве отстоящих друг от друга сечений 
(функция р(Е, х), 5, ), вызванных мгновенным изменением концентрации соли, конечны во времени (рис. 1), а 


фронт импульса плотности сильно ослабляется при движении в пространстве. 
р, кг/м 


120 130 


1010 


Рис. 1. График зависимости плотности от времени при резком изменении солёности на поверхности воды 
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Обсуждение и заключение. Представленная пространственно-распределённая модель пока не позволяет 
прогнозировать устойчивое перераспределение плотности и изменение градиента солёности при численных 
расчетах, потому что в ней не учтена плавучесть менее солёной воды и еб изменение при осолонении верхних 
слоёв. Математическое и программное обеспечение, значительно упрощающее моделирование процессов, 
приводящих к наблюдаемым эффектам галоклина, хемоклина и пикоклина, проверено и отлажено на множестве 
тестовых задач (транспорт импульса, соли, распространение Волны давления). Для высокоточного 
моделирования наблюдаемых физико-химических явлений в морях понадобится решать дополнительные 
задачи идентификации параметров модели, взяв в качестве исходной информации данные наблюдений и 
дистанционного зондирования Земли. 
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Н.В. Кудинов — дополнение метода контрольных ячеек алгоритмической моделью масок и индексов; 
программная реализация решения гидродинамических уравнений: алгоритм обновления значений в 
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